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Harvey's S values for sylvine at 295°K agree qualita- 
tively with the theoretical S values calculated from James 
and Brindley’s quantum-mechanical f and Ej; values. 
Quantitatively, the difference between experiment and 
theory is too great for comfort. Accordingly, we have made 
a more exact and careful determination of the experimental 
S values and have found a remarkable agreement between 
experiment and theory. This shows the excellence of the 
present theory of diffuse scattering and the accuracy of the 


f and E;; values which James and Brindley have calculated 
by approximate methods. Our results are not in accord 
with the idea that the probability distribution is the 
same for all electrons. Also experiment has shown that 
turning the crystal plate in its own plane has no effect on 
the diffuse scattering in a given direction. This indicates 
that the atoms of sylvine possess spherical symmetry as 
far as the scattering of x-rays is concerned. Spherical sym- 
metry is assumed in the theory. 


1. INTRODUCTION 


CCORDING to the most recent develop- 

ment of the theory, the diffuse scattering 

of x-rays by crystals consisting of two kinds of 
atoms (such as KCl) is given by 


S=5S,+5:/(1+a vers ¢)', (1) 


where 
S,= (f2+f2— Fi — (2) 
and 
S:=1 
(Ej7)1+ | Eal P-L | 


(3 


Z:+2Z2 


The various symbols are defined in the paper by 
Harvey, Williams and Jauncey,' the numerical 
subscripts referring to the two kinds of atoms. 
It has become the custom to use another symbol, 


* The senior author was aided in part by a grant from 
the Rockefeller Foundation to Washington University for 
research in science. 

1G. G. Harvey, P. S. Williams and G. E. M. Jauncey, 
Phys. Rev. 46, 365 (1934). 


Scinss, Which is defined by 
Scisss (4) 


On the basis of wave mechanics, James and 
Brindley” have calculated values of (their fo) 
and of f for the K* and Cl~ ions in sylvine. 
Values of F are obtained by multiplying the 
values of f by the Debye-Waller’ temperature 
factors, using 230°K for the characteristic 
temperature of KCI! and assuming the existence 
of zero point energy. The theoretical values of 
Sand Seisss can therefore be calculated. 

Harvey has obtained experimental S.:... values 
for sylvine at room temperature, first for an 
average wave-length of a band of x-rays* and 
second for monochromatic x-rays.* In both cases 
he obtained a maximum value of Szisss of 2.95 
at about (sin ¢/2)/A =0.365. Although the wave- 


*R. W. James and G. W. Brindley, Phil. Mag. 12, 81 
(1931). 

*P. Debye, Ann. d. Physik 43, 49 (1914); I. Waller, 
Zeits. f. Physik 17, 248 (1923). 

*G. G. Harvey, Phys. Rev. 38, 593 (1931). 

5G. G. Harvey, Phys. Rev. 43, 591 (1933). 
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mechanical Syisss curve has the same general 
shape as the experimental S.ias, curve, there is a 
distinct difference between the absolute values 
of the ordinates of the two curves as shown by 
the fact that the maximum theoretical Soeiss. 
value is 2.46 at (sin ¢/2)/A=0.59. The difference 
between the experimental and theoretical values 
of x(=(sin¢/2)/A) for the maximum Soeiss 
values might be explained as due to error in 
drawing the experimental curves since the points 
are somewhat scattered, but a difference of 
about 0.5 in the two maximum values of Seiass 
is too great for comfort. We have therefore 
made a careful investigation of the diffuse 
scattering from sylvine in order to find out 
whether the above discrepancy between theory 
and experiment can be reduced or removed. 

A second object of this research was to find 
out whether the diffuse scattering in a given 
direction depends on the orientation of the 
crystal with respect to the plane of scattering. 


2. EXPERIMENTAL METHOD 


Because we desired greater accuracy in our 
readings we used a band of wave-lengths from 
the continuous spectrum of x-rays from a 
tungsten target tube operated at 51.5 kv peak 
rather than monochromatic rays as in the 
experiments of Harvey’ and Williams.* The 
experimental method has been described in 
previous papers.’ Referring to the paper by 
Jauncey and May,® let us introduce 


cos @—cos (¢—@) 


where ¢ is the angle of scattering, @ is the angle 
between the normal to the crystal slab and the 
primary beam, and yu and ¢ are respectively the 
absorption coefficient and thickness of the 
crystal. Then, if D is the electrometer deflection 
corresponding to a setting @ and @¢ is kept 
constant, the graph of D/v against @ is a hori- 
zontal straight line excepting where Laue spots 
occur as shown in Fig. 1. As pointed out by 
. ja 43, 
. a . a 
(1924). 


» (5) 


v= 
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5° 10° 6 15° 20° 25° 


Fic. 1. Graph of D/v against @ for ¢=30°. The hori- 
zontal portion of the curve is important in diffuse scat- 
tering. Each experimental point of this portion is the 
average of numerous readings. The Laue spots occur at 
+(¢/2—6)=0°, 11.3° and 14°. 


Jauncey and Williams’ it is essential to separate 
the Laue spots from the diffuse scattering. 
Previously we have plotted D against @ and this 
gives rise to an approximate straight line with a 
small slope. Fig. 1 shows how well the experi- 
mental points fall on a horizontal line. For the 
case where a principal axis of the crystal is 
perpendicular to the plane of scattering, the 
values of @ at which Laue spots occur are given by 


(6) 
where ) and & are integers and either but not 


both may be zero. Eq. (6) holds so long as the 
wave-length 


= (2d sin 6)/(h?+k*)! 


6= ¢/2+tan" h/k, 


(7) 


is contained in the continuous spectrum. The 
positions of the Laue spots according to Eqs. 
(6) and (7) are indicated by the vertical lines in 
Fig. 1. 

Because v=1 when @=¢/2, the average value 
of D/v for points not at or near Laue spots is 
the value assumed for D, in the right side of 
Eq. (17) of the paper by Jauncey and Pennell.® 
The same ionization chamber and thickness of 
paraffin were used as are described in this paper 
and so the numerical value 0.314 included in 
the formula for S.x, is the same as that used by 
Jauncey and Pennell. For clarity it might be 
well to point out that 0.314 stands for 


Z'p't’ 1 Koo°7 90° 
x x 
W’ cos45° Ko(i+a)? 


(8) 
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The value of pt for the crystal is 0.325 g/cm’. 

Jauncey and Pennell determined the intensity 
distribution in the spectrum by means of their 
Eq. (7a). This, however, assumes that the 
Lorentz factor'® is proportional to \*/sin 28. 
This is the factor for monochromatic rays when 
the method of integrated reflection is used. 
However, the ‘‘white’’ Lorentz factor is pro- 
portional to */sin? @ and would seem to apply 
to our determination of the spectral distribution. 
In addition to the Lorentz factor there is some 
uncertainty concerning the proper value of yu 
to be used since » should be increased by the 
extinction coefficient. This coefficient probably 
becomes a greater proportion of the ordinary 
absorption coefficient at shorter wave-lengths. 
We therefore determined three curves—the first 
curve was the graph of the actual electrometer 
reading against \, while the second and third 
curves were the graphs of the readings corrected 
for the. “‘monochromatic” and “white” Lorentz 
factors, respectively, both the second and third 
curves being also corrected for the ordinary 
absorption coefficient. By experiment we deter- 
mined that 0.75 g/cm? of aluminum was neces- 
sary to reduce the primary rays to half value 
when the rays had already passed through the 
paraffin (p’t/=0.493 g/cm?), the sylvine crystal 
and 0.170 g/cm? of aluminum which was perma- 
nently placed in the primary beam during the 
scattering experiments. From each of the above 
three curves a modified curve for the spectral 
distribution of the rays getting through the added 
0.75 g/cm? of aluminum was calculated, using 
the values of » for various wave-lengths in 
aluminum as given by Compton’s tables."" Now 
comparing the area under each modified curve 
with that under the corresponding original curve 
we find that the fractions to which the primary 
rays are reduced are: First curve, 0.463; second 
curve, 0.433; third curve, 0.332. Obviously the 
first curve agrees best with the experimental 
value 0.5 and this is the curve, Fig. 2, we have 
used for the spectral distribution curve. 

At first we compared the scattering from KCl 
with that from paraffin and used Jauncey and 
Pennell’s formula for calculating S.x». However, 
in this formula it is assumed that the scattering 


See F. C. Blake, Rev. Mod. Phys. 5, 169 (1933). 


1 A. H. Compton, X-Rays and Electrons, p. 184. 


~— 


Fic. 2. Distribution-in-wave-length of the intensity in 
the continuous spectrum. Half-value absorption in alumi- 
num of this spectrum corresponds to an average wave- 
length of about 0.38A. 


from paraffin at 90° is entirely incoherent so 
that the Breit-Dirac” formula holds. However, 
from Fig. 2 it is seen that there is some intensity 
of \=0.7A present and at this wave-length some 
coherent radiation occurs at ¢=90°. Since this 
part of the spectral distribution curve is some- 
what uncertain,’* we measured the absolute 
value of the spatial scattering coefficient per 
unit solid angle at ¢=90° from paraffin. This 
necessitates the comparison of two intensities 
whose ratio is in the order of 1 to 10,000. First, 
we made the comparison by a photographic 
method. A photographic film was cut into 24 
rectangular pieces and these pieces were arranged 
in two piles of 12 pieces each. One pile was 
exposed to the primary x-rays for several seconds 
and the other to the scattered rays for several 
hours. Each piece was run through the micro- 


photometer'* and the photographic density 


2G. Breit, Phys. Rev. 27, 362 (1926); P. A. M. Dirac, 
Proc. Roy. Soc. All1, 405 (1926). 

13The minimum wave-length as given by the Duane- 
Hunt relation is 0.24A so that second order reflection from 
the rocksalt occurs at @ ter than that corresponding 
to first-order reflection of \=0.48A and so a correction in 
the curve had to be made. 

See W. D. Claus, Phys. Rev. 38, 604 (1931). 
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determined.'® The purpose of the pile of films 
was to determine the absorption coefficients of 
the primary and scattered rays in the films 
according to the method of Brindley and Spiers.'* 
Also the area of the spot produced by the 
primary rays was determined. The result can 
be expressed by saying that the numerical factor 
0.314 in Jauncey and Pennell’s Eq. (17) should 
be replaced by 0.352. Second, we made the 
comparison by means of an ionization chamber 
and electrometer. We cut off a flash of primary 
rays of duration 1/60 sec. by means of a pendu- 
lum" as described by Jauncey and May.* The 
numerical factor is now found to be 0.358. In a 
third method we compared the intensity of the 
rays diffusely scattered from the crystal at 
¢@ =40° with that of the primary rays and found 
the factor to be 0.355. We therefore replaced 
the factor 0.314 in Jauncey and Pennell’s Eq. 
(17) by the average value 0.355. 


3. COMPARISON WITH WAVE-MECHANICAL 
VALUES 


Values of the theoretical quantity defined by 
the left side of Jauncey and Pennell’s Eq. (17) 
can be calculated from James and Brindley’s* 
wave-mechanical values of f and E;; together 
with the values of =| E;.|* given by Waller and 
Hartree'* and Harvey, Williams and Jauncey.' 
The values of S.x, and Sy, are given in Table I. 
The results are also shown by the solid curve of 
Fig. 3. It is seen how well the experimental 
values agree with the theoretical. The maximum 
Sexp of 2.31 corresponds very well with the 
maximum Sy, of 2.36 at ¢=25°. Assuming an 
average wave-length of 0.38A (obtained from 
half-value absorption) this corresponds to x 
=0.57. The reason for the slight difference 
between the values given here and in the 
introduction to this paper is that Sviass was 


See G. E. M. Jauncey and H. W. Richardson, 
J. O.S. A. 24, 125 (1934). 
(1983) W. Brindley and F. W. Spiers, Phil. Mag. 16, 686 

1% The approximate distance by which the pendulum 
should be pulled out before being released can be calculated 
but the exact distance can be found experimentally by 
petting the electrometer deflection against the distance. 

ere is a point of inflection at the correct distance. Also 
at the correct distance the readings repeat very well but 
as the distance departs from the correct value the readings 
scatter more and more. 

8]. Waller and D. R. Hartree, Proc. Roy. Soc. Al24, 
119 (1929). 
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TABLE I. Diffuse scattering from sylvine. 


Sexp Sn | @ Sexp Sth 
10° 1.50 1.65 | 37° 1.99 2.06 
15° 1.95 2.05 | 40° 1.90 1.95 
20° 2.23 2.28 | 45° 1.70 1.79 
25° 2.31 2.36 | 50° 1.56 1.62 
30° 2.27 2.28 | 60° 1.32 1.37 
33° 2.14 2.21 | 75° 1.14 1.12 
35° 2.01 2.13 | 90° 1.01 0.97 


used there and monochromatic rays were as- 
sumed instead of the band of wave-lengths 
actually used. 

The excellence of the agreement shown in 
Table I indicates the correctness of the present 
theory of diffuse scattering and also the high 
accuracy of James and Brindley’s* f and £;; 
values which were calculated by approximate 
methods. The broken curve of Fig. 3 is based 


d 


i 
2 6 7 


"3 
sin ¢/2 —> 


Fic. 3. Diffuse scattering of x-rays from sylvine. Black 
circles, experimental points; solid curve, electrons having 
different probability distributions; broken curve, electrons 
having equal probability distributions. 


upon the assumption that the probability 
distribution is the same for all electrons so that 
LE;7=f/Z, ie., =f in the notation used by 
Jauncey.'* Our results show that this assumption 
is incorrect. The omission of the negative term 
— | Ej, |*/Z from S; makes a maximum increase 
of 0.08 in Sy at ¢=29°. The inclusion of the 
negative term thus gives slightly better agree- 
ment with S,x, at this value of ¢. 


& G. E. M. Jauncey, Phys. Rev. 42, 453 (1932). 
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We wish to notice that the difference between 
Sexp and Sy, at ¢=35° is somewhat larger than 
our experimental error which we place at about 
0.05 in S.x, in this region. We made repeated 
trials in this region and always found S.x, 
somewhat less than S,,. A rather large difference 
also occurs at ¢=10° and 15°; but at these 
small angles the experimental difficulties increase 
and the difference is probably not greater than 
the experimental error. Considering all the values 
of ¢, there is a tendency for all the experimental 
points to fall below the theoretical curve. This 
is probably due to error in the determination of 
the numerical factor mentioned at the end of 
the last section. 


4. Errect or CRYSTAL ORIENTATION 


In all the work done hitherto in this laboratory 
x-rays have for convenience been scattered from 
a crystal with one axis perpendicular to the 
scattering plane. We have now tried the effect 
of turning the crystal plate in its own plane 
upon the S,x, values. When the crystal is 


oriented so that no principal axis is perpendicular 
to the scattering plane, Eqs. (6) and (7) for 
the positions of the Laue spots become more 
complicated and contain a third integer 1. The 
S.xp values for the axis making angles of 0°, 
22.5° and 45° with the perpendicular are shown 
in the second, third and fourth columns of 
Table II. It is seen that turning the crystal in 


Taste II. Effect of crystal orientation—Sexp values. 


0° 22.5° 45° 
15° 1.95 1,99 1.93 
25° 2.31 2.31 2.29 
35° 2.01 2.00 2.01 
50° 1.56 1.54 1.54 


its own plane has no noticeable effect. The 
theory of diffuse scattering has been developed 
on the assumption of spherical symmetry in the 
atoms.” Our results are in accord with this 
assumption in the case of sylvine. 


2°G. E. M. Jauncey, Phys. Rev. 37, 1193 (1931). 
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The Verdet Constant of Heavy Water 


Francis G. Stack, Department of Physics, Vanderbilt University 
(Received October 1, 1934) 


The Verdet constants of samples of water containing 0, 31.1, and 99.7 percents of heavy 
hydrogen atoms have been measured for the wave-lengths \5893 and \5460.7A at 20°C. The 
resuits are given in tabular form. From these data the Verdet constant of pure heavy water 
(H?*,O) is found to be 3.93 percent less than that of ordinary water for both wave-lengths. Non- 
linearity of the Verdet constant-concentration relation leads to the calculation of an apparent 
Verdet constant for H'H*O which is 0.85 percent less than the Verdet constant of ordinary 
water. Values of the relative molecular rotatory power and of the magnetic rotatory dispersions 


are also given. 


HE Verdet constants of ordinary and of 
heavy water have been measured for two 
wave-lengths of light. 

The measurements were made by use of the 
Schmidt and Haensch Faraday effect apparatus. 
The optical system of the instrument consists of 
a triple field polarizer with the Lippich half-shade 
adjustment and an analyzing Nicol whose 
position may be read to 0.01°. The magnetic 
field was produced by the attached electromagnet 


with an iron core, drilled parallel to the magnetic 
field. The mean axial field strength was deter- 
mined by calibration against a standard solenoid 
by the ballistic galvanometer deflection method. 
The test coil used between the magnet pole 
pieces in calibration was of approximately the 
same dimensions as the cell used to contain the 
liquids whose Verdet constants were determined. 
The current for the electromagnet was supplied 
from storage batteries, and was determined by 
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the measurement of the potential drop across a 
standard resistance by means of a Wolff po- 
tentiometer. Current measurements were made 
before and after each observation and the 
average of these values used to determine, from 
the magnet calibration curve, the actual field 
strength at the time of observation. Two different 
field strengths were used in each determination. 
The light sources used were the sodium and 
mercury arcs with appropriate filters to transmit 
only the lines 45893 and \5460.7A. 

The liquids under observation were sealed in 
water-jacketed Pyrex cells approximately 5 cm 
long with thin cover glass ends cemented on with 
picein. The rotations due to these glass ends were 
measured and corrected for, the corrections 
being of the order of one percent of the total 
rotations. The temperature of each specimen 
was held constant to within 1.0°C by a constant 
flow of water through the jacket. The tempera- 
ture of the liquid is taken as the average of 
thermometer readings at the inlet and at the 
outlet from the jacket. These readings differed 
at most by a few tenths of a degree. The two 
cells, made by Eck and Krebs of New York, 
were the same length to within one part in 4000. 
Measurements of the Verdet constant of ordinary 
water made in the two different cells agreed 
within the estimated precision. 


Samples of heavy water of two different 
concentrations and one sample of ordinary water 
were kindly loaned by Dr. H. C. Urey and his 
assistants. These samples were freshly distilled 
and were sealed in the cells before delivery. 
The concentrations of these samples were also 
determined by Dr. Urey’s staff. 

The results of the measurements are given in 
Table I. Column 1 gives the concentration of the 
solutions in mol fractions of H*,O0. Column 2 
gives the specific gravity as computed from the 
known concentrations by the equation of Luten,! 
and using the value D**,; = 1.1079 given by Taylor 
and Selwood? for 100 percent H*,O. The indices 
of refraction in column 3 are taken from the 
data of Luten' by interpolation from his tables. 
Columns 4, 5 and 6, respectively, give the wave- 
length, temperature and field strength at which 
the measurements were taken. Each value in 
columns 5, 6 and 7 is the average of at least 
two independent sets of observations. Columns 
7 and 8 give the Verdet constants as measured. 
The values in column 8 are to be taken as the 
final values and are the weighted means of the 
values shown in column 7, corrected* to 20°C. 
By extrapolation the values 0.012554+0.000030 


1D. B. Luten, Jr., Phys. Rev. 45, 161 (1934). 

?H. S. Taylor, and P. W. Selwood, J. Am. Chem. Soc. 
56, 998 (1934); H. S. Taylor, J. Frank. Inst. 218, 1 (1934). 

* Int. Crit. Tab. 6, 425 (1929). 


TABLE I. The Verdet constants of solutions of H*,O in 1,0. 


Concentration Specific Refractive Wave- Temp. of Field Verdet Verdet 
mol fraction gravity index length liquid strength constant constant 
(H*,0) (20°C) (°C) (gauss) (min. /g-cm) (20°C) 
0.000 1.0000 1.33449 5460.7 20.89 2205 0.015426 0.015397 
-000 1.0000 1.33449 5460.7 20.40 1953 -015333 +.000030 
-000 1.0000 1.33300 5893 19.80 2207 -013069 .013067 
.000 1.0000 1.33300 5893 19.80 1954 .013065 +.000030 
311 1.0333 1.33302 5460.7 20.00 2202 015293 .015274 
311 1.0333 1.33302 5460.7 20.00 1945 .015256 +.000030 
311 1.0333 1.33158 5893 20.00 2206 .012979 .012963 
311 1.0333 1.33158 5893 20.00 1952 .012948 +.000030 
.997 1.1076 1.32977 5460.7 20.23 2196 .014798 .014793 
.997 1.1076 1.32977 5460.7 19.90 1960 .014782 +.000030 
.997 1.1076 1.32845 5893 19.90 2191 .012583 .012556 
997 1.1076 1.32845 5893 19.90 1960 .012530 +.000030 
1,000 1.1079 1.32976 5460.7 20.00 .014791* 
1.000 1.1079 1.32844 5893 20.00 .012554* 


* These values obtained by extrapolation. 
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and 0.014791+0.000030 min./gauss-cm are ob- 
tained for the Verdet constants of 100 percent 
H?,0 at 20°C for the wave-lengths \5893 and 
45460.7 respectively. 


DIscussION OF RESULTS 


These measurements show the Verdet constant 
of heavy water (H*,O) to be 3.93 percent less 
than that for ordinary water (H',O) at each of 
the two wave-lengths investigated. Thus the 
rotatory power relative to H',O is independent 
of wave-length in this region. The indices of 
refraction for the two waters at \5893 differ by 
approximately 0.34 percent and the dispersions 
(dn/dd) between the wave-lengths \5893 and 
5461 differ by more than 11 percent.'! The 
failure of the linear relations between Verdet 
constant and concentration and between mag- 
netic rotatory dispersion and concentration is 
striking. This does not necessarily indicate a 
violation of the Schénrock rule of mixtures, 
however, for we must recall the possibility of 
H'H°O molecules in the solution. Assuming the 
equilibrium constant K = (H'H*O)*/(H*,O) 
=3.27 as found by Topley and Eyring* we may 
compute from the known ratio of the number 
of heavy to light hydrogen atoms the percent of 
each type of molecule in a solution. This calcu- 
lation results in the following concentrations for 
the 0.311 solution: 10.61 percent H*,O; 48.38 
percent H',0; 40.98 percent H'H*O. Assuming 
the linear relation between Verdet constant and 
concentration we may then calculate the ap- 
parent Verdet constant of a liquid consisting 
only of the H'H*O molecules. The values so 


obtained are 0.012955 and 0.015266 min./gauss- 


*B. Topley, and H. Eyring, J. Chem. Phys. 2, 217 
ang A. Farkas, and L. Farkas, J. Chem. Phys. 2, 468 


cm for the wave-lengths \5893 and 5460.7, 
respectively. These values are both 0.85 percent 
less than the corresponding values for ordinary 
water, indicating approximate independence of 
wave-length for the relative rotatory power as 
was seen above in the case of H*,O. The assump- 
tion of any other equilibrium constant (K =4, 
for example) does not lead to this result. This 
thus indicates a means of checking the equi- 
librium constant of such solutions by measure- 
ment of the Verdet constant at different wave- 
lengths. Further checks by measurements on 
several different concentrations are suggested. 

The relative molecular rotatory power (Q) of 
H?,0 with respect to H',O at 20°C for 45893 and 
for \5460.7 is 0.9640 assuming molecular weights 
of 20.027 and 18.014 for the respective molecules. 
Lack of data on the density of H'H®O prevents 
calculation of its molecular rotatory power. 

The magnetic rotatory dispersion in the wave- 
length region used may also be obtained from 
the data. Expressed in min./g-cm per A the 
values for H',O, H*,O, and H'H"O are, respec- 
tively, 5.39 10~*; 5.18 10~*; 5.35 

Although the value measured for the Verdet 
constant of ordinary water (A5893) agrees within 
the precision of measurements stated in the 
table with the value given in the International 
Critical Tables, this precision is stated for 
comparative purposes only and not for the 
absolute values. 

The experimental work was done in the Ernest 
Kempton Adams laboratory of the Columbia 
University Department of Physics. We wish to 
express our great appreciation to Professors G. 
B. Pegram and H. W. Webb for the use of this 
laboratory and equipment. We are of course also 
greatly indebted to Dr. H. C. Urey and his 
assistants for the loan of the specimens of 
heavy water. 
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Atomic Energy Relations. I* 


R. F. Bacnert anv S. Goupsmit, University of Michigan 
(Received August 10, 1934) 


A simple method for the calculation of approximate 
energies of atomic levels is presented in this paper. It is 
based on the derivation of linear relations which express 
the unknown energy in terms of observed energy values of 
the atom and its 1ons. It is shown that the degree of approx- 
imation increases with the amount of experimental data 
available for use in the calculation and also how the best 


formulas can be obtained for each case. Several tables are 
given containing formulas for configurations involving s 
and ? electrons. They are applied to the spectra of carbon, 
nitrogen and oxygen and the energy values so determined 
are compared with those known from observations. In an 
appendix the method of approximation is compared with 
the quantum mechanical perturbation method. 


I. INTRODUCTION 


HE energy states of a system consisting of 
several electrons in a central field can be 
found approximately by the use of the perturba- 
tion method in the solution of the wave equation. 
The problem of determining the energy states 
with an accuracy approaching that with which 
they are known experimentally is sufficiently 
complicated of solution, that at present it has 
been carried out for the two electron case only. 
For many electrons the problem has been solved 
in the first order by Slater' and calculations 
have been made for several elements with wave 
functions determined by the method of Hartree. 
These calculations, while remarkably good when 
the complexity of the problem is considered are 
still not accurate enough in general to be of much 
use in the prediction and the correlation of 
atomic energy states. For this purpose it is 
desirable to have an easier and more accurate 
way of obtaining the energy states. 

In this article we shall develop approximate 
relations between the energy states of an atom 
and its ions. An unknown energy level can then 
be calculated by expressing it in terms of 
observed energy values. This is much less 
difficult than the direct quantum mechanical 
calculation because the relations derived here 
are simple linear expressions. Furthermore, the 
relations possess the favorable property that the 


* This is the first of a series of articles on this subject. 
Later papers will contain relations in isoelectronic spectra, 
formulas for configurations involving other than s and p 
electrons, a special treatment of the cases in which electrons 
are “‘lacking”’ from a complete shell and other extensions of 
the method presented here. 

t At nt at Columbia University. 

1j rol Slater, Phys. Rev. 34, 1293 i929). 


degree of accuracy increases rapidly with the 
amount of experimental information used. In 
almost all important cases the results are suffi- 
ciently accurate to be of value for spectroscopic 
purposes whereas such accuracy can _ hardly 
ever be reached with other methods used at 
present for the calculation of energies. 

The method by which observed energies are 
used for the prediction of unknown levels is 
based on principles which are also very simple. 
In a one-electron atom or ion (that is one 
electron outside of a core of closed shells), the 
observed energy levels give direct information 
about the binding energy of that outer electron. 
If we add an electron the observed energies of 
the two-electron ion compared with those of the 
one-electron ion give information about the 
change in energy due to the interaction between 
the two electrons. From the observed energies 
in three- and more-electron spectra further 
details about the interaction can be obtained 
and it is this information which is used for the 
calculation of unknown energies. 


Il. THe MANy PARTICLE PROBLEM 


(a) Let us consider the one-particle problem, 
for example one electron outside a core of closed 
shells. We shall denote its quantum states by 
A, B, C, etc., and the energies to remove the 
particle by W(A), W(B), W(C), etc. We add a 
particle and consider next the configuration AB 
of the two particle system.? If there were no 


? For the present it is assumed that all degeneracies have 
been removed. It should further be noted that we assume 
that it is possible to ascribe each energy level to a definite 
configuration. This condition is fulfilled in all cases in which 
this method of approximation can be successfully applied. 
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interaction between the two particles the total 
energy to remove both would be 


‘o(AB)= W(A)+ W(B). (1) 
The energy actually observed, W(AB), differs 
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The “triple energy” represents the change in 
the total energy due to the influence of the 
third particle on the interaction of the other two. 

In the system of four particles we can set up 
three successive approximations W)(ABCD), 


s from this by an amount w(AB) which we shatt W,(ABCD), and W,(ABCD) which we shall call 
mn, call “pair energy’’ and which is defined by the “zeroth,” “first,” and “second approximation,” 
ed. relation respectively. W, uses only the one-electron 

W(AB)= W(A)+W(B)+w(AB). (2) energies, W, includes the “pair energies” and 


W, the “pair” and “triple energies.” The 


The “pair energy” represents the influence of &XPressions are 


the interaction of the two particles on the total 


W (ABCD) => W(a), (7) 
i were no interaction the tota => W(as)-25 Wa), (8) 


at Wa). (9) 


This will of course be a very rough approximation 


to the energy actually observed, W(ABC). By The difference between the last approximation 


rd using the values of the “‘pair energies’’ one is and the energy actually observed gives the 
le led to a next approximation “quadruple energy,” defined by 
W,(ABC)= W(A)+ W(B)+ W(C) W(ABCD)= > W(a)+ 
mn. 4 The energy of a system containing five 
of _ WilABC)= WAB)+W(AC)+W(BC) particles can be built up in much the same way 
‘e —W(A)—W(B)—W(C). (4) as those of systems containing fewer particles. 
e . 
en The notation will be simplified by writing 
ies W(as) to indicate the sum of two-particle the 
ler energies and W(a) for the sum of one-particle if P od 
ed energies, the sums to be extended over the states ae 
he of the configuration such as in Eq. (4). In this W,)(ABCDE)=>W(a), (11) 
notation Eq. (4) can be written W,(ABCDE) =. W(a8) W(a), (12) 
W,(ABC)= > W(as) — > Wa). 5 
)= W(a8) — W(a) (5) W2(ABCDE) = W(a8) 
The right sides of Eqs. (4) of (5) contain only 
observed energies. This expression is often a +32 W(a), (13) 
good approximation to the actual energy W;(ABCDE)=> W(apy) 
W(ABC). The difference between the actual rs " 
energy and this appronimation we hall call the +2 W(ab)—2W(a). (14) 
- “triple energy,” w(ABC). It is defined by the It is clear that the procedure which has been 
ve equation followed above can be continued and thus 
W(ABC)= W(A)+W(B)+W(C)+w(AB) extended to any desired number of particles. 
as Indeed the approximate energy including the 
me +w(AC)+w(BC)+w(ABC) interaction in groups one less than the total 
ute al number of particles may easily be written for 


ed. W(ABC)= > W(a8)—-X W(a)+w(ABC). the case of particles. 


(6) 
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W,-2(A1A2- :-A,)= * *Qn—1) 
— > 
+(—1)""D War). (15) 


The fact that in all‘important cases the ‘‘pair,”’ 
“triple,” “quadruple energies,” etc., form a 
rapidly decreasing series is the reason for the 
usefulness of this method of approximation. 

It has been tacitly assumed in the relations 
which have been developed above that it is 
possible to determine on which states of the 
(f-—1), (f—2), etc., particle systems, a particular 
state of the f particle system is built. It is 
necessary to know for each state into what 
state of the single particle problem it will be 
transformed if the interaction between the 
particles is removed. Although this is always 
uniquely determined theoretically, there are 
cases where it is not known with certainty. For 
atoms, this happens when levels of different 
configurations lie so near together that the large 
perturbations make it difficult to assign con- 
figurations to the individual levels. For the 
important low energy states of atoms, this 
difficulty occurs infrequently. 

A further difficulty which appears in the 
atomic problem when only the electrostatic 
interaction is considered, is the high degree of 
degeneracy. Due to this degeneracy, it seems at 
first impossible to determine for a multiplet of a 
configuration of equivalent electrons, on which 
multiplet of the ion it is built. This is purely a 
formal difficulty which can be overcome and it 
will be discussed later (section III) when the 
atomic problem is presented in greater detail. 

(b) It has been shown in the preceding section 
how the energy states of the many particle 
problem can be related to those of the cases of 
successively less particles down to and includ- 
ing the one particle problem. It is frequently 
the case, however, that some of the energies 
which are needed for the determination are not 
known. In such a case it will be shown that the 
energy of the given many-particle problem can 
still be found (somewhat less accurately, to 
be sure) in terms of those which are known. 
Suppose that one wishes to express the energy 
of a state of the five particle problem in 
terms of those of four, three and two particles, 
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without using the energies of the one particle 
problem. For this purpose we consider Eq. (9). 
In the five particle case there are four other 
equations similar to (9) which may be written 
for W:(ABCE), W:(ABDB), W:(ACDE) and 
W.(BCDE). With these five relations it is 
possible to eliminate the quantities W(A), W(B), 
W(C), W(D), and W(E£) from relation (14). In 
the result it is necessary to replace the approxi- 
mate energy W.(ABCD) by the exact energy 
W(ABCD), and similarly for the others, because 
the approximate energies involve a knowledge 
of the single particle energies. The result is 
given in (16), the energy of the five particle 
problem now being written W,'(ABCDE) due 
to the approximations that have been made. 


W2'(ABCDE) = 3/40 
—1/2)) W(aBy)+1/4 Was). (16) 


It is not necessary to go through this lengthy 
algebraic elimination to obtain (16). Since there 
are three coefficients which may be adjusted, we 
can fix them such that the Oth, 1st and 2nd 
approximations of both sides are identical. The 
number of approximations which can be balanced 
in this way depends entirely upon the number of 
coefficients on the right side of the equation and 
that may change according to the specific case 
considered. This observation leads to the fol- 
lowing practical method of determining the 
coefficients. Let us first rewrite (16) with 
unknown coefficients as in (17). 


W2'(ABCDE) =x 
Was). (17) 


If the Oth order approximations are to be the 
same on both sides of (17) then the coefficients 
must be so adjusted that a particular state, say 
A, occurs the same number of times on each 
side. In the first sum there are four terms 
containing A, in the second sum six and in the 
last four, and A occurs once on the left side. 
This gives the numbers needed for (18). 


(A) 1=4x-6y+4s. (18) 


Proceeding similarly with a pair of states, for 
example (AB) and next with’a group of three, 
(ABC), we get two more equations. 


or 
e, 
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(AB) 1=3x—3y+s, (19) 
(ABC) 1=2x—y. (20) 


The solution of Eqs. (18), (19) and (20) gives 
the coefficients of (16). 

The same method used above to determine 
the coefficients can be used to get an idea of the 
relative approximation of the formulas. Let us 
compare (13) and (16) which give W2(ABCDE) 
and W,’(ABCDE). 


W:(ABCDE) = W(aBy) — W(a8) 
+30W(a), (13) 
W2'(ABCDE) = 3/4 W(aBy8) 
— 1/200 Was). (16) 


Each of these relations is a second approximation, 
that is, the interaction in groups of three as 
well as the interaction in pairs has been taken 
into account. W2'(ABCDE), however, is a some- 
what better approximation. If we consider the 
interaction in groups of four, for example in a 
particular group of four (ABCD), it is seen that 
the interaction of this group of four occurs once 
in W(ABCDE). This interaction does not ap- 
pear at all in the expression given above for 
W.(ABCDE) while it occurs with a factor 3/4 in 
the expression for W,'(ABCDE), so the latter is 
a better approximation to the exact energy 
W(ABCDE) than W,(ABCDE). 

In the practical application of these relations 
to atoms, one seldom knows the energies of all 
ions down to the single electron case. We shall 
now show that this method can be used when 


the energies are only those of the electrons’ 


outside the core. Let us take a group of three 
electrons in states A, B and C outside a core of 
two electrons in states D and E. The symbol * 
will be used to indicate an energy relative to the 
core, thus W*(ABC) is defined by (21). 


W*(ABC)=W(ABCDE)—W(DE). (21) 

From (4) we can write a relation for W,*(ABC). 
W,*(ABC)= W*(AB)+W*(AC)+W*(BC) 

—W*(A)—W*(B)—W*(C). (22) 


In terms of total energies including the core (22) 


. may be written as (23). 


W,"(ABCDE) = W(ABDE) + W(ACDE) 
+ W(BCDE) — W(ADE) — W(BDE) 
—W(CDE)+W(DE). (23) 


Since the energies of the singles and pairs cancel 
in ‘this equation, it is a first approximation. In 
addition, all those groups of three which involve 
the core states D and E also cancel and the 
only one which is not considered is the triple 
w(ABC). This triple energy is just the difference 
between the exact energy and the first approxi- 
mation in the three electron problem, so it 
appears that, including terms of the second- 
approximation, Eq. (22) which expresses the 
energies of the three electrons relative to the 
core, neglects the same interaction which is 
neglected in the problem of three electrons 
without a core. There are, however, in the 
problem including the core, several quadruple 
energies (those not containing both D and E£) 
and of course the quintuple energy which are also 


neglected. Since these are in general considerably _ 


smaller than the triple energy which was neg- 
lected it is apparent that the energies may well 
be considered relative to the core, and ‘the 
resulting equation may be expected to hold 
nearly as well as the similar equation for the 
same particles without core. . 

(c) In order to get an idea of the accuracy to 
be expected from the relations presented above, 
the method of approximations will be compared 
with the usual perturbation theory. Table I 


TABLE I. Relation of the energy corrections in the recurrence method to those ; 


in the perturbation 


Energy 
correc - 
tions 


d- in- 
Pertur- Pairs | Triples Ete. 
ion 


energies 


Ist order 
2nd order 
3rd order 
4th order 
etc. 


shows the relation between the two methods. 
The columns in this table give the energy 
corrections as they have been applied in this 
paper, first the pair energy, next the triples, etc. 
The rows list successively the energy corrections 
given by the ordinary perturbation theory. The 
energy correction which arises from’ the source 


le 
In 
| 
ze 
5) | 
Ly 
ve 
id 
1e 
d 
of é 
id 
e 
h 


— 


indicated at the head of a row or column is 
split up into all the individual parts, designated 
by crosses in that row or column. Thus, there is 
a contribution to the total pair energy from each 
order of the perturbation starting with the first. 
It is, indeed, immediately apparent that if the 
relation considered is a first approximation or 
better and if the exact energies of the two 
particle problem are used in that relation, then 
all the higher orders of the pair interaction have 
been included. It appears from Table I that the 
triple energy enters for the first time in the 
second order of the perturbation theory. It also 
appears that if the pair energy is considered 
then the complete first order of the perturbation 
theory is included. Similarly if the triple energy 
is considered then the complete second order of 
the perturbation theory is included, etc. The 
proof of these statements will be found in 
Appendix I. It will appear that the first approxi- 
mation is expected to hold better than a first 
order calculation by the perturbation theory 
since the part of the second order which is 
included in the former is considerable. 


Illa. ABsoLuTE Atomic ENERGIES 


The recurrence method described in the pre- 
ceding section will next be applied to the problem 
of atomic energy states. The present treatment 
is limited to those states involving s and p 
electrons; the states involving d electrons will 
be treated separately, later. It is essential for 
the application of the recurrence method to be 
able to correlate states of successive stages of 
ionization. The great degeneracy of the problem 
offers a serious complication. To overcome this 
complication it is necessary to introduce quan- 
tum numbers which describe the system with all 
degeneracy removed, and second to correlate 
these non-degenerate states with the degenerate 
states which are found from experiment. This 
involves no fundamental difficulties and can be 
carried out by a procedure which has already 
been used several times for atomic energy 
problems and which originated with the well- 
known sum relations of Pauli. 

Let us consider the problem of an atom with 
electrostatic and spin-orbit interaction of the 
electrons neglected. If a magnetic field is intro- 
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duced the degeneracy is completely removed and 
each electron is in a particular quantum state 
which is characterized by the quantum numbers 
n, l, s, m,, and m, representing for each electron, 
respectively, the total quantum number, the 
orbital quantum number, the spin quantum 
number and the projections of the spin and 
orbital momenta on the magnetic field. If we 
now introduce the electrostatic interaction be- 
tween the electrons but still neglect the spin- 
orbit interaction (Russell-Saunders coupling), 
each state is characterized by a particular value 
of the resultant spin S, a particular value of the 
resultant orbital moment L and their magnetic 
projection Ms and M,. It is necessary to find 
how the states characterized by the quantum 
numbers n, 1, s, m, and m;, of the individual 
electrons are related with those denoted by S, 
L, Ms and M;. 

Since Ms and M, are the total projections 
of the spin and orbital momenta, we have the 
relations 


Ms= dim, (24) 


The sums are taken over all the electrons. In a 
few cases where there is only one level with a 
given Ms and M,, the sum relations (24) are 
sufficient to correlate these with the individual 
electron quantum states. In general, however, 
this correlation is somewhat more involved. 


Two electrons 

Let us first consider two equivalent p electrons. 
In Table II are gathered all the states of p* for 
which Ms+M,=0. It is known that this 
configuration gives rise to multiplets *P (S=1, 
L=1), 'D (S=0, L=2) and 'S (S=0, L=0). 
The first two columns give the various states of 
the individual electrons omitting the quantum 
numbers n, / and s. If one of these states has 
energy W(A) and W(B) for the individual 
electrons then the particular combination of 
multiplets found in the last column has the 
energy W(AB) as is indicated at the bottom of 
the table. The third column gives the total 
projections Ms and M,. The last column gives 
the multiplet or combination of multiplets which 
can be assigned to each individual non-degenerate 
state. It is clear that each of the first four rows 
are states belonging to *P since only for *P is 


| 
| 
| 
— 
| | 
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TABLe Il. Configuration p*. 


TaBie III. Configuration sp. 


there a projection Ms s=+1. Similarly the fifth 
state belongs to 'D since only for 'D is there a 
projection M, =2. The assignment of the combi- 
nation of multiplets to the states with Ms=0, 


M,=1 and those with Ms=0, M,=0 cannot ° 


be made in this way. To the first set belong *P 
and 'D and to the second *P, 'D and 'S. The 
coefficients in the linear combinations which are 
correlated with each state can be found by 
standard quantum mechanical methods, for 


. example by the famous & 7 method. They can 


also be found directly from the recurrence 
method, as is demonstrated in Appendix 2. 
Since we neglect the spin-orbit interaction, it is 
not necessary to distinguish between the various 
J-levels of each multiplet. In order to find the 
energies it is not necessary to know the coeffi- 
cients for the Ms=0, M,=1 group since *P and 
'D are determined elsewhere, but for the determi- 
nation of the coefficients in more complicated 
configurations directly from the recurrence 
method, it is necessary to know the proper linear 
combination for each individual state. It should 
be noted that the sum of the coefficients for each 
state, as well as the sum of the coefficients of a 
particular multiplet in a group with given Ms 
and M Ly is 1. 

Table III gives the configuration sp, showing 
the linear combination. of *P and 'P for each 
non-degenerate state with Ms+M,2=0. All 
those states with Ms=+1 are *P alone, while 
the two sets with Ms =0 are *P and 'P together. 
The coefficients in the linear combinations can 
be found in a way siniilar.to that used above for 
?* or they ¢an be written directly as 1/2, since 
one state with Ms=0 and M,=1 van be ob- 


™, mi ms m Ms ML Multiplets m miMsg ML Multiplets 
1/2 1 1/72 1 1 1/2 (1 i 1 
1/2 1 yz -1 1 0 1/2 172 1 0 ap 
1/2 —1 1 1/2 72 -17 1 -1 
—172 -1/2 O| 1 —1/2 —1/2 1 
72 1 | —1/2 1 2 1p 1/2 —1/2 1 
72 1 | 1/24P+1/2'D 
172 O| —1/2 1 1/23P+1/2'D 1/2 -172 0 1/2*P+1/2'P 
—1/2 172 1/2*P+1/2'P 
72 1] -1 0 | 1/23P+1/6'D+1/3'S 
172 O| -172 O 2/3 W(A) W(B) Ms ML W(AB) 
72 -1| -1/72 1 1/2 *P +1/6'D+1/3 1S 
W(A) W(B) Ms ML W(AB) . 
as tained from the other by changing the signs of 


the spins. The same is true for the two states 
with Mgs=0 and M,=0. It is not necessary to 
give a table for s* since this gives but one 
multiplet, 'S. With the help of these tables for 
the two electron systems, we can now build 
configurations of more than two electrons, 
restricting ourselves at present to s electrons and 
to equivalent p electrons.” 4 


Three electrons 

The first three electron system to be con- 
sidered is the configuration »* which gives the 
multiplets ‘S, 2D and *P. The first three columns 
of Table IV give the individual electron states 
A, B and C. The fifth column gives the linear 
combination of multiplets assigned to each state. 
The first state is uniquely *S and the second *D. 
The assignment of linear combinationg to the 
other states is discussed in detail in Appendix 2. 
The last column of Table IV gives the sum of 
the two electron states (AB)+(AC)+(BC) on 
which the three electron states are built. For 
example, if we consider the first state we find 
that each of the three possible two electron 
states belong to p? *P in Table II. 

The configuration sp* gives the multiplets ‘P, 
°P,*D and *S. The first three columns of Table V 
give the states of the individual eleetrons. The 
fifth column gives the linear combination of 
multiplets corresponding to each state. Each of 
the first three states is uniquely ‘P and each of 
the next two is uniquely *D. The proper linear 
combination to be assigned to the other states 
is discussed in Appendix 2. The last two columns 
give the multiplets of ~* and sp, respectively, on 
which the three electron states of sp* are built. 

Table VI gives the non-degenerate states of 
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Tasie IV. Configuration 


mM, mi Ms mi ‘ms mi Msg ML Multiplets Multiplets from pairs 
1/2 1 1/2 0 1/72 —1 3/2 0 4s 3°P 
1/2 1 1/2 0 —1/2 1 1/2 2 3/2 +3/2'D 
1/2 1 1/2 0 —1/2 0 1/2 1 1/2*D +1/2*P 3/2*P+7/6'D+1/3'S 
1/2 1 1/72 —1 —1/2 1 3/2 *P+7/6'D+1/3'S 
1/2 1 —1/2 1 —1/2 0 —1/2 2 2D 3/23P+3/2'D 
1 1 —1/2 1 —1/2 —1/2 1 1/2*D+1/2°*P 3/2 *P +7/6'D+1/3'S 
1 0 —1/2 1 —1/2 0 1/2*D +1/2*P 3/2§P +7/6'D+V¥3'S 
1/2 1 1/2 0 -1/2 1/2 0 1/3 4S +1/6°*D +1/2°*P 24P+2/3'D+1/3'S 
1/2 1 1/2 —1 —1/2 0 1/3 4S +2/3*D 2 
1/2 0 172 —1 —1/2 1 1/3 4S +1/6*D +1/2*P 2§P+2/3'D+1/3'S 
W(A) W(B) wc) Ms ML W(ABC) W(AB) +W(AC) +W(BC) 
TABLE V. Configuration sp*. 
s sp? Multiplets f i 
m, | Me m| m My Multiplets sp 
1/2 | 1/2 1 1/2 3/2 1 ‘P ap 2°P 
1/2 | 1/2 1 172 3/2 0 ‘P 1p 25P 
1/2 | 1/2 1 |-1/2 1/2 2 2D ip 3/24P+1/2'P 
—1/2 | 1/2 1 |—1/2 1} —1/2 2 Ip 3/2°P +12 'P 
1/2} 1 1 |—1/2 1/2 1 1/3 +1/2°D+1/6°P 1/2*P+1/2'D 3/2°P+W/2'P 
1/2 | 1/2 0 |-1/2 1 1/3 +1/22D +1/6°P 1/23P+1/2'D 
—1/2 | 1/2 1 1/2 1/3 +P +2/3°P P+ ip 
1/2 | 1/2 ij-1/72 1/2 1/3 +1/6 *D +1/6 *P +1/3 *S 1/2°P+1/6'D+1/3'S 3/24P+1/2'P 
1/2 | 1/2 0 2/3*D +1/3 2/3 'D+1/3'S 3/2°P+1/2'P 
172} 1/2 —1 |—1/2 1 1/3 +1/6 *D +1/6 *P +1/3 2S 3/24P+1/2'P 
—1/2 | 1/2 1 172 1/3 +P +2/3°P P+ ip 


Ms ML W(ABC) 


W(BC) W(AB) +W(AC) 


s*p and is very simple since the two s electrons 
form a closed shell. 

With the information available in Tables IV, 
V and VI it is now possible to write the energies 
of the three electron problem directly from Eq. 
(4). From the first state of Table IV, for example, 
we may write 

Wi 4S)=3W(p? *P)—3W(p?P). (25) 


It is important to note that these energies must 
be measured either with respect to the naked 
nucleus or to the core, if there is one inside the 
p electrons. From the second state of Table IV 
we have 

Tasie VI. Configuration s*p. 


s s sp Multiplets from 
pairs 
ms Ms m, Ms multiplets} s* sp 
172} —1/2 172 1 172 1 13/2§P+1/2'P 
172} —1/2 1/72 0 172 0 IS |3/2*P+1/2 
172 | —1/2 |-1/2 1/-1/72 1 |3/2*P+1/2 


W,(p* *D) = 3/2 W(p? *P) 

+3/2 W(p?'D)—3W(p?P). (26) 
W:(p*?P) can easily be found by elimination 
from the other states which are known in Table 
IV. The energies of the various multiplets of ’, 
sp’, *p and s*s’ found in this way are presented 
in Table VII. The W’s are omitted in Table VII 
for convenience. 


TABLE VII. Results for three electrons. 


sp 
—3°P 
2D= |3/24P+3/2'D 
*P+5/6'D+2/3'S —3°P 
= —2°P| —*S 
iP= 1/2*P+3/2'P| —2°*P | —*§ 
1p 3/24*P+1/2'P| | 
1S —2°P | — 
| ss’ sp sand s’ 
Ss 3/2 —238 
sts’ 2S = —*§ 


| 

| 
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TABLE VIII. Configuration p*. 


‘Mm, Mm, mi mim, mi MsMr Multiplets 
1/2 1 172 OF 1/2 —1}—1/2 1 1p 
1/2 1 1/2 1/2 1 
1/2 1 172 OF 1/2 —1)—-1/2 —1) 1 1p 

—1/2 Oj}—1/2 —1} 1/2 1 
1/2 1 12 1} —1/2 2 ip 
1/2 1 1/2 1)}—1/2 0 1} 1/24P+1/2'D 
1/72 1 172 O;—1/2 1;}-1/2 —1 1/2*P+1/2'D 
1/72 1 1/2 —1;/—1/2 1)—1/2 0} 2/3'D+1/3'S 
1/72 1 1/2 O}—1/2 -1/2 -1 1/24*P +1/6'D 

ol +1/3'S 

172 1/2 1)-1/2 same 
(A) (B) (C) (D) (ABCD) 


Four electrons 


The configuration p* gives rise to the multi- 
plets *P, '\D and 'S. Table VIII shows the 
various states of the individual electrons and 
the linear combination of the multiplets which 
represents jeach state with Ms+Mz,=0. Table 


-IX is a similar table for the configuration sp\ 


The energies of the various multiplets of p*, sp* 
and also s*p? in terms of the three, two and one 
electron energies are gathered in Table X. 
These relations include the single, pair and triple 
energies and should therefore be rather good 
approximations to the actual energies of the 
states. In the notation of the preceding section 


these relations are secdnd approximations and 
should be written, for example, W2(p* *P). 

The relations given in Table X are suitable 
for finding the energy of a state of either the - 
four, three, two or one electron system in terms 
of the others. It may happen that there is not 
sufficient information and it is desired to get a 
relation, say, between the four, three and two 
electron states. Relations which give the relative 
energies for this case can be obtained directly 
from Table X, by subtraction. If we wish’ the 
absolute energies we can get them in the fol- 
lowing way for example for p**P. From Table 
X we have . 


We(p* *P) = 4/3 (p* 4S) +5/3 (p* *D) 
+(p* *P) —4(p* *P) —5/3 (p* 1D) 
(p* +4(p*P). (27) 


The energies given here are all relative to the 
closed shell beneath. If we rewrite Eq. (27) 
giving the energies relative to the ion p*P we 
find that the terms (p*P) cancel with the 
exception of a sirtgle one. 
W2*(p* *P)=4/3 (p* 4S*) +5/3 (p* *D*) 
+ (p* *P*) —4(p* *P*) —5/3 (p? *D*) 
—d/3 (p? *S*)+(p*P). (28) 


TaBLe IX. Configuration sp’. 


? 

Ms mM, mi Ms mi Ms My sp® Multiplets 

1/2 v2 1/2 0 2 -1 | 2 6s 

1/2 1/2 1 1/2 0 —1/2 1 1 2 1D 
—1/2 1/2 1 1/2 0 —1/2 1 0 2 1/24D+1/2'D 

1/2 2 1 0 -1/2 1 . 1/24D +1/2'D 

1/2 2 1 1/2 0 -1/2 Of 1 1 1/24D +1/2°P 

1/2 2 1 1/2 #1 1/2°D +1/2°P 

1/2 oO | i 0 1/4*S +1/12 *S +1/6*D +1/2"P* 

1/2 12 0 12-1 

1/2 1/2 1 1/2 —1/2 0 1/4 +1/12 4S +2/3*D 
—1/2 2 1 1/2 1/2 1/4*S+3/44S 
—1/2 1/2 1 1/2 0 —1/2 0 0 1 1/4*D+1/4'D+1/4*P+1/4'P * 
—1/2 2 1 1/2 —1 -12 1 same 

1/2 1/2 1 —1/2 1 

1/2 1/2 0 1 -1/72 0 
—1/2 1/2 1 1/2 0 -1 0 0 1/6 *S + 1/6 *S +1/12 *D+1/12 +1/4*P +1/4'P 
—1/2 1/2 0 12-1 -12 1 same 

1/2 —1/2 1 —1/2 0 172 —1 

1/2 -172 0 —1/2 68 
—1/2 1/2 1 1/72 —1 —1/2 0 1/6 *S +1/6 4S +1/3 *D+1/3 

1/2 1 -1/2 -1 1/2 0 same 

1/2 1/2 0 2-1 -1 1/2"D +1/2"P 

1/2 —1/2 1/2"D +1/2°P 
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Taste X. Results for four electrons. 
sp? sD ? 5 
4/3 4*S+5/3*D+ —5/3 'D—1/3'S8 +4°P 
ID= 3°*D+ —8/3 'D—1/3'S +4°P 
iS = 4°P —5/3 'D—4/3'S +4°P 
4s +3 °P +728 
4s +1/3 +8/3 2'1P | +3°P +28 
2p +4/3 4P +1/6°P +3/2°D —3/2§P —3/2'D —5/2*P —1/2'P | +3°P +28 
ID= +3/2 °*P+3/2 2D —3/2*P —3/2'D —3/2§P —3/2'P| +32P| +28 
IiP= +4/3 +P + 1/6 +5/6*D +2/3 —3/2 *P —5/6'!D—2/3'S —5/2*P—1/2'P | +3°*P +28 
+3/2 *P +5/6*D +2/3*S —3/2 *P —5/6 —2/3'S —3/2§P —3/2'P| +3°P +38 
stp | sp? | sp | st | | s 
stp? = 2°P +4/3 +2/3°P +23P +228 
2°P + -1§ +2°P +238 
TaBLe XI. Configuration sp‘. 
Ms ms, mI Ms mi Ms mi Ms mi Ms ML Multiplets 
1/2 1/2 1 1/2 0 172 —1 —1/2 1 3/2 1 ‘P 
1/2 1/2 1 1/2 0 —1/2 1 -1/72 Oo 1/2 2 2D 
—1/2 1/2 1 1/2 0 1/2 —1/2 1 1/2 1 1/3 +2/3°P 
; 1/2 1/2 1 1/2 7) —1/2 1 -1/2 —1 1/3 +1/62P +1/22D 
1/2 1/2 1 1/2 -1 —1/2 1 same 
—1/2 1/2 1 1/2 72 -1 —1/2 1/2 0 1/3 +P +2/3 °P 
1/2 1/2 1 1/2 -1/2 O -1/2 -1 1/3 +P + 1/6 *P +1/6 2D +1/3 
; 1/2 1/2 0 1/2 -1 —1/2 1 —1/2 0 same 
1/2 1/2 1 1/2 -1 —1/2 1 -1/2 -1 2/3 +1/3 2S 
| (A) (B) (C) (D) (E) Ms ML (ABCDE) 


Now it is possible to eliminate this state by 
using the relations for p* from Table VII. This 
can be done in a variety of ways since p* gives 
three multiplets’ but if we use the relation for 
*S we find 


W:2*(p* *P)= 1/3 (p* *S*) +5/3 (p* *D*) 
+ (p? *P*) — (p? *P*) —5/3 (p* 
—1/3 (p? 1S*). (29) 
This expression cannot be expected to be as 


accurate as (27) since (27) is a second approxi- 
mation and we have obtained (29) by elimination 


with a first approximation. The accuracy of (29) 


will be about that of a first approximation 
relation such as those in Table VII. No table 
of the relations of this lower approximation is 


given, since they can be obtained directly from 
Tables VII and X by elimination. Furthermore, 
several relations can be obtained for each state 
and the most desirable one will depend upon 
the information available in the particular 
example at hand. 


Five, six and seven electrons 


In the set of five electron systems we consider 
the configurations sp*, p> and s*p*. Table XI 
gives the individual electron states for sp* and 
the linear combination of the multiplets ‘P, *P, 
2D and *S which characterizes each state. 
Several states which give no additional informa- 
tion about the energies are omitted. The energy 
relations for the various states of sp*, p> and 
s*p are given in the following equations. 


Ws(sp* #P) = + (sp®)[5/4 5S+5/3 8S] —(p*)[4/3 *S+5/3 *D+2P] 
— (sp®)[11/3 *P+1/3 *P+5/3 2D+1/3 2S] +(p2)[4 *P+5/3 ‘D+1/3 'S] 


* If we use any of the other relations the resulti ua- 


tions, though 


ifferent, are quite as good as (23). 


+(sp)(7/2 *P+1/2 1P]—4(p 2*P) —(s 


relative accuracy of these different expressions depends on 
the relative accuracy of the relations used from Table VII. 
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W;(sp* (pt *P) + (sp*)[4/3 *S+5/12 *D+5/4 1D+1/4 *P+3/4 
— (p*)[4/3 4S+5/3 *D+*P]—(sp*)[2/3 *P+-10/3 #P+5/3 *D+ 1/3 *S) 
+ (p*)[4 *P+5/3 'D+1/3 'S}+(sp)[2 *P+2 'P]—4(p *P) —(s *S), 
W;,(sp* *D) = (p* 'D) +(sp*)[9/4 *D+3/4 'D+3/4 *P+1/4 'P]—(p*)[3 *D+*P 
—(sp*)[2 #P+*P+8/3 *D+1/3 *S]+(p*)[3 *P+8/3 1D+1/3 
+(sp)[3 *P+'P]—4(p *P) —(s *S), 
W3(sp* *S) = (p* 4S) +(sp*)[3 *P+'P]—4(p* *P) — (sp*)[2 *P+2P+5/3 *D+4/3 2S] 
+ (p*)[3 *P+5/3 'D+4/3 'S}+(sp)[3 *P+'P]—4(p *P) —(s *S), 
p> W;(p' *P) = (p*)[3 *P+5/3 'D+1/3 'S]—(p*)[2 4S+5 *D+3 *P] 
+(p*)[6 *P+10/3 1D+2/3 'S]—5(p *P), 
W,(s*p* 4S) =3(s*p* *P) + (sp*)[5/4 °S+3/4 *S]—3(s*p *P) —(sp*)[4 *P+2 *P] 
— (p* *S) + (s? 4S) + (sp)[9/2 *P+3/2 'P]+3(p* *P) —2(s *S)—3(p *P), 
W;(s*p® *D) = (s*p*)[3/2 *P+3/2 1D ]+(sp*)[3/2 *D+1/2 ‘D]—3(s*p *P) 
—(sp*)[2 *P+3 *D+?P]—(p* *D) 1S) +(sp)[9/2 *P+3/2 'P] 
+(p2)[3/2 *P+3/2 ‘D]—2(s *S) —3(p *P), 
W;(s*p* *P) = (s*p?) [3/2 *P+5/6 'D+2/3 'S]+(sp*)[3/2 *P+1/2 'P]—3(s*p *P) 
—(sp*)[2 *P+2P+5/3 *D+4/3 *S]—(p +(s2 1S) 
+(sp)[9/2 *P+3/2 'P]+(p*)[3/2 *P+5/6 'D+2/3 'S]—2(s *S) —3(p *P). 
The lower approximations can be found also for these configurations but since they can be obtained 
in the straightforward way described above for p* they are not given here. 
Of the six electron configurations we consider sp*, p* and s*p*. The energies of the various multi- 


plets of these in terms of the energies of the higher ions, with all energies measured to the core 
beneath, are given in the following relations. 


°P) = (p**P) + [8/3 *P+1/3 2P+5/3 *D+1/3 2S] —(p)[3 *P+5/3 'D+1/3 1S] 

— (sp)[5/3 *S+1/3 *S+25/6 *D+5/6'D+5/2*P+1/2'P] 
+ (p*)£2 (sp*)(14/3 *P+4/3 *P+10/3 *D+2/3 *S] 
— (p*)[6 *P+10/3 '‘D+2/3 

W, (sp? 'P) = (p> *P) + (sp*)[3 *P+5/3 *D+1/3 *S]—(p)[3 *P+5/3 'D+1/3'S] 
—(sp*)[2 *S+5/2*D+5/2 'D+3/2*P+3/2'P]+(p)[2 
+ (sp*)[2 *P+4 2P+10/3 *D+2/3 *S]—(p*)[6 *P+10/3 ‘D+2/3 'S} , 

— (sp)[3*P+2 'P]+5(p*P)+(s*S), 


W.(p* 'S) = 6(p* *P) — (p")[9 *P+5 'D+'S]+(P*)[4 ‘S+10*D+6*P] 
— 


m . 
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Spt Wi(s*pt *P) = (s*p’)(4/3 *S+5/3 *D+*P]+ (sp*)[4/3 *P+2/3 *P] 
— (p* *P) —(s*p*) [4 *P+5/3 'D+1/3 1S] —(sp*)[5/3 §S+5S+5/2 *D+5/6'D 
+3/2*P+1/2 'P]+4(s*p *P) + (p*) [4/3 *S+5/3 ] 
+(sp*)[16/3 *P+8/3 *P+10/3 2D +2/3 *S]—(s* —(p*)[4*P+5/3 'D+1/3'S]} 
—(sp)[6*P+2'P]+2(s*S)+4(p*P), 


W.(s*p* 'D) = (s?p*) [3 *D+2P]+2(sp* 2D) — (p* — (s*p*)[3 *P+8/3 ‘D+1/3 'S] 
— (sp*)[9/2 *D+3/2'D+3/2*P+1/2 'P]+4(s*p*P) + (p*)[3 *D+*P] 
+ (sp*)[4 *P+22P+ 16/3 2D+2/3 2S] —(s? 1S) —(p*)[3 *P+8/3 'D+41/3 'S] 
— (sp)[6*P+2 'P]+2(s*S)+4(p?P), 
= 4(s*p? *P) +2(sp* 2S) — (p* — (stp) [3 *P+5/3 ‘D+4/3 'S]—(sp*)[6 *P+2'P] 
+4(s%p *P)+4(p* + (sp*)[4 *P+22P+10/3 2D+8/3 2S] 
—(s* 1S) — (p*)[3 *P+5/3 1D+4/3 'S]—(sp)[6 *P +2 'P]+2(s*S)+4(p *P). 


For the configuration (s*p*) we have 
s*p’ W,(s*p® *P)=(s*p*)(3 *P+5/3 'D+1/3 'S]+(sp*)[3/2 *P+1/2 1P]—(s*p*)[2 *D+3 ] 
— (sp*)[44P+2*P+ 10/3 *D+2/3*S]— (p> *P) + (s*p*) [6 *P+10/3 *D+ 2/3 *P] 
+ (sp*)[5/2 §S+3/2 *S+15/2 *D+5/2 1D+9/2 *P+3/2 'P] 
+(p*)[3 *P+5/3 1D+1/3 'S]—5(s*p *P) — (sp*) [8 *P+4 *P+ 20/3 2D+4/3 2S] 
—(p*)[2 *D+3 *P]+(s? 1S) +(sp)[15/2 *P+5/2 'P] 
+(p*)[6 *P+10/3 'D+2/3 'S]—2(s *S)—5(p *P). 


IIIb. INTER-MULTIPLET SEPARATIONS 


Relations have been developed above to give the absolute values of energy states in terms of those 
of higher ions. It is frequently the case among known energy states that the absolute values of 
several ions in succession are not known. Still it is desirable to know the relative separations in 
such cases. These can be obtained easily and directly for all the configurations which are studied 
above by subtraction of relations for states of the same configuration. Due to the equality of the 
sums of the coefficients for each configuration which occurs, these difference relations can always 
be written in terms of the energy differences between multiplets in the higher ions. 

It is also useful to have these difference relations in lower approximation. For example, the direct 
differences of the relations given above for sp* give the separation of two multiplets in terms of 
multiplet differences in sp’, p*, sp’, p°, sp and p*. A lower approximation can be found which gives 
the same separation in terms of the same differences in sp*, p*, sp? and p*. We find the coefficients 


as in section IIb. 


(ABCDE) = 2/32) W(aBy6) — 1/3) W(aBy) +1/3L W(a). (30) 


If we subtract two such expressions for different multiplets of the same configuration we find that 
the actual differences are the same as in the next higher approximation but all are multiplied by 
the coefficients in (30). The difference relations for the various configurations discussed above are 
collected below. Since the differences are all given from the lowest level of the configuration any 
other difference can be obtained by one subtraction. Some of the other differences give much simpler 


relations. 


> 
| 
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= (p*) [*S—*D ]= (6")3/2 
(p*)[*S—*P]= (p*)[S/6 (*P—'D)+2/3 (*P—'S)], 
sp? (sp*)[*P—*P]= (sp)3/2 [*P—*P], 
(sp*)[*P —*D]= (p*)[*P —'D]+(sp)1/2 [*P—'P], 
(sp*)[*P —*S]= (p*)[*P —'S]+(sp)1/2 (*P—'P] 
(p*)[*P —'S]= (p*)[4/3 (*S—*P)+5/3 
Next lower approximation : 
(p*)[*P— =1/2 (p)4/3 [*S—*D], 
—1S]=1/2 (p*)[4/3 (*S—*P)+5/3 (*D—*P)], 
(sp)[*S—4S]= (sp?)8/3 [4P—*P]— (sp)2[*P 
(sp*)[*S—*D ]= (p*)[*S—*D]+ (sp*)[1/6 (*P—*P)+3/2 (*P—*D)] 
— (p*)3/2 [*P—'P], 
(sp*)[*S—'D (p*)[*S—*D] + (sp*)[3/2 (*P—*P)+3/2 (*P—*D)] 
— (p*)3/2 [*P—'D]—(sp)3/2 
(sp*)[*S—*P]= (p*)[*S—*P]+(sp*)[1/6 (*#P —*P) +5/6 (*P—*D) +2/3 (*P—*S)] 
— (p*)[5/6 (*#P—'D)+2/3 [*P—'P], 
(sp*)[*S—'P]= (p*)[*S—*P]+(sp*)[3/2 (*P —*P)+5/6 (*P—*D) +2/3 (*P—*S)] 
— (p*) [5/6 (*P—'D) +2/3 (*P—'S)]—(sp)3/2 [(*P—*P]. 


The next lower approximation is obtained from the above by dropping those terms in two electrons 
(sp and p*) and multiplying the terms in three electrons by a factor 1/2. For example, for *S—*D 
(sp*)[*S—*D]= 1/2 (p*)[*S—*D]+1/2 (sp*)[1/6 (*P —*P) +3/2 (*P—*D)], 
s*p? (s*p*)[*P—'D (sp*)[4/3 ((P—*D)+2/3 (#P—*D) ]—(p*)[*P—'D], 
— 1S] = (sp*)[4/3 (*P—*S) +2/3 
For the next lower approximation drop the terms in ~* and multiply those in sp* by 1/2. 
sp* (sp*)[*P—*P]=(sp*)[S/4 (*S—*S)+5/4 (*D—'D)+3/4 (*P—'P)] 
— (sp*)3[‘P —*P}+(sp)3/2 (*P—'P], 
—*D]= (p*)[*P — 1D] + (sp*)[1/4 (*P—"P) +1/12 (*S—*D)+3/4 
+1/2 +2/3 (*P—*P)] 
+(p*)[*P—'D]+ (sp)1/2 [*P—'P)}, 
(sp) [*P—*S]= (p)[*P— 1S] + (sp*)[5/3 (*D—*P) +1/3 (*S—8P)+11/12 (6S—P) 


+1/12 (S—*P)+5/3 @D—*P)] 
(sp*)[2/3 (*P—*P) + (*P—*S) ]+ (p*)[*P — 1S] + (sp)1/2 [*P—'P}, 
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To obtain the next lower approximation drop the terms in sp and p* and multiply those in p* and 
sp* by 2/3 and those in p* and sp* by 1/3. For the lowest approximation drop the terms in p* and 


sp* as well and multiply those in p* and sp* by 1/3. 
s*p* —*D]= (s*p?)3/2 —'D] +(sp*)[5/4 (*S—*D) +1/4 (*S—*D) +1/2 (*S—'D)] 
— (sp*)[2(4P —*D) + —*D) ]—(p*)[*S—*D]+ (p*)3/2 'D), 
(s*p*)[*S—*P]= (s*p*)[5/6 (*P—*D) +2/3 (*P—'S) ]+(sp*) [5/4 (8S—*P)+1/4 — 
+1/2 (*S—'P)]—(sp*)[5/3 (“P—*D) +1/3 (*P—*S) + (?P—*S) ] 
— @P—*D)+2/3 (*P—'S)]. 


Follow the same procedure as for sp‘ to obtain the lower approximations. 
sp* (sp*)[*P—'P]=(sp*)8/3 [(*P (*S—*S)+5/3 

+ (sp*)8/3 [*P —*P]—(sp)[*P—'P]. 
To obtain the next lower approximation drop the two electron terms and multiply all five electron 
terms by 3/4, all four electron terms by 1/2 and all three electron terms by 1/4. For the next approxi- 
mation drop the three electron terms as well and multiply the five electron terms by 1/2 and the 


four electron terms by 1/6. For the lowest approximation drop the four electron terms also and 
multiply the five electron terms by 1/4. : 


s*p* [“S—*D] + (sp*)[4/3 (4P—*D) +2/3 
— (s*p*)[*P +1/3 @S—*D)+2/3 
+ (p*)4/3 [*S—*D]+ (sp*)[4/3 (*#P —*D) +2/3 
(s*p*)[*P —'S}= (s*p*) [4/3 (*S—*P)+5/3 (*D—*P) ]+(sp*)[4/3 (*P—*S)+2/3 ?#P—*S) ] 
— —'S]— (sp*)[5/3 @S—*P) + 
+3/2 @D—'P)+(*D—*P)+5/6 (‘D—*P) ]+(p*)[4/3 (*S—*P)+5/3 ?D—*P) ] 
+(sp*)[4/3 (*P—*S]+2/3 ?P—*S) ]—(p*)[*P—'S]. 
To obtain the lower approximations follow the same procedure as for sp*. 


IV. EXAMPLEs states. The examples considered here will be 
. ‘ . confined to the elements O, N and C largely 

doped because the beautiful work of Edlén® has made 
energy states of these elements so fully and 


of an atom or ion in terms of the energy states 
of the higher ions. In these relations we may use 
the values of the calculated energies‘ or we may 
use experimentally known states and from them 
determine other states. It is the latter application 
which is most useful since the energies have, for 
the greater part, been calculated only to the 
first order by the perturbation method and in 
general do not agree very well with the observed 


‘For example, if the Slater relations are known for two 
electrons they can be found for several electrons by using 
the relations which give the many electron energies in 
terms of those of two electrons. 


accurately known. Table XII gives the energies 
of various states of O, N and C relative to the 
1s? core, as taken from the work of Edlén. The 
spin-orbit separations have been reduced to the 


center of gravity. 
Three electrons 


Let us consider the state 2s2p??P. From 
Table VII we find the following relation 


*B. — Zeits. f. Physik 84, 746 (1933). This is a con- 


densed work. For complete information see Nova Acta 
wee Societatis Scientiarum Upsaliensis, Series IV, Vol. 9, 
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Taste XII. Energies to the 1s* shell for oxygen, nitrogen and carbon. 
Configuration Configuration 
and State Oo N c and State Oo N c 
1p 3,476,859 1p 2,801,355 
1s 3,458,934 is 2,756,342 
sp 3,369,401 stp *Pe 2,656,580 1,796,546 1,102,954 
sp? +P 2,585,463 1,739,779 1,059,945 
stp’ 4S° ; 3,382,889 2,152,548 i 2,476,194 1,650,703 992,344 
De 3,356,073 2,133,325 2p 2,529,895 1,695,636 1,028,064 
3,342,422 2,123,708 as 2,492,471 1,665,659 
3,262,959 2,064,415 4se 2,425,562 1,610,257 
‘ 3.170.240 2pe 2,401,669 1,593,583 952,531 
2D 3,216,898 ape 2,367,817 1,566,255 25 
2s 3,187,179 
stis 2,032,441 1,414,037 906,337 
ip 3, ‘079 374 2,020,019 1,183,619 sp *P* 1,950,515 1,347,114 853,969 
1s 3,056,461 2,002,648 1,172,164 ip 1,873, 1,283,342 803,825 
3,039,750 1,988,667 (1,159,500) 1,819,006 1,238, 768,912 
age 2,902,558 1,880,205 ly 013 ip 1,800,719 1,225,152 
ape 2,979,606 1,943,091 1,129,724 is 1,744,532 1,178,667 723,818 
ipe 2,912,595 1,891,146 
ape 2,957,258 1,926,115 1,118,557 s*S 1,114,009 789,538 520,178 
ipe 2.889, 186 1,868.569 
p*Pe 1,017,280 708,901 455,523 
Taste XIII. Three electron configurations. Absolute values to normal state of next ion. 
° 
Oo N c 
and State 
2s*2p °*P 617,295 624,139 6,844 374,365 382,509 Bide 187,067 196,617 9,550 
2s2p? +P 539,026 553,022 13,996 . $11,645 325,742 14,097 139,289 153,608 14,319 
2P 423,717 443,753 20,036 215,987 236,666 20,679 64,073 86,007 21,934 
2p 482,302 497,454 15,152 266,117 281,599 15,482 105,768 121,728 15,960 
2s 426,115 460,030 33,915 219,632 251,622 31,990 69,123 iT 31,042 
2p°4S 372,737 393,121 20,384 175,642 196,220 20,578 33 54,615 20,785 
2p 345,305 369,228 23,923 155,179 179,546 24,367 21,155 46,194 24,039 
2p 307,850 335,376 27,524 124,189 152,218 28,029 —3,272 27,915 31,187 
(2s2)2p? 4S 261,376 283,244 21,868 94,803 117,214 22,411 
231,271 256,428 25,157 71,961 97,991 6. 
2p 216,001 242,777 26,776 60,382 88,374 27,992 


(sp?) *P=(p*) *P+(sp)[1/2 *P+3/2 'P] 
—2(p) *P—(s) *S. 


-If we use this relation to determine the state 


(sp?) ?P, then we insert the correct energies on 
the right side from Table XII. For O IV we 
find (sp*) = 2,456,158 cm™ with respect to the 
1s? shell while the observed energy is 2,476,194 
cm; The discrepancy between the two is 20,036 


-cm™'. This appears to be a rather small error 


compared to the total energy to the shell beneath 
and it is even fairly small compared to the 
interaction energy of the three electrons, 672,375 
cm. Under the usual criteria for a perturbation 
problem this would be declared to be a good 
approximation. Unfortunately, however, we are 
seldom interested in the total energy but rather 


in the ‘‘absolute energy” to the normal state of 
the next ion. This absolute value calculated is 
423,717 and observed 443,753 cm. This error 
in the calculated value is due to the omission of 
the triple energy and we may therefore expect 
the discrepancy for the three electron cases to 
be greater than for more electrons where the 
triple energies are included. The results for 
configurations 2s*2p, 2s2p*, 2p* and (2s*)2p* (i.e., 
2p outside a closed 2s? shell), giving the calcu- 
lated and observed absolute values to the normal 
state of the next ion and their differences, are 
presented in Table XIII. The difference between 
the calculated and observed values (the triple 
energy) for a given state is nearly constant® with 


* The ities of the pair and triple energies as well 
as the total electrostatic interaction as a function of Z will 
be discussed in a later paper. 
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Tase XIV. The four electron configurations. Absolute values to normal state of next ion. 
N Cc 
State 
Calc. Obs. Obs.-Calc. Cale Obs. Obs.-Cale. Calc. Obs. Obs.-Calc. 
2s*2p? *P 440,345 442,602 2257 235,965 238,584 2619 88,155 90,786 2631 
420,342 422,536 2194 220,707 223,357 2650 77,910 80,622 2712 
1s 401,681 399,623 —20538 207,238 —1252 71,431 69,167 —2264 
2s2p°*S 382,399 382,912 513 191,449 192,005 556 55,894 
as 244,756 245,720 964 81,456 83,543 2087 —24,412 — 14,984 9428 
sp 322,816 322,768 48 146,062 146,429 367 26,053 26,727 674 
ip 253,993 255,757 1764 91,066 94,484 3418 — 14,027 a —_—-- 
301,469 300,420 — 1049 129,740 129,453 —287 17,827 15,560 —2267 
ip 232,648 232,348 —300 74,744 71,907 —2837 —22,326 — — 
158,228 158,943 715 18,191 
ip 144,657 144,517 —140 9, —_— 
1s 99,288 99,504 216 — 25,896 — 
108,516 109,761 1245 
1p 92,827 93,970 1143 
1s 74,787 76,045 1258 


Z. Thus if an adjacent element is known the 
actual correction (triple energy) can be extra- 
polated with considerable accuracy. 


Four electrons 

Since the triple energies are included in the 
relations for four electrons we expect them to 
agree better with the experimentally determined 
states and indeed this proves to be the case. 
Table XIV gives the absolute energies (to the 
normal state of the next ion) of the various 
multiplets of the configurations 2s*2p*, 252)’, 
2p* and (2s*)2p* as determined from the relations 
developed above and as observed, together with 
their differences. The discrepancy for s*p? in the 
absolute values is about 2500 cm~. The relative 
value *P—'D is, however, within 100 cm~ for 
all three elements. The '‘D—'S separation has a 
larger discrepancy than the absolute values. It 
is interesting to note how the differences (ob- 
served minus calculated) decrease with increasing 
Z, with the exception of 1S for N. The configura- 
tion sp* fits better than s*f* except for one or 
two multiplets in N and C which are not far 
from the ionization limit. In general we can 
expect these relations to hold best for the low 
terms of high multiplicity since the higher 
orders are not as important.’ In C, the sp**S 
state, which is important as the lowest tetra- 
valent state, has been predicted at 55,894 cm. 

7 It is possible to apply these relations to the deter- 
mination of the energy states of negative ions. The higher 
orders, however, may be important and the results there- 
fore may not be very accurate. For example, one finds that 
O should have an electron affinity in the 2p* state of about 


8000 cm or about 1 e. volt, whereas Lozier (Washington 
meeting 1934) gives 2.2+0.2 e. volts. 


This is 34,994 cm or 4.32 electron-volts above 
the normal state of C. The multiplets of p* 
(O III) agree somewhat better than either of 
the other configurations and for (2s?)2p* (normal 
state of O I), the absolute values are all about 
1200 cm too small, thus making the separations 
check with a maximum error of about 100 cm™. 
The error of only 15 cm in the *P—'S separa- 
tion is certainly fortuitous. 


Five and six electrons 

The ordinary absolute values with respect to 
the normal state of the next ion, of the multi- 
plets of configurations 2s?2p* of O and N and of 
2s2p* and 2s?2p* of O, determined by the 
methods above, the values obtained experi- 
mentally, and the differences between these are 
presented in Table XV. It has been pointed out 


Taste XV. Five and siz electron configurations. Absolute salues to normal state of 
next von. 


Obs.- 
Cale. | Obs. Cale. | Cale. | Obs. Cale. 


640 | 117,214 | —426 
2D 256,019 | 256,428 409 po} 97,991 | —592 


: 88.374| 1816 
«P 162,775 | 163,314 
68.964 | 70. 1631 
2D 115,860 | 117,253 | 1393 
28 85,804 | 87. 1730 
111,088 | 109,761 | —1277 
1D 95.472 | 93.970 | —1502 
1s 78,465 | 76,045 | —2420 


that the methods used here are more accurate 
the more electrons there are considered: thus 
p* is better than p*. This does not go on indefi- 
nitely, however, since as more electrons are 
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added the Oth approximation recedes furthér 
from the actual problem and higher orders 
become larger. The percentage error on the total 
energy measured to the closed shell beneath 
will continue to decrease but the actual error 
may not. There are also other perturbing effects 
due to spin-orbit interaction which become 
apparent in these higher approximations. The 
absolute energies calculated for 2s*2p* of O I are 
enough greater than the energies obtained when 
only the 2* group was considered that the 
differences between the observed and calculated 
values have negative signs, and indeed the 
deviations are not as small as they are for 2p*. 
This is a case in which the actual size of the first 
term not considered in the approximation is 
higher for s*p* than for p*, in spite of the fact 
that the former is a higher approximation. 


The atomic energy relations developed here have been 
used in the examples given above to calculate the absolute 
energies of atoms and ions in terms of the energies of 
higher ions. The spectroscopist is generally interested in 
the reverse process and these relations can be used as well 
for this purpose. Let us consider (2s*)2p* of O I and suppose 
that we wish to find the absolute value of 2s*2p*P in O IV. 
The energies of all the states which were used in the rela- 
tions for 2p* were with respect to the 2s* shell but if we do 
not know the O IV states we would know the states of the 
lower ions only with respect to 2s*2p *P,. Using energies to 
this limit (indicated by a prime) we can rewrite the relation 
for (2s*)2p**P in the following way: 
4(s*p)*P —3(s*p) *Py = (s*p) *P—24 *P = (s*p*) *P’ 

— (s*p*) [4/3 4S’+5/3 *D’+*P"] 

+(s*p*)[4 *P’+5/3 1D’+1/3 4S’], (31) 
where A*P is* the separation *Py—*Py, assuming the 
interval rule. If we neglect this correction we find 

(s*p) *P : 624,611 cm™, observed : 624,139 cm™. 

If we know the *P separation (as we do here) we would find 
a value (s*p) *P =625,384 cm™ which does not agree so 
well. The error in the corrected absolute value (1245) is 
about the same as the error for the normal state of O I, but 
the percentage error is very much smaller due to the large 
ionization potential. 

Let us consider the configuration s*p*. If we rewrite the 
relation for s*f**S so that the energies (indicated with a 
prime) are all measured to (s) *S; for the case in which 
absolute values in the highest ion are unknown, we find 
3(p) *P— 2(s) *Sy=(p) *P— +4 

== — (s*p*) 4S’ +3(s*p*) *P’ +(sp*)[5/4 §S’+3/4 

+(sp)[9/2 *P’+3/2 *P’]+3(p*) *P’, (32) 

*This complication is removed, of course, if A*P is 

known. 


where »(?S;—*P,4) is one of the resonance lines. Thus if 

we know the states in the lower ions and this one spectral 

line we can (neglecting 4 *P) immediately get the absolute 

values. For O we have 

(p) +4 *P = 823,693 

vobs. = 96,907.7 cm™, (p)*P+A*P = 1,017,508; 

observed : 1,017,280 cm™. 

Using the observed A*P we find 

*P =1,016,975, error : 305 

If we did not know this one line it would still be possible 

to eliminate either (s)*S or (p)*P between relations for 

two different configurations, e.g., s*p* and sp*, and solve 

finally for both. 

For the examples in O, N and C considered here, only 
those states have been considered which involve electrons 
with total quantum number n=2. For these elements the 
states which have an electron with higher value such as 
2s*2p*4s are not suited to this method since the electrostatic 
interaction of the inner electrons may easily be larger than 
the energy of the outer electron in the central field. In such 
a case the multiplets which arise from the same state in the 
ion lie very close together. The separation of the different 
groups built on different states in the ion can be found 
directly from the distance of these basis states in the ion. 


(b) In section III(b) it was shown that ex- 
pressions could be found for the differences of 
energy states within a configuration, in terms of 
the differences found in higher ions. To use these 
expressions absolute values are not necessary. 
The various multiplet separations of 2p?, (2s*)2p* 
and 2s2p* determined from these relations and 
observed for the elements O, N and C, are 
presented, together with their differences, in 
Table XVI. The differences between the ob- 
served and calculated values are roughly con- 
stant with Z, as they were for the absolute 
values. The errors for the three electron con- 
figurations are rather large but their regularity 
for different Z makes a much closer prediction 
possible if an adjacent element is known. 

The differences calculated and observed for 
the four electron configurations 2s*2p*, 252)’, 
2p* and (2s*)2p* are presented in Table XVII. 
For four electrons we have two approximations 
for the differences. The higher one (marked II) 
determines the separations in terms of differences 
in the three and two electron ions, and the lower 
one uses only the three electron ions. It is 
interesting to notice that certain separations 
[for example 2s*2p*(*P—'D)] are quite close in 
the lower approximation while other differences 
2s*2p?(‘\D—'S) are much more in error. The 


ve 

p* 

of 
ial 
ut 
ns 

to 

ti- 

of 

he 
ri- 
ire 
ut 
of 
2 
A 
te 
us 
fi- 
re 


964 R. F. BACHER AND S. GOUDSMIT 
TaBie XVI. Three electron configurations. Inter-multiplet separations. 
Oo N Cc 
Configuration 
and States 
: Cale Obs Obs.-Cale Cale. Obs Obs.-Cale Cale. Obs. Obs.-Cale. 
27,430 23,893 —3537 20,463 16,674 —3789 12,673 8,421 —4,252 
32p—tP 37,458 33,852 — 3606 30,997 27,328 — 3669 24,430 18,279 —6,151 
4S —1P 64,888 57,745 —7143 51,460 44,002 —7458 37,103 26,700 — 10,403 
(2s%) 4S 30,099 26,816 —3283 22,840 19,223 —3617 
15,275 13,651 —1624 11,581 9.617 — 1964 
q 4S —2P 45,374 40,467 —4907 34,421 28,840 —5581 
; 2s2p?+P —*P 115,308 109,269 —6039 95,658 89,076 —6582 75,216 67,601 —7,615 
J ‘P—2D 56,723 55,568 —1155 45,528 44,143 — 1385 33,521 31,881 — 1,640 
7 *?p—*P 58,585 53,701 —4884 50,130 44,933 —5197 41,695 35,720 —5,975 
4 1D —*S 56,187 37,424 — 18,763 46,485 29,977 — 16,508 36,645 21,562 — 15,083 
Tasie XVII. Four electron configurations. I nter-multiplet separations. 
oO N ¢ 
Configuration 
and States 
Calc Obs. Obs.-Calc. Calc. Obs. Obs.-Cale. Cale. Obs. Obs.-Cale. 
Il —'D 20,002 20,066 15,260 15,227 —33 10,245 10,164 —81 
Iip—'§ 18,661 22,913 4252 13,469 17,371 3902 6,479 11,455 4976 
a p—is 38,663 42,979 4316 28,729 32,598 3869 16,724 21,619 4895 
7 1 2s?2p24P —!D 19,145 20,066 921 14,452 15,227 775 9,347 10,164 817 
. Iip—is 37,424 22,913 —14,511 29,977 17,371 — 12,606 21,562 11,455 — 10,107 
i. sp—'s 56,569 42,979 — 13,590 44,429 32,598 —11,831 30,809 21,619 —9190 
4 Il 2s2p°4S—'S 137,640 137,192 —448 109,992 108,462 —1530 
59,590 60,144 554 45,388 45,576 188 
80,933 82,492 1559 61,708 62,552 844 
s3D—*P 21,343 22,348 1005 16,320 17,076 756 8,224 11,167 2943 
MH ID—'ID 68,820 67,011 — 1809 54,996 $1,945 —3051 
sp—ip 68,820 68,072 —748 54,996 57,546 2550 
I 145,692 137,192 —8500 118,768 108,462 — 10,306 
if ‘S—§D 62,728 60,144 —2584 48,869 45,576 —3293 
92,129 82,492 —9637 72,524 62,552 —9972 
29,401 22,348 —7053 23,655 16,976 —6679 
4 72,846 67,011 —5835 59,384 $1,945 —7439 
sp 72,846 68,072 —4774 59,384 57,546 — 1838 
13,570 14,426 856 
q 45,369 45,013 —356 
58,939 $9,439 500 
I 2p°*P—'D 15,928 14,426 — 1502 
ip—1s§ 50,778 45,013 —5765 
ap—is 66,706 59,039 —6767 
II (2s*)2p48P —!D 15,688 15,791 103 
Iip—1s 18,040 17,925 —115 
sP—1§ 33,728 33,716 —12 
I 17,877 15,791 
ip—1s 20,476 17,925 —4551 
38,353 33,176 —4637 


II denotes highest approximation; I next highest. 


The multiplet differences for the configuration, 


higher orders are apparently more important 


for the latter, or put another way, the ‘S is 
probably perturbed by adjacent levels. It is also 
interesting to notice that there are certain 
ratios of the separations holding for the lower 
approximation which still hold for the higher 
approximation. For example in 2s2p*, (§S—*S) 
is twice as large as either (*7D—'D) or (@P—'P) 
in the lower approximation.’ The same relation 


holds in the higher approximation. 


* This same ratio is obtained (See M. H. Johnson, Phys. 
Rev. 39, 209 (1932)) from the Slater-Condon relations. 


2s*2p', of O and N and for 2s2p* and 2s*2p* for 
O are presented in Table XVIII. For the five 
electron configurations there are three approxi- 
mations, III being the highest. Again (4S—*D) 
of 2s*2p° is good in the lowest approximation 
but the other separations of that configuration 
are not. Agreement closer than 200 or 300 cm™ 
must be regarded as fortuitous and it is doubtful 


One cannot conclude from the apparent validity of these 
ratios in this and similar configurations, that first order 
calculations should be expected to hold. 


| 
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TaBLe XVIII. Five and six electron configurations. Inier-multiplet separations. 


Element Configuration and States Obs. Cale. II | Cale. I | Cale. | Cale. II Cale. I Cale. I 
2s? 4S —2D 26,816 27,102 —286 26,891 —75 27,000 —184 
2D 13,651 14.209 —558 12,418 1233 20,167 16 
‘S —2P 40,467 41,311 —844 39,309 1158 47,167 —6700 
N 25293 4S —2D 19,223 19,054 169 19,403 — 180 19,540 —317 
9.617 12,026 —2409 9.790 -173 16,111 —0494 
‘S 28,840 31,080 —2240 29193 —353 35.651 —6811 
Oo 2s2p* +P —?P 92,719 93,808 — 1089 94,935 —2216 102,103 —9384 
46,061 46.913 —852 47.218 —11$7 —4093 
—t*S 29,719 30. —341 811 37,618 —7899 
PS 75,780 76,973 —1193 75,793 —18 87,772 | —11,992 
2s*2p43P —1D 15,791 15,907 —116 15,973 — 182 16,516 —725 
17.925 17,592 333 16.261 1664 25,008 —1173 
33,716 33,499 217 $2,254 1482 41.614 —7898 
Cale. IV @P 15,568; —1S) 17,002; @P —1S) 32,560 
Obs.-Cale. 223; 923; 1146 


For s*p? and sp‘, III denotes the highest approximation, II next highest, I lowest. 
For s*p* there are four approximations: IV highest, III next highest, etc. 


whether the greater error in the highest approxi- consider only the 2‘ outside the closed 2s 
mation for ‘S—*D of 2s?2p* has any significance. shell. There are four approximations available 
For 2s?2p* the greater error in the highest for 2s*26* and (*P—'D) is not bad even in the 
approximation indicates that it is better to lowest. 


APPENDIX 


The relations developed here will now be examined by using the methods of the perturbation theory. For three electrons 
with coordinates x,, x2: and x3, the Hamiltonian has the following form: 


H(x,X2%3) = Ho +0(x,%2%3) = Ho +-0( +0( x, (33) 


where Hp includes the central field due to the nucleus and v(x,x2x3) is the interaction of the three electrons. 

Let us consider the three electron state which is characterized by the quantum numbers A, B and C for the three 
electrons. Let ¥(ABC) represent the normalized wave function of this state including the first order terms. We can 
expand this ¥(A BC) in terms of the infinite set of one electron wave functions which are themselves properly normalized . 


B(x) C(x) a(x) (x) 
1 1 
W(ABC)=—|A(x2) B(x2) C(x2)| + «(abc)——|a(x2) b(x2) (34) 
v6 a,d.¢ v6 
A(xs) B(xs) C(xs) a(x3) 
sabe) VA BClabe) 


BC) — E,(abc) 


and a(x,) is the wave function of the state a of the one electron problem which has the same central field and the prime 
on the sum indicates that the term ABC is omitted. Eo(A BC) is the energy of the state A BC of the three electron problem 
if the interaction of the electrons is completely neglected. 


E(ABC) = E(A)+E(B)+E(C). 


The change in the energy of the state ABC due to the elect ostatic interaction, including terms of the first and second 
order, can be written 


E(ABC) = f ¥*°(A BC) dx, dx2dx3. (35) 
There will be no terms in E(A BC) in which the quantum numbers of the antisymmetric combination of one electron wave 


functions will differ for more than two electrons. On this basis the energy E(A BC) may be written in three parts, 
E(ABC) = (36) 


where E,(ABC), E:(ABC) and E,(ABC) contain only terms which differ in zero, one and two electrons, respectively. 


- — 
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For E,(ABC), E:(ABC) and E;(ABC) in terms of the expanded wave functions, we find 


E,(ABC) = V(ABC/ABC) (37) 
E,{ABC)= {x*(A Bn) V(ABC/ABn) +«(ABn) V(ABn/A BC) +x«*(AnC) V(ABC/AnC) +«(AnC) V(AnC/ABC) 
+«*(nBC) V(ABC/nBC) +«(nBC) V(nBC/ABC)} 
= {2«(A Bn) V(ABC/ABn) +2«(AnC) V(ABC/AnC) +2«(nBC) V(ABC/mBC)}, (38) 
n 
(39) 


E,(ABC) = {2«(Amn) V(ABC/Amn) +2«(mnC) V(ABC/mnC) +2«(mBn) V(ABC/mBn)}, 
m,n 


A(x;) m(x) n(x) 
A(x) B(x:) C(x2) A(x2) m(x2) (x2) 
A(x3) B(xs) C(xs) A(xs) m(xs) (xs) 


Since we wish to examine certain relations between the energies of the three electron problem and those of the two 
electron problem we must now study the latter. The central field will be taken as the same as that for the three electron 
problem so that if the interaction of the electrons is neglected, the energy states for a single electron are the same for 
both cases. The change in energy of the state AB due to the electrostatic interaction may again be written in three parts. 


E(AB) = E,(AB)+Ex(AB)+E,(AB), (40) 


where £,(AB), E;(AB) and E;(AB) contain only terms which differ in zero, one and two electrons, respectively. In terms 
of the expanded wave functions of the two electron problem (the \’s are the expansion coefficients), we find 


A(x) B(x:) C(x) 


v(x) 


where V(ABC/Amn) = 1/6 f 


E,(AB) = V(AB/AB), (41) 
E,(AB) == {2\AB/An) V(AB/An)+2\(AB/nB) V(AB/nB)}, (42) 
(43) 


E,(AB)= {2(AB/mn)V(AB/mn)}, 
m,n 
B(x) m(x,)  m(x,) 
Aas) | n(x) 
V(AB/mn) 
ond MAB/mn) B) — Eolmn) 


We find that the V's of the three electron problem can be written in terms of those of two electrons due largely to 


dx,dx- 


Eq. (44). 
= 0(x1x2) +0(x1x3) +0(x2x5), (44) 
V(ABC/ABC) = V(BC/BC)+ V(AC/AC)+ V(AB/AB), (45) 
V(ABC/ABn) = V(BC/Bn)+V(AC/An), (46) 
V(ABC/Amn) = V(BC/mn). (47) 

Using the above relations between the V's we can find the «’s in terms of the X's. 

x(Amn) = \(BC/mn), (48) 
«(A Bn) = \(BC/Bn)+(AC/An). (49) 


These relations enable us to write the electrostatic energies of the three electron problem in terms of those of the two 


electron problem. 
E\(ABC) = V(AB/AB) + V(AC/AC) + V(BC/BC) = E,\(AB)+E,(AC)+£,(BC). (50) 


This is a relation between the first order terms of the electrostatic interaction. If we include the energy in the central 


field we find 
W\(ABC) = W(AB)+W(AC)+ WBC) — W(A)— W(B)— W(C). (51) 


Since this is the same as Eq. (4), this shows that (4) is exact to the first order terms. For E,(A BC) and E;(ABC) we find 
E(ABC) = {2[\(BC/Bn) +(A C/An) V(BC/Bn) + V(AC/An)]+2[MA B/An) +X(BC/nC) 


X[V(AB/An) + V(BC/mC)}+2[MAB/mB) + (A C/nC) V(AB/nB) + V(AC/nC) J}, (52) 


E,(ABC) = E;(AB)+E£2(AC)+E:(BC) +cross terms in (52), (53) 
E,(ABC)= {2\(BC/mn) V(BC/mn) +2(A B/mn) V(AB/mn) +2d(AC/mn) V(AC/mn)}, 


(S4) 


E,(ABC) 


If it were not for the cross terms in (52) we would again get Eq. (4) when (53) and (54) are added to (51). Except for 
the cross terms in (52), Eq. (4) would hold up to the second order. We conclude that Eq. (4) which relates the energies 


| 
i 
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of the three, two and one electron problems, is exact up to the first order and holds also for parts of the second order. 
Let us now consider four electrons in the same central field. The Hamiltonian of the problem will have the form 


where the central field is included in H». The wave function (including first order terms) of the state which is characterized 
by the quantum numbers A, B, C and D for the four electrons may be written 


A(x;) B(x:) C(x) D(x) m(x:) o(x:) p(x) 


1 | A(x) Bs) Cle) Dies)! 1 | m(xs) (xs) plas) 
A(x,) Bix) C(x) D(xd) (xa) pled) 
V(ABCD/mnop) 


where yu(mnop) = =0 if mnop differs from ABCD in the quantum numbers of more than two 


BCD) — Eo(mnop) 

electrons. 

The change in the energy E(ABCD) due to the electrostatic interaction, including second order terms, may be split 
into three parts as before. 

E(ABCD) = E,(ABCD)+E(ABCD)+E,ABCD). (57) 

For £,(ABCD), E;(ABCD) and E;,(ABCD) we find 
E,(ABCD) = V(AB/AB)+ V(AC/AC)+ V(AD/AD) + V(BC/BC) + V(BD/BD) + V(CD/CD), 
E,\(ABCD) = E,(ABC)+E,(ABD) +£,(A CD) +E,(BCD) — E,(A B) — E,(AC) — E,(AD) — E,(BC) E,(BD)—E,(CD), (58) 
E,(ABCD) == {2¢(ABCn)[V(AD/An) + V(BD/Bn) + V(CD/Cn)] 


+24(ABnD)(V(AC/An) + V(BC/Bn) + V(CD/nD)}+2u(AnCD)[V(AB/An) + V(BC/nC) + V(BD/nD) 
+2u(nBCD)[V(AB/nB) + V(AC/nC) + V(AD/nD))j, 
E,(ABCD) = E,(ABC) + E,(A BD) + E,(A CD) + E,(BCD) — B) — E,(AC) — E,(AD) — E,(BC)— E,(BD)— ECD), (59) 
E,(ABCD)= {2u(ABmn) V(CD/mn) +2u(AmCn) V(BD/mn) +2u(mBCn) V(AD/mn) 
+2u(mBnD) V(A C/mn) +2y(AmnD) V(BC/mn) +2u(mnCD) V(AB/mn)}, 
E,x(ABCD) = E,(ABC) +E,(A BD) +E,(A CD) E,(BCD) — (AB) — E,(AC) — E,(AD) — E,(BC)—E,(BD)—E,(CD). (60) 
If we include the energy of the unperturbed problem we find for the total energy in second order. 


W*(ABCD) = W*(ABC) + W*(ABD) + CD) +W*(BCD) — WAB)— WAC) — WAD) — 
— W*(BD)— W*(CD)+W(A)+W(B)+W(C)+W(D). (61) 


This relation is the same as Eq. (7) so we conclude that (7) is exact up to and including the second order. It is very 
probable that (7) holds for parts of the third and higher orders. Up to the second order one can verify the relations for 
more than four electrons in the same way, although the equations become somewhat lengthy. Inclusion of the third and 
higher orders, however, greatly increases the complication. 


APPENDIX II 


In the configuration p* consider the two states with Ms=0, M,=1. We know that these two states together give *P 
and 'D, each being some linear combination of the two. In order to find this linear combination, let us consider in Table 
IV for p* the two states which are *D alone. These two states must have the same pair energy. If we write the pair energy 
of two electrons as w(m, m; m,' m,’) and the total pair energy for *D as w(?D) we have the following relations 


w(?D) =w(1/2 1, 1/20)+w(1/2 1, —1/2 1)+w(1/20, —1/2 1), (62) 
1D 
w(?D) =w(1/2 1, —1/2 1) +w(1/2 1, — 1/2 0)+w(—1/2 1, —1/2 0). (63) 
ip 


In (62), the first pair energy is that of *P of ~* (compare Table II) and the second that of 'D. In (63), the first is ‘D and 
the last *P. Since the pair energies must be equal in both cases we have (64) 
w(1/2 0, —1/2 1) =w(1/2 1, —1/2 0). (64) 


If we go back now to the two states of p? having Ms=0 and M,=1, we see that (64) shows they have the same pair 
energy and hence the two states must have the same linear combination of *P and 'D. The two states together give *P +'D 
so each state is 1/2*P+1/2 D. 
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Let us consider the three states with Ms=0 and M,=0 from the configuration p*. We know that these three states 
together give *P, 'D and 'S. If we go to Table VIII for configuration p* there are four states which are uniquely p* *P. 


Since these must have the same pair energy Eqs. (65), (66) and (67) result. 


*P) =w(1/2 1, 1/2 0) +w(1/2 1, 1/2—1) +w(1/2 1, —1/2 1) +w(1/2 0, 1/2—1) +w(1/2 0, —1/2 1) 
+w(1/2—1, —1/21), (65) 


w(p* *P) = w(1/2 1, 1/2 0)+w(1/2 1, 1/2—1) +w(1/2 1, — 1/2 0)+w(1/2 0, 1/2—1)+w(1/2 0, —1/2 0) 
+w(1/2—1, —1/20), (66) 


w(p* *P) =w(1/2 1, 1/2 0) +w(1/2 1, 1/2—1) +wl1/2 1, —1/2—1)+2(1/2 0 1/2—1)+w(1/2 0, —1/2—1) 
+w(1/2—1, —1/2—1). (67) 


Utilizing knowledge of the pair energies of other states of p*, and noting that if the signs of all the quantum numbers in 
a pair are changed the pair energy is unchanged, we get the following relations: 


w(1/2 0, —1/2 0) =1/2 w(*D)—1/2 w@P)+w(1/2 1, —1/2—1), (68) 
w(1/2—1, —1/2 1) =w(1/2 1, —1/2—1). (69) 

In addition the three states together give *P, 'D and 'S. 
w(1/2 0, —1/2 0)+2w(1/2 1, —1/2—1) =w(*P)+w('D)+w('S). (70) 

Hence we have (71) and (72) 

w(1/2 1, —1/2—1) =1/2 w(@P)+1/6 +1/3 w('S), (71) 
w(1/2 0, —1/2 0) = 2/3 wD) +1/3 w('S). (72) 
This means that in the set Ms=0, M,=0 the first and third states are 1/2*P+1/6'D+1/3'S and the second is 


2/3 'D+1/3 3S. 

Eqs. (64) and (69) indicate that the division into multiplets is the same for two states which differ either in the signs 
of all the m, or all the m;. The state which differs from another state in the signs of all the m, and all the m; is simply 
the corresponding state from the half of the Pauli table with Ms+M ,<0 and also has the same division into multiplets. 


From the configuration #* there are multiplets 4S, *D and *P. In the set with Ms=1/2 and M,=1 (see Table IV), 
we find that each of the states has the same pair energy and hence each state is given by 1/2 *D+1/2*P since the two 


together give *P+*D. The same is true for the set with Ms=—1/2, and M,=1. 
In the set with Mgs=1/2 and M,=0 there are three states which together give ‘S, *D and *P. Since we know that the 


pair energies of a state (ABC) must be the same as for (AB)+(AC)+(BC), we have 
w(A BC) =w(AB)+w(AC)+w(BO), (73) 


where w(ABC) is used here to represent the pair energies of (ABC). Utilizing (73), we find from the states which are 
already known in Table IV the following relations 


w(p* 4S) = 3w(p* *P), (74) 
w(p? *D) = 3/2 w(p? *P)+3/2 w(p* (75) 
*P) =3/2 w(p* *P)+5/6 +2/3 w(p? (76) 
From the three states with Ms=1/2 and M,=0 we have 
w(1/2 1, 1/20, —1/2—1) =2w(p* *P) +2/3 w(p? !D)+1/3 w(p* 4S), (77) 
w(1/2 1, 1/2—1, —1/2 0) =2w(p? *P)+w(p? (78) 
w(1/2 0, 1/2—1, —1/2 1) = 2w(p? *P)+2/3 w(p* 'D) +1/3 w(p? (79) 


We wish to find w(1/2 1, 1/20, —1/2—1), w(1/2 1, 1/2—1, —1/2 0) and w(1/20, 1/2—1, —1/2 1) in terms of w(p* 4S), 
w(p?*D) and w(p*?P) so that all the Eqs. (74) to (79) are satisfied simultaneously. If we solve (74), (75) and (76) for 
w(p* *P), w(p? DD) and w(p* 4S) and substitute these values in (77), (78) and (79), we get the following 


w(1/2 1, 1/2 0, —1/2—1) =1/3 w(p* 4S) +1/6 w(p**D) +1/2 w(p* *P), (80) 
w(1/2 1, 1/2—1, —1/2 0) = 1/3 *S) +2/3 w(p* *D), (81) 
w(1/2 0, 1/2—1, —1/2 1) =1/3 w(p* *S) +1/6 w(p* 2D) +1/2 w(p?*P). (82) 


The relations (80), (81) and (82) for the pair energies give the proper linear combinations found in Table IV for the 


three states of the set Ms=1/2, M,=0. 

It may seem that this is somewhat involved, but it is nearly always possible, with a little practice, to determine the 
coefficients immediately from inspection. For example, in the last case it is apparent since ~* *P occurs (see Table IV, last 
column) only with ~'S, that it will occur in the first and third states of the set with coefficient 1/2. The second state 
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has a 'D in the multiplets from pairs but no *P in the »* multiplets, so *D must occur with the coefficient 2/3. The first 
and third states must have the same linear combination and the sum of the *D coefficients is unity so *D occurs in the 
first and third states with coefficient 1/6. Since for each state the sum of the coefficients must be unity, the coefficient of 
*S is 1/3 for each state. 

It is not necessary to know the exact linear combination of all the states of Table IV in order to find the energies of 
the various multiplets but the information is useful in finding the linear combinations in more complicated configurations. 


From the configuration sp* there are multiplets ‘P, *P, *D and *S. Each of the states in Table V with Ms =3/2 is 
assigned to ‘P alone and those with M,=2 belong to *D. In the set with Ms=1/2 and M;,=1 three of the possible 
four states will be represented, ‘P, ?P and *D. Using (73) and the information available about ‘P and *D we find 


w(sp? =w(p* *P)+2a(sp *P), (83) 
w(sp* *D) =w(p*? 'D)+3/2 w(sp *P)+1/2 w(sp'P), (84) 
Considering the whole set Ms=1/2, My, =1 we have 
w(sp* *P) +w(sp* *D) +w(sp* *P) = 2w(p* *P) +w(p* +4w(sp *P) +2w(sp *P). (85) 
With (83) and (84) this gives : 
w(sp**P) =w(p* *P)+1/2 w(sp *P) +3/2 w(sp *P). (86) 


From the individual states of the Mg=1/2, My, =1 set, we have 
w(1/2, 1/2 1, —1/2 0) =w(1/2, 1/20, —1/2 1) =1/2 w(p* *P) +1/2 w(sp *P)+1/2 w(sp'P), (87) 
w(1/2, 1/2 1, 1/2 0) =w(p* *P) +w(sp *P) +w(sp *P). (88) 
If now we solve (83), (84) and (86) for w(p?*P), w(sp*P) and w(sp'P) in terms of w(sp* *P), w(sp**D), w(sp**P) and 
w(p?'D) and put them into (87) and (88) we find 
w(1/2, 1/2 1, —1/2 0) =w(1/2, 1/20, —1/2 1) = 1/3 w(sp* *P) +1/2 w(sp* *D) +1/6 w(sp* *P), (89) 
w(1/2, 1/2 1, 1/2 0) =1/3 w(sp*? *P) +2/3 w(sp* *P). (90) 
These relations for the pair energies give the linear combinations for the three states with Ms=1/2 and M;=1 found 
> pr vd of states with Mg=1/2 and M,=0 all four possible multiplets, ‘P, *P, *D, and *S, are represented. Con- 
sidering the whole set together we have 
w(sp? +w(sp* *P) +w(sp* *D) +w(sp? *S) = 2w(p* *P) +w(p? +w(p? 4S) +11/2 w(sp*P)+5/2 w(sp'P). (91) 
With the help of (83), (84) and (86) we find from (91) 
w(sp* *S) =w(p? *S)+3/2 w(sp *P)+1/2 w(sp (92) 
For the four states we get the following relations for the pair energies 
w(1/2, 1/2 1, —1/2—1)=w(1/2, 1/2—1, —1/2 1) =1/2 w(p* *P)+1/6 *D) 
+1/3 w(p? *S)+3/2 w(sp*P)+1/2 w(sp'P), (93) 
w(1/2, 1/20, —1/2 0) =2/3 w(p* +1/3 w(p* 1S) +3/2 w(sp *P)+1/2 w(sp *P), (94) 
w(—1/2, 1/2 1, 1/2—1) = *P)+w(sp *P) +w(sp 'P). (95) 
Solving (83), (84), (86) and (92) for w(p**P), w(p 'S), w(sp*P) and w(sp'P) in terms of w(sp**P), w(sp**P), 
w(sp?*D), w(sp?*S) and w(p*?D) as above and putting them in (93), (94) and (95) we find 


w(1/2, 1/2 1, —1/2—1) =w(1/2, 1/2—1, —1/2 1) =1/3 w(sp* *P) +1/6 w(sp? *P)+1/6 w(sp* +1/3 w(sp* *S), (96) 
w(1/2, 1/2 0, —1/2 0) =1/3 w(sp**S)+2/3 w(sp? 2D), (97) 
w(—1/2, 1/2 1, 1/2—1) =1/3 w(sp* *P)+2/3 w(sp* 2P). (98) 


The relations (96), (97) and (98) give the proper linear combinations found in Table V. For this configuration the in- 
spection method gives the coefficients almost immediately. 

No further particular examples are worked out here in detail as we believe the method is amply illustrated above. The 
coefficients can nearly always be found most easily by the inspection method for which a certain amount of familiarity 
with the Pauli tables is indeed helpful. There are numerous regularities to guide one which unfortunately cannot be 
pointed out here in detail. 
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Certain new infrared absorption bands of CH;Cl were 
mapped in the region from 0.74 to 7.0u. In each case the 
bands show the structure predicted by the theory of the 
symmetry properties of molecules having three equal 
atoms. It is also predicted that every third line in the 1 
bands shall be enhanced. This was observed in the present 
investigation as well as previously by Bennett and Meyer. 
Several new intervals between the lines in the | bands 
were observed which may be of value with regard to a 


theoretical interpretation of this anomalous spacing. Two 
absorption regions, observed here, appear to give evidence 
of an interaction between the levels »; and 2» thought to 
exist in some of the methyl compounds. These bands are 
v2 2x4) at 5900 and »3+(», 2») at 4200 cm. A 
few lines of »; in the P branch near the center were resolved. 
The Q branch proved to be much more pronounced than 
Bennett and Meyer's work shows. 


I, INTRODUCTION 


HE infrared vibration-rotation spectrum of 

the methyl! halides was measured by Ben- 

nett and Meyer' who found that the | bands 
showed an anomalous spacing of the lines. This 
effect, they suggested was due to an interaction 
between the vibration and the rotation of the 
molecule. A pair of || bands, separated by 
85 cm™, was also found at 3.24 which have 
recently been interpreted as (», 2»,).2 Normally 
the levels y; and 2», would lie very close to- 
gether. It is thought that in a number of methyl 
compounds they interact as in the case of the 
CO; * molecule and that the observed separation 
is a result of this interaction. This resonance also 
serves to increase the intensity of the || com- 
ponent of 2», so that it is comparable with ». 
Only in the case of »; for CHsF was resolution 
obtained by Bennett and Meyer so that the 
moments of inertia of the remaining molecules 
are unknown. The line of further investigation 
pointed out by the previous observations is as 
follows: (1) to acquire more experimental data on 
the anomalous spacing of the lines in the 1 
bands; (2) to obtain more experimental evidence 
of the supposed interaction between »; and 2»,; 
(3) to determine, if possible, the moments of 
inertia A =B from the resolution of one of the | 
bands; (4) to obtain experimental data on the 
complex structure of the harmonics and com- 
bination bands of | vibrations predicted by the 

1 Bennett and Meyer, Phys. Rev. 32, 888 (1928). 

2 The bracket (», 2»,) indicates that the wave functions 
representing the two perturbed levels are linear combina- 


tions of the wave functions of the levels » and 2»,. 
*D. M. Dennison, Phys. Rev. 41, 304 (1932). 


theory of molecules having three equivalent 
particles.* 

Two infrared spectrometers, one of the Pfund® 
type designed by J. D. Hardy, and another 
designed by C. F. Meyer® were available for 
this investigation. Three gratings ruled with 
7200 lines per inch, 4800 lines per inch and 2400 
lines per inch were used for the 1.5u—3.0y, 4.0u, 
and the 7.0u regions, respectively. In addition, 
a concave grating of 15,000 lines per inch in a 
Paschen mounting was used in the photographic 
region. For the bolometric measurements a 25.4 
cm cell closed with rocksalt windows, and a 1.5 
meter cell having glass windows were used. In 
the photographic region the equivalent length of 
the absorption cell was 100 meters at atmos- 


pheric pressure. 


II. EXPERIMENTAL RESULTS 


The 1.6y region 

Two | bands and one | band were found as 
shown in Fig. 1. The lowest frequency | band 
occurs alone while the other | one overlaps the 
|| band. The predicted structure of the second 
harmonic of a | vibration is a pair of close 
bands of | and | character. This is just what is 
observed in curves 5, c, Fig. 1, hence this pair 
has been designated as 2¥2. The spacing of the 
lines in the | component is 12.0 cm™. The zero 
branch of the || component is quite sharp and 


*D. M. Dennison, Rev. Mod. Phys. 3, 338 (1931). 
‘A. H. Pfund, J. Opt. Soc. Am. and Rev. Sci. Inst. 15, 


69 (1927). 
*E. F. Barker and C. F. Meyer, Trans. Farad. Soc. 25, 


913 (1929). 
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5800 cm 5900 cm* 6000. 6100 6200 


Fic. 1. Curve (a) shows the perpendicular structure of 
the band »,+(, 2»,). Curves (b) and (c) show the splitting 
up of the parallel, and perpendicular components of 2». 


more intense than any found by Bennett and 
Meyer (with the exception of »;). The other L 
band is thought to be ve+(v:, 2») involving 
the two resonance levels. This identification 
would indicate that the two || components of 


Taste I. (Data on the || bands.) 


P branch Q branch R branch »P—¥R 
4025.7 cm™ 4047.1 cm™ 4056.3 30.6 
4078.0 4089.5 4106.3 24.3 
4213.8 4230.6 4245.9 32.1 
6003.4 6015.2 6034.0 30.6 


TABLE II. (Data on curves a and c, Fig. 1.) 


Observed » Av Observed » Avy 
5855.1 cm™ 5920.4 cm™ 
6.7 9.1 
5861.8 $929.5 
4.5 9.0 
5866.3 5938.5 
6.7 7.0 
5873.0 5945.5 
8.9 9.2 
5881.9 5954.7 
9.0 6.9 
5890.9 5961.6 
68 6.9 
5897.7 5968.5 
9.0 9.3 
5906.7 5977.8 
6.0 
5912.7 
7.7 
6048.3 6110.7 
11.9 12.2 
6060.2 6122.9 
12.0 11.0 
6072.2 6133.9 
12.0 8.6 
6084.2 6142.5 
14.4 
6098.6 
12.1 
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1967, 185% 
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Fic. 2. Perpendicular structure to which no interpretation 
can at present be given. 


vo+2v, are much weaker than the two | com- 
ponents since the region shows only | structure. 
The data on this region are given in Tables I 
and II. For this region the 25.4 cm cell was 
used, and, since the background was free from 
atmospheric absorption, galvanometer deflec- 
tions were plotted against cm. 


The 1.8y region 

The absorption in this region is shown in 
Fig. 2. The band is essentially of the 1 variety, 
but it is so complex that no interpretation can at 
present be given to it. The absorption is con- 
siderably weaker here than at 1.6u necessitating 
the use of the 1.5 meter cell. The curve is plotted 
in percent absorption because the background is 
complicated by the strong water vapor absorp- 
tion. The complex appearance of this region is 
no doubt partly due to distortion by the water 
vapor lines. 


The 2.4y region 

A. Fig. 3 shows the highest frequency absorp- 
tion of the group of bands in the region of 2.4y. 
Its structure is | throughout. This absorption is 
believed to be composed of two overlapping L 
bands, namely, »:+¥; having its center at 4383.8 
and with its center at 4478.5 cm™. 
Theoretically both and would be 
single | bands, and this seems to agree with 
observation in that the edges of the region have 
sharply defined lines but the middle is distorted 
from overlapping. The average spacing of the 
lines in Table III is seen to be 8.5 cm™. The 
curves show clearly that every third line is 
enhanced as predicted by the theory of systems 
having three equivalent nuclei each with a spin 
angular momentum of }(4/27). As in the 
region the 25.4 cm cell was used and galva- 
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TABLE III. (Data on the 1 bands of Figs. 3, 4, 6.) ee 
Observed » A» Observed » AN 
4323.9 4441.6 cm™ V 
9.5 85g ‘ 
4333.4 4449.1 ly J 
4341.7 4457.9 3 /¥\ 
8.4 8.8 Center Center 
4350.1 4406.7 88 
6.0 7.7 Fic. 3. Two overlapping perpendicular bands, »:+¥; and 
4364.4 4483.2 va 
7.2 9.0 
4371.6 4492.2 
. . 2439, 2409 2.38) 2352 
7.3 9.0 
4387.4 4510.1 
9.8 10.4 Se 
4397.2 4520.5 
7.3 10.3 
4404.5 4530.8 de a 
8.6 7.9 
4413.1 4538.7 3 
10.5 6.5 
4423.6 4545.2 & 
10.5 3 
4434.1 Rad 
7.5 "20 
Mean spacing of lines in low frequency side of band 
=8.4 cm", ° 
Mean spacing of lines in high frequency side of band 
=8.6 
4123.9 cm™ 4182.3 cm™ cm" 42 00cm" 4250 
12.8 4.9 
4136.7 4187.2 Fic. 4. Perpendicular and parallel structure showing the 
11.3 5.0 separation of the bands »;+(m, 2»). 
4148.0 4192.2 
9.7 6.6 
4157.7 4198.8 nometer deflections were plotted against cm~. 
3.3 6.7 The slit widths included 1.6 cm~ and readings 
4161.0 4205.5 
43 g.3 were taken at intervals of 1.2 cm™. 
4165.3 4213.8 B. At somewhat longer wave-lengths a pair 
4169.1 3.8 of bands were observed of which one is || and 
13.2 the other 1. These bands, shown in Fig. 4, 
33807 cmt have been designated as 274) which is 
81 9.5 the second observed group thought to contain the | 
2397.8 is 2466.4 i resonance pair of levels. This identification is 
2407.6 ; 2476.2 “made because the separation of these bands is | 
10.6 10.2 just the separation expected for a combination of 
2418.2 2486.4 97 some other level with the pair 2»,). Theo- 
2427.9 2496.1 “ retically v3+» is a || band, and v3+2», is com- 
9.5 10.1 posed of a || and a | component lying close 
2437.4 96 2506.2 together. This is in agreement with observation 
2447.0 ' since one || band is found, and the | structure | 
9.9 shows such distortion as might be expected from ‘ 


overlapping. Again in this case the zero branch 
of the parallel band is very sharp and strong, in 


Mean spacing of lines=9.7 cm™. 
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Fic. 5. One perpendicular and two parallel bands 
comprising »%. 


contrast to Bennett and Meyer's observations. 
The absorption was observed with the aid of the 
1.5 meter cell and the curve was plotted in 
percentage absorption because of some strong 
atmospheric lines. The slits included 1.5 cm~ 
and readings were taken at intervals of 1.6 cm™. 
Table III gives the frequencies of the lines in the 
1 band, while Table I contains the data on the 
|| one. 

C. The lowest frequency portion of the group 
at 2.4u is shown in Fig. 5 and is composed of a 
1 and two | bands all overlapping to some 
extent. This group has been identified as v2+ 
since the character of the observed band agrees 
with that predicted for such a combination. 
The three bands are all of about the same in- 
tensity, but since the perpendicular part was 
plotted in percent absorption the ordinates are 
not directly comparable. Due to intense water 
vapor lines in the region, the authenticity of 
this band was for a time doubted. However, the 
gas was dried so carefully and such precautions 
were taken as to remove this doubt. Table IV 


TABLE IV. (Data on curve b of Fig. 5.) 


H,O H,O 
Observed lines Observed lines 
(em~) Av Av 


3963.0 3962.8 —0.2 3994.1 3995.7 1.6 
3968.9 3970.3 0.14 4003.1 


3976.3 3976.6 0.3 4010.6 40093 —1.3 
4020.7 4.0 
3983.7 39836 —0.1 4016.7 
4012.2. —4.5 
3991.1 3991.2 0.1 4024.2 4026.7 2.5 


408, 4000. 

1 l 1 i 

245 


Fic. 6. The perpendicular structure of »+2»;. 


shows the close agreement of some of the 
CH;,CI lines with water vapor lines.’ Another 
evidence indicating that the lines are due to 
CH, Cl is the fact that every third line shows the 
characteristic enhancement. The zero branches 
of the parallel bands chosen to lie at 4047.1 cm= 
and at 4089.5 cm are not very well defined. 
The region was measured with the 25.4 cm cell 
and the slit widths included 1.4 cm™, readings 
being taken at intervals of 1 cm™. 


The 4.0u region 

Only the one | band shown in Fig. 6, was 
found in this region. It is very clean throughout 
giving clearly an average spacing between lines 
of 9.7 cm~'. The 1-1-2 intensity characteristic 
is shown more clearly in this band than in any 
other found. This band is believed to be vg+ 2,5 
which according to theory must be a single 1 
band and should lie more free from overlapping 
than any other band in this region. The fre- 
quencies of the lines are given in Table III. 
Mapping of this band was accomplished with 
the 25.4 cm cell and the slits were set to include 
2.5 cm~. Galvanometer readings taken at in- 
tervals of 1 cm were plotted against cm™. 


The 7.0u region 


Resolution of vs; would give the moments of 
inertia A=B of the molecule. As Fig. 7 shows, 
this was not wholly successful despite the refine- 
ments in modern apparatus. Only a few lines in 
the negative branch near the center of the band 
were resolved giving a spacing of 1.3 cm™. 
These provide insufficient information for a 
reliable determination of the moments of inertia. 


7 Plyler and Sleator, Phys. Rev. 37, 1493 ¢1931). 
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1310 cme’ 1330 cm’ 1350 cm 1370 co’ 
Fic. 7. This curve shows the result of the new attempt to 
resolve v5. 


It is believed that the failure is due partly to 
the isotope effect in this vibration, partly to 
strong water vapor absorption and partly to 


TABLE V. (Data on photographic bands.) 


Observed Ap Observed » Ap 
9992.6 cm™ 10061.4 
5.8 16.8 
9998.4 10078.2* 
3.3 6.9 
10001.7 10085.1* 
7.4 6.8 
10009.1 10091.9* 
9.2 8.7 
10018.3 10100.6* 
7.0 6.3 
10025.3 10106.9* 
7.8 4.0 
10032.1 10110.9 
7.7 2.4 
10039.8 10113.3 
4.0 7.4 
10043.8 10120.7 
8.5 3.9 
10052.3 10124.6 
10057.4 
4.0 


Mean spacing of lines = 6.7 cm™. 
* Lines most accurately measured. 


13778.6 13827.6 

28.0 9.7 
13806.6 13837.3 
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inadequate resolving power. The zero branch, 
however, has a very different appearance from 
that found by Bennett and Meyer, being very 
sharp and as high as the positive and negative 
branches of the band. Its position agrees very 
exactly with that determined in the earlier work. 
The slits include 0.67 cm™ and galvanometer 
readings were taken every 0.35 cm~'. The final 
curve was plotted in percent absorption. 


The photographic region 

Two bands were found photographically at 
0.724 and at 1.04, both apparently of the L 
type. Prints of these do not reproduce well and 
are omitted here, but the frequencies of the 
lines are given in Table V. The lines marked with 
an asterisk were most accurately measured and 
given a spacing of 6.7 cm™'. No interpretation 
of the 0.724 band can satisfactorily be given at 
this time, but the one at 1.0u is believed to be 
a single | band. 


TABLE VI. 
Character 
observed Predicted 
Observed » Identification bands character 
3047.2 A. 
1355.3 v3 
1459.6 a. 
732.3 il 
1019.9 
vet2vs 
4047.9 votre 1, 1, 
4383.8 
6050.0 
13800.0 ? L 


In conclusion we submit Table VI giving the 
centers of the observed bands, and the present 
interpretation of them. The character of these 
bands as observed is compared with that pre- 
dicted by the theory of such molecules which 
have three equivalent particles. 
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PHYSICAL REVIEW 


VOLUME 46 


Energy Levels of the Sm*** Ion as It Exists in Monoclinic 
Crystals of Sm,(SO,) 3° 8H.O 


FRANK H. SpEDDING AND RICHARD S. Bear, Chemical Laboratory, University of California 
(Received October 11, 1934) 


The absorption spectrum and energy levels indicated therein for the samarium ion in crystal- 
line monoclinic Sm,(SO,);-8H,O are shown to be similar to those of SmCl,-6H,O, which is also 
a monoclinic salt. Absorption lines of the sulfate at temperatures between 20°K and room 
temperature are used to prove the existence of several low-lying levels separated from the 
basic one by 160, 188 and 225 cm™'. These levels are of particular interest in view of heat 


capacity and magnetic studies of this salt. 


N previous papers of this series' we have de- 
scribed the absorption spectra of the sama- 
rium ion for monoclinic crystals of SmCl;-6H,O 
and for hexagonal ones of Sm(BrO,);-9H,O, 
and SmI,-9H,O. The chlo- 
ride spectrum was found to be distinctly different 
from those of the hexagonal group in that its 
individual multiplets were spread out somewhat 
more and frequently contained more lines. These 
results were in agreement with those of Spedding 
and Nutting® for the corresponding gadolinium 
compounds, i.e., the extension within multiplets 
and the number of components of each of the 
multiplets of the rare earth ions seem to depend 
more on the positions of surrounding atoms or 
ions in the crystal and Jess on the chemical 
composition of the compound, except insofar as 
the composition afiects the crystal structure. 
This paper describes the spectrum of a second 
monoclinic salt, Sme(SO,);-8H,O. The magnetic 
susceptibility’ and heat capacity‘ of this solid 
have already been measured over a wide temper- 
ature range, and a spectroscopic determination 
of the energy levels will be useful in attempting 
to understand the results of these other investi- 
gations. Here we present only the experimental 
facts. In a later paper we shall introduce a 
discussion of the relation between the physical 
and chemical properties of samarium salts and 
their energy levels. 


Spedding and Bear, Phys. Rev. 42, 58, 76 (1932); 
44, 287 (1933); 46, 308 (1934). 

* Spedding and Nutting, Phys. Rev. 38, 2294 (1931): 
J. Am. Chem. Soc. 55, 496 (1933). 

3 Freed, J. Am. Chem. Soc. 52, 2702 (1930). 

* Ahlberg and Freed, Phys. Rev. 39, 540 (1932). 
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EXPERIMENTAL PART 


Because of the difficulty of obtaining the 
moderately soluble sulfate as large single crystals, 
conglomerate absorption has been used exclu- 
sively. However, the constituent crystals were 
of sufficient size to escape the description 
“powdered” and were of approximately the same 
individual dimensions as the material used in 
preparing the former conglomerates. Groth® 
states that the crystals are monoclinic in 
external symmetry, with axial ratios a:b :c 
= 3.0030 : 1 : 2.0022, 8=118° 16’. 

Two different preparations (however, from 
the same ultimate source) gave identical results. 
One was secured by recrystallization ef sulfate 
of high degree of purity as originally supplied by 
the Wellsbach Company ; the other was obtained 
from the same material after precipitation of the 
samarium as the oxalate, ignition to the oxide, 
solution in HCl, evaporation and reconversion 
to the sulfate by the method described by 
Freed* (double decomposition between H,SO, 
and the chloride in alcoholic solution and 
digestion of the solid sulfate with water). 
Crystals finally used were grown by slow 
evaporation of saturated solutions. 

Even more than in the chloride investigation 
was it difficult to choose the proper thickness of 
absorbing layer. Because of the magnitude of 
the separation between the basic state and the 
excited lower levels, rather thick conglomerates 
are required to bring out the high-temperature 
lines, i.e. those increasing in intensity as the 
temperature is raised. The extreme thickness 


* Groth, Crystallography, Vol. 2, p. 458. « 
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976 F. H. SPEDDING 
in turn offers trouble with increased complexity 
of spectra and makes longer exposures necessary. 
The sulfate has the additional inconvenience of 
showing rather diffuse lines. In obtaining the 
photographs with the sulfate, however, the same 
thicknesses, from 5 to 10 mm, have been found 
most satisfactory. 

Investigation of only the visible region of the 
spectrum has been done. The equipment de- 
scribed before was used to obtain the photo- 


graphs at low temperatures. 


RESULTS 


In Fig. 1 are shown the spectra of the sulfate 
as transmitted at 20°, 78°, 169° and 298°K by a 
conglomerate of 5 mm thickness. In addition is 
shown a 10 mm layer photograph as it appears 
at 78°. Table I gives measurements and roughly 
estimated intensities of the absorption lines and 
bands observed at those temperatures. 

Close examination of the photographs at low 


ta 


u 


AND R.S. BEAR 

temperatures of the thin- and thick-layer spectra 
shows a marked tendency for the thicker layer 
to produce doublets from lines which are ap- 
parently single in the thin-layer case. In some 
instances the lines can be observed to be double 
for both thin and thick layers and then the 
effect is noticed as a widening of the separation 
between the two peaks of intensity. 

This behavior suggests that the double lines 
themselves are not completely resolved. If they 
were composed of very fine lines of varying 
intensities, with the strongest located on the 
inner side (with respect to the doublet) of each 
component, one would understand how the 
widening occurred. Increase in layer thickness 
would enlarge the apparent separation of the 
doublet components by increasing the number 
of fine lines that are able to absorb strongly the 
light. The fine unresolved lines composing the 
doublet components could possibly be results of 
Stark splitting of the resolvable lines by the 
crystal field, but one is more inclined to suppose 


Fic. 1. Conglomerate absorption of Sm_.(SO,),-8H,O. The substances at whose boiling points the spectra were photo- 
graphed are indicated on the right, the corresponding temperatures on the left. H is placed above regions wherein the 
most prominent high-temperature lines are located. 
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that they are due to magnetic interactions 
between samarium ions because of the extreme 
smallness of the effect. This phenomenon is not 
entirely absent in the chloride spectrum but it 
is interesting that this peculiarity should be 
most pronounced with the sulfate, whose spec- 
trum in other ways seems to be that resulting 
from a crystalline field which is even stronger 
than that of the chloride. 

It seems reasonable to assume that the easily 
distinguishable lines composing a multiplet group 
are the ones that reflect the separations between 
levels arising from electric splitting of a highly 
excited level which would be single and de- 
generate in the absence of the crystal field. 
Smaller separations between lines within a 
multiplet may indicate small electric separations 
(of the order of 10 cm~ or less) of the well- 
shielded basic */J;,2 level, while the variability 
in apparent position of individual lines indicates 
that a further, possibly magnetic, perturbation 
affects the levels. This division of the experi- 
mentally observed effects cannot be proved at 
present, however, since both electric and mag- 
netic phenomena would produce qualitatively 
identical results, i.e., with stronger or less homo- 
geneous fields the magnitudes of both types of 
splitting would be greater. If the above explana- 
tion of the behavior of the lines is correct, one 
would expect that greater spectroscopic resolu- 
tion, and perhaps lower temperatures, will reveal 
important facts concerning the nature of energy 
levels in solids. The longer exposures required 
for this procedure (they are already from 30 to 
60 hours in length) would raise extreme diffi- 
culties at 20°K or lower. 

In spite of the trouble in locating the exact 
position of the individual components of the 
doublets, it has been found possible to determine 
pretty definitely the situation of the lower levels 
that are most important spectroscopically. This 
has been done by working with centers of 
doublets and diffuse lines or bands. The results 
for the most clear multiplets (the same ones 
previously diagrammed for the chloride) are 
given in Fig. 2, though many of the less im- 
portant faint low-temperature lines of these 
regions are not included. All of the most intense 
high-temperature lines in these multiplets can be 
explained by the use of the levels 0, 160, 188 


and 225 cm~. Evidence from other more complex 
portions of the spectrum supports these values 
in general. Again both constant energy separations 
and proper temperature behavior of intensity have 
been made the requirements for the location of 
the lines in the diagram. 

There is some evidence at higher temperatures 
for the existence of another level at approxi- 
mately 245 cm= from the basic one. This is 
most apparent at about 23,600 to 23,800 cm™, 
where no other simple explanation of the red- 
most ends of the bands observed at ethylene 
and room temperatures seems possible. The 
faint bands, which extend 300 and more cm 
from the nearest low-temperature line in the 
chloride spectrum and which probably require 
levels of corresponding separations, are not 
visible in the sulfate spectrum. This may be 
related to the general displacement of the sulfate 
levels to higher energy values, a shift perhaps 
just able to eliminate possibility of observing 
transitions from these levels. 

In the foregoing consideration it is assumed 
that the basic levels of the ions of the two salts 
are similar (*H5.) and that the 145, 160, 204 
and 217 levels of the chloride are related, 
respectively, to the 160, 188, 225 and 245 ones 
of the sulfate. In this connection it is interesting 
that though the chloride lines arising from the 
levels 160 and 217 seem to be most conspicuous, 
in the sulfate the 160 and 225 lines are most 
noticeable. This suggests a reversal of the rela- 
tion just given within the doublets of levels. 

The multiplet on the long wave-length end of 
the spectrum at 17,000 to 17,300 cm™ is one 
that offers some rather puzzling features. In 
the chloride, ethylsulfate and bromate spectra 
nothing of certainty could be gathered from it, 
hence it could not be used to determine the 
energy levels. The corresponding sulfate lines, 
however, do seem to indicate the presence of 
another level by the existence of two high- 
temperature satellites 78 cm=' from their re- 
spective low-temperature companions. It is not 
justifiable to base the existence of a level, whose 
Boltzmann population should allow it to be 
made more evident elsewhere in the spectrum, 
upon such isolated cases. Though the evidence 
at hand is inconclusive, it would not be surprising 
to have such a level be demonstrated éventually. 


| 
| 
| 
| 
| 


982 F. H. SPEDDING AND R. S. BEAR 


22600 
ig 
q 
am 
| WH 
228 


Fic. 2. Energy levels of Sm2(SO,);-8H:O. Shaded areas represent bands. Figures 


in parentheses denote 


band edges measured at ethylene temperature. All other figures are from plates taken at 78°K. Levels 
which are actually probably double are indicated by “2.” 


It may be anticipated that the bromide spectrum 
also has a possible level at 58 cm™', evidence 
for which is based on similar scanty evidence in 
this region, beside the proper “monoclinic” 
levels somewhat higher. 

Table II shows how well the lines agree with 
the level diagram assumed. In many other places 
on the plates the existence of the new lines 
brought into evidence at low temperatures by 
the lessened thermal agitation of the crystal and 
the higher population in the basic state makes 
it impossible to pick out satisfactorily the 
location of the high-temperature satellites. How- 
ever, at ethylene temperature (169°K), when the 
fainter low-temperature lines have been ‘‘washed 
out,” one is often able to locate additional 
confirmation of the postulated levels. 

The most important relations between the 
sulfate and chloride spectra in other types of 
behavior are the following. The sulfate shows 
practically the same amount of displacement in 
the positions of lines upon changing the temper- 
ature from 298° to 20°K. At any given tempera- 


ture the sulfate multiplets are shifted to shorter 
wave-lengths from the corresponding ones in the 
chloride spectrum, though only by about 10 to 
20 cm~!. On the other hand it will be recalled 
that the bromate acted in a different manner, 
manifesting slightly greater shifts with temper- 
ature change and possessing multiplets located 
50 to 100 cm~! toward longer wave-lengths than 


TABLE II. Prominent low-temperature lines and their high- 
temperature satellites (78°K). 


L.-t. line Satellite A Satellite B Satellite C 
v(cm™) Avg v(cm) Avg’ v(cm™") Ave 


17702 17667 224 
17732187 17692 227 
18949 18787 162 18762 187 18722 227 
20041 19882 159 Too diffuse 19815 226 
20056 19895 161 19868 188 19831 225 
20130 19968 163 19942 188 

20522 20361 =161 

20621 20456 165 

22190 22030 160 21966 224 
22239 22074 22052 = 187 22014 225 
22255 22094 161 

23975 23814 = 161 23789 =186 23750 


17891 17732 159 
17919 17758 = 161 
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those of the chloride. The ethylsulfate of sa- 
marium resembled the bromate in these respects. 


CONCLUSIONS 


Again, this time in the sulfate spectrum, is 
proved the existence of more than three lower 
levels which are of importance in determining 
the absorption spectrum of the samarium ion. 
It seems impossible that crystal splitting of the 
basic °/7s,2 level could account for them all, 
since only three new levels (J+4) can be 
expected when electric perturbations remove 
the degeneracy of such a state. A possible 


explanation of these levels has been advanced 
by one of us.*® 

Magnetic interactions between samarium ions 
could increase the number of levels resulting 
from the *//;,. one, but it is unlikely that these 
could be of sufficient magnitude since the 
magnetic moment of a samarium ion js so small. 
We have indicated that magnetic effects are 
probably responsible for the diffuseness and 
variability in appearance of certain of the 
absorption lines, an effect so small as to be 
unresolvable on the present plates. 


* Spedding, J. Am. Chem. Soc, 54, 2593 (1932). 
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Elastic Electron Scattering in Potassium 


J. Howarp McMILLEN,* Palmer Physical Laboratory, Princeton, New Jersey 
(Received July 2, 1934) 


A method was developed for the measurement of slow 
electron scattering in vapors. Angular distribution curves 
from 25° to 160° were obtained for 5, 10, 15, 25, 50, 75, 
100 and 150 volt electrons elastically scattered by potas- 
sium atoms. All curves contained a large peak near 100°. 
The scattering curves fell off rapidly for angles up to about 
60° and rose rapidly for angles greater than about 130°. 
The position of the peak was found to shift to larger and 
then smaller angles as the energy of the colliding electrons 


INTRODUCTION 


HE angular distribution curves for electrons 
scattered elastically by gaseous atoms 
generally contain a set of maxima and minima. 
These are referred to as diffraction patterns and 
are attributed to the diffraction of the associated 
de Broglie waves in their passage through the 
field of the atom." 

If exchange effects are neglected the nature of 
these patterns is seen to depend almost entirely 
upon the electrons’ distance of closest approach 
as taken from the classical theory and the 
magnitude of the field at that distance. Con- 
siderable progress has been made in our under- 


* National Research Fellow. 
1 For general discussion of theory, see Mott and Massey, 
Theory of Atomic Collisions, Oxford Press (1933). 


was decreased. Henneberg's calculations which employ the 
Thomas-Fermi atomic field were found to give scattering 
curves whose general shapes agree quite well with the 
experimental curves. Best agreement occurred for the 5 
and 50 volt electrons. When the potassium and argon 
experimental curves were compared it was found that the 
curves were quite similar for electron energies of 50 volts 
or more. For slower electrons the curves became more 
dissimilar. 


standing of the relations between these factors 
by the comparison of diffraction patterns for 
various gases and for various energies of the 
scattered electrons. 

So far scattering measurements have been 
made on several important groups of atoms, 
namely, the inert gases,” * the halogens,‘ and 
three elements from the third column of the 
periodic table, zinc,’ cadmium,* and mercury.” * 
The longest unbroken sequence of increasing 


*F. L. Arnot, Proc. Roy. Soc. A133, 615 (1931). 
asst McMillen and Webb, Phys. Rev. 41, 154 

*F. L. Arnot, Proc. Roy. Soc. A144, 360 (1934). 

5E. C. Childs and H. S. W. Massey, Proc. Roy. Soc. 
Al41, 473 (1933). 

*E. C. Childs and H. S. W. Massey, Proc. Roy. Soc. 
A142, 509 (1933). 

7 F. L. Arnot, Proc. Roy. Soc. A140, 334 (1933). 
*E. B. Jordan and R. B. Brode, Phys. Rev. 43, 112 
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atomic number is that made up of phosphorus, 
sulphur, chlorine and argon.° 

It was deemed advisable to exten the data 
pertaining to electron diffraction patterns by 
obtaining curves for one of the alkalis. Of this 
group potassium was considered the most im- 
portant, as it extends the phosphorus-argon 
series mentioned above. It was also hoped that 
because of its hydrogen-like structure the experi- 
mental data would be of considerable theoretical 


interest. 
APPARATUS 


Electron scattering in a vapor was first 
accomplished by Childs and Massey.* They 
concentrated the vapor at the scattering center 
by utilizing an undefined jet. A similar arrange- 
ment was attempted by the author but aban- 
doned because of the following difficulties. It 
was found first of all that the background 
scattering (scattering current observed with the 
oven cooled) formed too large a portion of the 
scattering from the vapor.'® Increasing the 
effective pressure within the jet by raising the 
temperature of the oven so reduced the life of 
the charge as to prohibit any improvement in 
this direction. Consideration of the possibility 
of lengthening the collision path was abandoned 
because of the non-uniformity of pressure 
throughout the jet included in this path. Over 
the collision path the pressure within the jet 
varies as cos* ¢ where ¢ is the angular spread of 
the jet. This variation in pressure is negligible 
as long as the collision path is very short or the 
oven far away. Both these latter requirements 
seriously reduce the magnitude of the scattered 
current. 

These difficulties were overcome by con- 
structing a collision chamber that was practically 
molecule-tight for low pressures and which was 
capable of being raised to about 150°C, a 
temperature which gives a pressure of about 
0.001 mm of mercury. A small slit in the chamber 
permitted the passage of the scattered electrons 
to the collector and at the same time accounted 


*C. B. O. Mohr and F. H. Nicoll, Proc. Roy. Soc. A138, 


469 (1932). 
1 This may be the surface ionization effect observed by 


Childs and Massey, reference 5, when they attempted to 


scatter electrons from potassium in the apparatus which 
they subsequently used for cadmium and zinc. 
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for very little loss of the vapor. Within this 
chamber was rotated the electron gun. This was 
accomplished by suspending the gun from the 
lid of the chamber and then rotating the lid 
while the lid’s own weight kept it in place. 
With this method as high a pressure as was 
desired could be maintained and with.a much 
less expenditure of energy than in the jet method. 

Details of the apparatus are shown in Fig. 1. 
A large brass base plate P, a cylindrical wall W 


Fic. 1. Apparatus. 


and a disk K formed a chamber in which was 
placed the oven-collision chamber A and the 
electron energy filter F. Electrons from a 
tungsten filament leave the electron gun G and 
enter the collision chamber A. Those scattered 
at a particular angle @ pass out of the chamber 
through S to the electron energy filter F and 
from there to the collector C. 

The collision chamber was a thick walled 
chamber, made of nichrome IV alloy. This 
alloy is non-magnetic and gave no evidence of 
reacting with the potassium. The lower part of 
A held the charge of potassium and fitted into 
a lavite furnace B which was wound with heating 
wire. The chamber was supported in a glass 
tube which was waxed into a base plate P. The 
chamber was well insulated from losses by heat 
conduction by using small wire stays and glass 
as supports. To reduce the radiation losses the 
chamber was polished. Water jackets served to 
cool the rest of the apparatus in the proximity 
of the oven. 

The electron gun was made of thin wall 
nichrome IV and fastened rigidly to the top of 
the chamber A. By means of a ground glass 
joint the top of A could be rotated. 

The electron energy filter separated the elec- 
trons of various energies through the use of a 
radial electrostatic field. The details of this 
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method have been discussed in earlier reports." ™ 
In this experiment the parts were contained in 
a brass box. The box was supported on a brass 
tube which ran out through the base plate. The 
filter was electrically insulated from the rest of 
the apparatus and it was possible to lower or 
increase the energy of the electrons before 
entering the filter. The resolution of the electron 
energy filter was such as to filter out the electrons 
suffering the first inelastic loss (1.6 volts) if the 
energy of the elastically scattered electrons was 
50 volts or less. Hence in studying the scattering 
of electrons of energy greater than 50 volts 
they were retarded to 50 volts energy, by 
appropriately adjusting the potential of F, before 
the electrons entered the electron energy filter. 

The pressure in the oven was determined from 
thermocouple readings of its temperature. The 
pressure was kept low enough to insure single 
collisions. At the temperatures used the vapor 
was 99.95 percent monatomic.'* 

Readings were taken before and after the oven 
had been heated to measure the scattering from 
the walls of the apparatus. These readings were 


( a L. Hughes and V. Rojansky, Phys. Rev. 34, 284 
1929). 
( A. L. Hughes and J. H. McMillen, Phys. Rev. 34, 291 
1929). 
143F, W. Loomis and R. E. Nusbaum, Phys. Rev. 39, 
89 (1932). 


120 160° 
Scattering Angle 
Fic. 2. Angular distribution curves for electrons of various energies: continuous line for potassium, broken line for argon. 


negligible in all cases except that in which 5 volts 
electrons were used. For this case the back- 
ground scattering had to be subtracted. Every 
set of readings was checked by repeating on 
both sides of the zero angle. 

The following dimensions may be of interest : 
Slit S, 8X4 mm; slits S,, S: and S;, 1X10 mm; 
diameter of holes in gun, 1 mm; distance between 
S, and S:, 10 mm; diameter of deflection plates 
M and N, 24 mm and 30 mm; cross section of 
electron gun, 3 mm by 11 mm. 


EXPERIMENTAL RESULTS 


Elastic scattering data was obtained for 5, 10, 
15, 25, 50, 75, 100 and 150 volt electrons over 
an angular range extending from 25° to 160°. 
(See Fig. 2.) 

All curves show a peak or maximum in the 
vicinity of 100°. This peak shifts to larger angles 
and then retreats to smaller angles as one goes 
from higher colliding electron energies to lower 
ones. The peaks for electrons of energies less 
than 100 volts are. symmetrical about their 
center and their positions quite easily ascer- 
tained. In the 150 volt curve, however, the peak 
is not well defined and its position is somewhat 
difficult to determine. The position of the peak 
is given as a function of colliding electron energy 


\ | 
vas \ 
lid \\ \ 
ice. \ / \ 25 7s | 150 
W / | 
\\ / 
y f 50 100 
the 
a 
ind 
red 
ber | 
ind 
led 
his 
of 
of 
ito 
ing 
ass 
‘he 
eat 
ass 
the 
to 
ity 
all 
of | 
ass 
ec- 
fa 
his 


986 J. 
m 


s is so 100 150 
Electron Energy (votts) 


Fic. 3. Showing how the position of the principal 
maximum varies with the scattered electron energy: 
circles are experimental for potassium, rectangles theo- 
retical, ical dashes are experimental! for argon. 


in Fig. 3. One notes that the points lie on a 
smooth curve. The position of the peaks could 
be determined in most cases to within 23°. At 
larger angles beyond the maximum the curves 
rise rapidly, indicating large probability for 
backward scattering. On the 10, 15 and 25 volt 
curves there is an indication of a second peak 
somewhere in the neighborhood of 40°. Its exact 
position is very difficult to determine, however, 
as it is superimposed on a very steep part of 
the curve. At small angles the scattering falls 
off less rapidly for slower electrons, as has been 
found to be the case in the investigations of 


other gases. 
DISCUSSION 


The wave theory of elastic electron scattering 
is embodied in the following expression :"* 


(21-+-1)(2I’ +1)P,(cos 8) 
Pr(cos @) sin sin yr cos (yi— 


where dw is the increment of solid angle P, is 
the Legendre polynomial of the /th degree, / and 
l’ the orbital angular momentum quantum 
numbers, and y; the phase shifts brought about 
by the scattering field of the atom. Both ex- 
change effects and the possibility of polarization 
have been neglected in arriving at this expression. 
Neither become important until colliding elec- 
tron energies approach magnitudes comparable 
with those of the ionization potentials. Just how 
necessary it is to take into account these effects 


is still very uncertain. 


14 P. M. Morse, Rev. Mod. Phys. 4, 577 (1932). 
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To evaluate Eq. (1) it is necessary to compute 
the phase shifts from the atom’s potential field. 
Faxen and Holtsmark'® have carried out an 
approximate computation for argon and krypton. 

Henneberg'’® has calculated the phase shifts 
for mercury, using the Wentzel-Kramers-Bril- 
louin method of approximation and the Thomas- 
Fermi field. In this method the phase shifts are 
obtained by numerically integrating (in atomic 
units) the following expression : 


v= f 
R 
1(1+1) 
Re 


The limit R is the positive root of the first 
integrand and corresponds to the distance of 
closest approach on the particle theory of scat- 
tering. The limit Ry is the positive root of the 
second integrand and in the particle theory is 
known as the impact parameter. Phase shifts for 
mercury are given by Henneberg as a function 
of J for various values of k°. Now a very im- 
portant feature of Henneberg’s theory lies in the 
fact that these phase shifts as calculated for 
mercury are applicable to any other scattering 
atom. To obtain a set of phase shifts for an 


| dr. (2) 


2 
Valves of 2 


Fic. 4. Showing potassium phase shifts and their sines for 
various / values and electron energies. 


16H. Faxen and J. Holtsmark, Zeits. f. Physik 45, 307 


(1927). 
16 W. Henneberg, Zeits. f. Physik 83, 555 (1933). 
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Fic. 5. Continuous line indicates experimental scattering 
curve, broken lines refer to the calculated curves of Henne- 


berg. 


atom of atomic number Z one first defines a 
variable g as: 


q= (2'/Z)!, 


where Z’ represents the atomic number of 
mercury. The new phase shift is then given by: 


yz(k)=q" (qk) 
and the new / values by: 
—}. 


Using these relations, it was possible to obtain a 
set of phase shifts for potassium. A few of these 
are plotted in Fig. 4. 

Now Henneberg also pointed out that to 
duplicate the main features of the slow electron 
scattering curves it was necessary to employ 
only the first few terms in expansion of Eq. (1). 
It turns out that for the higher terms the values 
of (2/+1)(sin y;) change so slowly with / that 
the sum of all these terms gives a curve which 
falls off monotonically. Hence the characteristics 
of the curves for slow electrons other than those 
of falling off monotonically will be given by the 
first few Legendre coefficients whose sin y values 
change rapidly with 7. These large changes in 
the values of sin y come about for values of y; 
greater than x/2 and from Fig. 4 we see that for 
potassium they include values of / equal to and 
less than 2. 


We have used these phase shift values and 
summed the six terms which involve Po, P;, Ps 
and their combinations. These curves are com- 
pared with the experimental ones in Fig. 5. 
One observes that the agreement is surprisingly 
good. The main characteristics of the experi- 
mental curves, namely, a maximum near 100°, 
are reproduced. The shift of the experimental 
maximum to large and then smaller angles as 
one reduces the energy of the colliding electron 
is also in evidence in the calculated curves. In 
Fig. 3 the position of the calculated maximum 
does not always agree any too well with the 
experimental, except perhaps in the 5 and 50 
volt curves. Another feature of the experimental 
curves that is entirely lacking in the theoretical 
curves is that of the second maximum in the 
vicinity of 40° which is present in the 10, 15 and 
25 volt curves. It is quite evident that Henne- 
berg’s calculations gives the correct form of the 
scattering curves and that the curves shift with 
colliding electron energies in the proper manner. 
However, these same calculated curves would 
fit the experimental argon curves almost as well. 
Whether this lack of precision is due to the 
crudeness of the Thomas-Fermi field employed 
or whether it is due to the approximation of the 
method is difficult to say. Deviation between 
theory and experiment for the slower electron 
energies can of course be safely attributed to 
the neglect of exchange considerations. 

As argon lies next to potassium in the periodic 
table, it is of interest to compare curves of 
argon!’ with those of potassium. This can be 
done by inspecting Fig. 2. In the first place one 
notices that their diffraction patterns are of 
the same type, namely, that of a maximum in the 
vicinity of 100°. Moreover, one notices that the 
50, 100 and 150 volt curves agree quite well, 
the only discrepancy being that the position of 
the maximum is at slightly smaller angles in the 
case of argon. For energies less than 50 volts, 
however, the curves become more discordant. 

This tendency of scattering curves for faster 
electrons to be similar, and slower electrons to 


17 The 25, 50, 100 and 150 volt curves for argon were 
taken from unpublished data furnished through the 
kindness of A. L. Hughes and G. M. Webb. The 5 and 10 
volt curves are those of Ramsauer and Kollath, Ann. d. 
Physik 5, 529 (1932). The 15 volt curve is that of Nicoll 
and Mohr, Proc. Roy. Soc. A142, 320 (1933). 
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be dissimilar, when curves of atoms differing 
very little in atomic number are compared, was 
first observed by Mohr and Nicoll by comparing 
molecular H,S and PH; curves with the atomic 
argon curve.* Their comparison was based on 
the assumption that the sulphur and phosphorus 
atoms were largely responsible for the shape of 
the molecular scattering curve. That we have 
verified this effect with atomic potassium bears 
out their assumption. The explanation for this 
effect is simply that fast electrons are scattered 
in the inner shells of the atom where the fields 
are very much alike and hence give rise to 
similar scattering curves. Slow electrons, how- 
ever, are scattered in the valence electron shells 
and there encounter very much different fields, 
thus giving scattering curves of different pat- 
terns. 

To show just what part of the atom is most 
influential in determining the form of the 
scattering curve, we have superimposed on the 
density distribution curve of argon'® the two 
distances R and Ro. See Fig. 6. On the particle 
theory it is between these two distances that 
the particle experiences the greatest deviation as 
it passes through the central field of the atom. 
On the wave theory one finds that the magnitude 
of the phase shift as computed by Eq. (2) is 
largely accounted for by the integration between 
these two distances R and Rp». Since our curves 
are predominantly of the P: type one is inter- 
ested in these radial distances for the case in 
which the / value is 2. We see from the figure 
that it is electrons of energy less than 50 volts 
that are influenced by the composition of 
potassium’s valence electron shell. One also 


" Hartree’s distribution curve for argon was used, since 
the one for potassium was not available. It is also im- 
portant to bear in mind that the R and R» values were 
calculated with the Thomas-Fermi field of potassium. 
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Electron Densit 


° 4 


2 3 
Radia! Distance (atomic units) 


Fic. 6. Electron density distribution for argon's L and 
M shell with arrows indicating the R and Rp distances for 
various colliding electron energies and / values. 


notes that the part played by this valence 
electron shell increases for the slower electrons 
just as the difference between the experimental 
curves for these two atoms increases for slower 
electrons. 

Here again, however, one has an alternative 
explanation of the relation between the curves 
for these two gases, namely, the effect of electron 
exchange. For it is just in this energy range, 
that is, for energies near that of the ionization 
potentials, that one expects the exchange effects 
to enter. Unfortunately there is nothing in the 
experimental curves to suggest that one or the 
other of these explanations is the more proper. 

The author wishes to express his appreciation 
of the facilities offered him while carrying on 
this work in the Palmer Physical Laboratory, 
and to thank Professor R. Ladenburg for the 
interest shown in this work. 
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The Spectrum and Ionization Potential of Radium 


Henry Norris Russet, Carnegie Institution of Washington, Mount Wilson Observatory 
(Received October 15, 1934) 


A new interpretation of the connection between the two groups of terms in Ra I found by 
Rasmussen is suggested and confirmed by new observations which he has courteously com- 
municated. Certain term-designations are altered. The first ionization potential is 5.252 volts. 


ASMUSSEN, in his admirable analysis of 
the arc spectrum of radium' has con- 
clusively identified two groups of lines—one 
resulting from transitions from *D and '‘D terms 
to higher odd levels and the other from *P° 
and 'P* to higher even levels and the ground- 
term ‘So. Absolute term-values for the first 
group are well determined from a series of four 
3F°® terms. The only series available in the other 
consists of two *S terms, and gives but a rough 
limit. Connection between the two groups should 
occur through the ultraviolet line 3101.80 
which must arise from the 'S level. Rasmussen 
assigns this as 1'S,;—3'P,.2, The combination 
3'D,—3'P, should lie in the red and he adopts 


the strongest unclassified line at 6336.90. All. 


terms are thus referred to an accurately known 
limit, giving 1'S)=41,949.19 and 5.176 volts for 
the ionization potential: but this arrangement 
is not confirmed by any other combinations. 
An independent determination of the term-values 
for the P® terms and those which combine with 
them may be made by combining the low 'S and 
the two higher *S terms in a single series. This 
apparently crazy method—first suggested to the 
writer some time ago by Professor A. G. Shen- 
stone in connection with another spectrum— 
works very well, as is illustrated by Table I 
which gives the Rydberg denominators n* for 
the other alkaline earths. 

The progression of these values is very regular 
and leads to the estimates for Ra I given in the 
corresponding line. The estimated tolerances 
represent the probable limit of error provided 
that the progression remains regular. The “‘ob- 
served’’ values are discussed below. With Ras- 
mussen’s value 'S)=41,950 the values of An* 


1 Rasmussen, Zeits. f. Physik 87, 607 (1934). 

2 These serial numbers are “‘Laufzahlen” beginning with 
1, 2, 3, 4 for the S, P, D and F series. They are retained 
here for comparison with Rasmussen's tables. 


TABLE |. Rydberg denominators for alkaline earths. 


Term 23S, 39S, 
Element n* An* n* An* n* 
Be 1.208 | 0.971 | 2.179 | 1.026 | 3.205 
Mg 1.334 | 0.981 | 2.315 | 1.031 | 3.346 
Ca 1.492 | 0.993 | 2.485 | 1.040 | 3.525 
Sr 1.546 | 1.003 | 2.549 | 1.044 | 3.593 
Ba 1.616 | 1.014 | 2.630 | 1.045 | 3.675 
Ra(est.) 1.025 + 0.005 | 1.050 + 0,003 
Ra(obs.) 1.605 | 1,028 2.633 | 1.047 | 3.680 


come out 1.070 and 1.145 (*S—8S). 
Increasing the assumed limit by 500 reduces 
these to 1.036 and 1.066; by 200 more, to 1.024 
and 1.036. The two estimated values of An* for 
Ra I therefore require the additions of 683480 
and 606+20 to Rasmussen's limit,—or, in the 
mean, 630+30. With the new value of ‘So, the 
term which combines with it to give \3101 
comes out 10,350+ 30. 

Now Rasmussen has already identified a level 
dp*D,° at 10,347.38. The corresponding level in 
Ba I combines strongly with ‘So, and it appears 
very probable that the radium line is homologous. 
Upon communicating this suggestion to Dr. 
Rasmussen the writer received a cordial letter 
stating he had found upon his plates another 
ultraviolet line which confirms the new arrange- 
ment; 2955.65, intensity 2. As the further 
investigation upon this spectrum, which he plans, 
may take some time, he generously suggested 
that this evidence should be presented in the 
present paper. The new combination, together 
with some others whose interpretation has been 
altered, are exhibited in Table II. The intensities 
are given before the wave numbers, and the 
differences of the observed and computed values 
below. It is evident that the true connection 
between the two parts of the spectrum has been 
found. Retaining Rasmussen's values for the 
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990 
TaB_e II. Term values in Ra I. 
3'D, 3°D; 34D; 
25495.95 27870.00 28583.38 28861.50 42577.35 

dp*D,° 10347.38 (3) 15148.55 (8) 18236.03 (10) 18514.07 (6) 32230.02 

—0.02 +0.03 —0.05 +0.05 
dp*P? 8194.44 (5) 17301.40 (6) 19675.58 (4) 20389.00 (3) 20667.06 
—0.11 +0.02 +0.06 0.00 

3pP,° 8753.65 (5) 19829.62 (5) 20107.96 (2) 33823.65 

—0.11 +0.11 —0.05 
8794.94 (4) 20066.56 
(0.00) 


first set (3*D, etc.) those of his second set (given 
in his Table IV) should be increased by 627.66. 
The corresponding values of m* for the S terms 
are given in the last line of Table I. 

The ionization potential corresponding to the 
new data is 5.252 volts. 

The adoption of the term dp*P,°, as here 
given, accounts for the unclassified line at 
44982 (*D,—*P,) and removes from the list 
Rasmussen's doubtful level at 5431.68 combining 
with 2°P°. His levels at 10,702.91 and 11,014.06 
may then be interpreted as 3°P,° and 3'P,°. 
The isolated line at \5755.45 (v= 17,370.03) then 
gives 3*P,° = 11,491.47. This interpretation differs 
from Rasmussen's mainly in the choice between 
alternatives which he mentioned as possible. 
It puts dp*P® higher than dp*D°, as in the case 
of Ca, Sr and Ba, and 3°P° between dp*F° and 
dp*D®, as is true in Ca and Sr (but not in Ba, 
where the limits of the normal and displaced 
series are closer together). Apart from some 
effects of this difference, the similarity of the 
spectra of the alkaline earths is remarkably 
detailed. For example, dp'D® is superposed on 
dp*F* in all four spectra, and the separation 
dp(*P,°—*P,°) is abnormally small in every case. 

Some revision of the assignments of the high 
level even terms appears also to be desirable. 
Table III gives the new and old term-values and 


TABLE III. Term values and designations. 


Old New 

28,871.25 23P,° 29,498.91 
25P, 27,950.31 23P,° 28,577.97 
25P; 25,261.15 2*P,° 25,888.81 
2'P, 21,233.98 2'P,° 21,861.64 
dd'*P, 9948.87 4°D, 10,576.53 
4°D, 9956.28 4°D, 10,583.94 
4°D; 9752.41 45D; 10,380.07 
45D, 10,701.08 *P; 11,328.74, 
9008.56 9636.22 

9734.85 10,362.51 


designations. The slight inversion of 4*D,; appears 
preferable to the very large interval demanded 
by the old assignment. The outstanding levels 
are attributed to the p*, rather than the @, 
configuration, since the former is lower in Ca 
and Sr (which resemble Ra more closely than 
Ba does). Which of the last two levels is *P2 is 
doubtful: the present assignment assumes that 
the difference *P,—*P, is likely to be large. 
Further identifications must be postponed till 
the fainter lines have been more fully observed. 

In conclusion, the writer desires to express his 
hearty thanks to Dr. Rasmussen for his generous 
communication of the new data and his per- 
mission to refer to them here. 
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The Effect of Hydrogen on the Afterglow in Mercury Vapor 


M. M. Mann* AND WALTER M. Niecsen, Duke University 
(Received July 6, 1934) 


Observations have been made on the change in intensity 
of the radiation in the mercury afterglow when hydrogen 
was introduced into the luminous vapor stream. It has 
been found that with no potential applied to a grid 
embedded in the vapor stream, the introduction of hy- 
drogen to a pressure of 4 10-* mm of Hg (at the point of 
observation) decreases slightly but definitely the intensity 
of the afterglow. With higher hydrogen pressures but with 
other conditions the same the lines increase in density. 


With grid potentials of from 3.5 to 4 volts (positive with 
respect to the anode) additional excitation takes place in 
the Vicinity of the grid and the introduction of hydrogen 
to a pressure of 4+ 10~* mm of Hg brings about a marked 
reduction in the intensity of the lines emitted in the 
afterglow. It is suggested that the intensity of the lines in 
the mercury afterglow when excited under these conditions 
is in part dependent upon the concentration of atoms in 
the 6 states. 


T is generally assumed that the nature of the 

spectrum emitted in the afterglow of mercury 
vapor under ordinary conditions of excitation, 
pressure etc., is typical of recombination spectra.' 
Such spectra are characterized by an abnormally 
high intensity of the lines originating in high 
energy levels compared with lines having their 
origin in lower energy states. The work of 
Miss Hayner! on mercury vapor and of Kenty* 
on the afterglow in argon indicate that recombin- 
ation between positive ions and electrons is of 
greater likelihood when the temperature of the 
electrons is small. The variation of recombina- 
tion light as a function of electron energy in 
mercury vapor as observed by Miss Hayner, 
and particularly by Webb and Sinclair, shows 
a much more rapid decrease with increase in 
electron energy than can be accounted for on 
the generally accepted views of the recom- 
bination process. From their work Webb and 
Sinclair concluded that their results could be 
accounted for by assuming that recombination 
takes place in two stages, (a) capture of an 
electron into a high energy level, and (b) transi- 
tions to lower levels with the emission of series 
lines. In order to account for the rapid decrease 
in intensity of the visually observed afterglow 
they assumed that such highly excited atoms 
might suffer a collision with a relatively high 
energy electron and be re-ionized, and from their 
results attempted to evaluate the critical energy 


* Based in part on a thesis submitted by the junior 
author for the M.A. degree at Duke University. 

‘Lord Rayleigh, Proc. Roy. Soc. Al08, 262 (1925); 
All2, 14 (1926); Miss Hayner, Zeits. f. Physik 35, 365 
(1926); Webb and Sinclair, Phys. Rev. 37, 182 (1931). 

? Kenty, Phys. Rev. 32, 624 (1928). 


necessary to carry through such re-ionization 
processes. 

We had initially planned to study the distri- 
bution of energy among the lines of the Balmer 
series of hydrogen recombination spectra and to 
compare it with the theoretical expectation from 
a knowledge of the recombination coefficients* 
and transition probabilities.‘ The observations 
reported in this note are an outgrowth of some 
observations made in an attempt to observe 
such spectra in the luminous stream diffusing 
from an arc struck in a stream of hydrogen and 
mercury vapor.’ In some of the initial experi- 
ments with no added hydrogen in the arc region 
or vapor stream, the intensities of some of the 
lines varied under different conditions (to be 
described later) in such a way as to suggest that 
processes other than recombination were play- 
ing an effective rdle. There are reasons® for be- 
lieving that the concentrations of excited atoms 
in the lowest *Po. states of mercury are under 
certain conditions intimately connected with the 
character of the spectral radiation in the mercury 
discharge, and so it seemed that the effect on the 
afterglow of hydrogen introduced into the lumi- 
nous stream might be of interest. 


EXPERIMENTAL 


The tube used was similar to and of approxi- 
mately the same dimension as that of Webb and 


* Stueckelberg and Morse, Phys. Rev. 36, 16 (1930). 

* Francis G. Slack, Phys. Rev. 31, 527 (1928). 

*The low concentration of undissociated hydrogen 
ey ey present in the arc region of the proton source 

Lamar and Luhr (Phys. Rev. 44, 941 (1933)); 46, 87 
(1934) would seem to provide a more satisfactory means 
of exciting such spectra and efforts on the initial problem 
are now being made in this direction by one of us. 


* Kessel, Zeits. f. Physik 71, 614 (1931). * 
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Sinclair.! Mercury vapor from a reservoir diffused 
upward in a vertical tube and through the arc 
region between an oxide-coated filament and a 
cylindrical ancde (placed upstream from the 
filament). Positive ions, electrons, and excited 
atoms were swept out of this arc region and into 
and along an approximately horizontal side tube 
of 4 cm diameter. The tube was mounted in a 
furnace and heat was applied in such a way that 
condensation was prevented at all points except 
at the end of the side tube outside the furnace, 
from which the condensed mercury was returned 
to the mercury reservoir. A grid of No. 20 P. 
and §. wire with 4 mm spacing was mounted in 
the horizontal tube at approximately 1 cm from 
the junction of vertical and horizontal tubes. 
The diameter of the circular grid was slightly less 


2537 


(a) 
(b) 
(c) 
(d) 
(e) 


os 2967 


(f 


(g) 
(h) 


Fic. 1. Typical spectra of mercury afterglow. Arc 
voltage 5 volts, and arc current 0.5 ampere. Other con- 
ditions as indicated. 

Exposures @ to ¢ inclusive showing effect of positive grid 
potential on afterglow with no foreign gas present. 


Grid voltage 0 Liv 2.2v 3.3 4.4v 
Electron Temperature 756°K 760°K 1344°K 3452°K 4681°K 
Exposures f, g, 4 showing effect of hydrogen with acceler- 
ating voltage of 3.6 volts on grid. 
Exposure f 
Hydrogen pressure in mm Hg 0 
Exposures #, j showing effect of hydrogen with accelerating 
voltage on grid of 4.0 volts. 


4X10™ mm 16x10? 


Exposure i j 
Hydrogen pressure in mm Hg 0 &§x10-" 


Exposure &, Mercury Spectrum. 
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than the inside diameter of the tubing. Immedi- 
ately behind this grid was a quartz observation 
window attached by means of a graded seal. 
A second window was attached 10 cm down the 
horizontal tube from the first window. Probes 
were placed in the center cf the tube opposite 
the observation windows. Hydrogen was intro- 
duced at a point approximately midway between 
the two observation windows. The approximate 
pressure of the gas in the region opposite the 
second observation window was measured on a 
McLeod gauge. The pressure of mercury vapor 
at the same point was fixed as to its order of 
magnitude by a method similar to that used by 
Webb and Sinclair. 

The electron temperatures and concentrations 
at the point of observation were measured in the 
usual manner’ immediately after the corre- 
sponding exposure had been made. The use of 
probes in the measurement of low electron tem- 
peratures was found to be difficult due to changes 
in its contact difference in potential, probably 
brought about by changes in the gas layer on the 
probe. Webb and Sinclair reported the same 
difficulty. We found that fairly reproducible 
measurements could be made by flashing the 
probe to a dull red temperature before each 
series of measurements. A few photographs at 
the second window (downstream window) are 
reproduced in Fig. 1 under conditions given in 
their captions. The exposure time in each case 
was 15 minutes. 


RESULTS 


In Fig. 1 a, 5, c, d and e, we have reproduced 
tvpical spectrograms obtained with no foreign 
gas present in the vapor stream. For this group 
of exposures the mercury vapor pressures was 
0.09 mm of mercury. The electron concentration 
decreased from 1.610" per cc in exposure a to 
6.0 10° per cc in exposure d. The variation in 
positive ion concentration was probably of a 
comparable magnitude. The electron concentra- 
tion in exposure ¢ was 6.5 X 10° per cc. The den- 
sities of the sharp triplet and other lines of 
longer wave-length are probably not recorded 
correctly in this group of exposures because of 


? Langmuir and Compton, Rev. Mod. Phys. 3, 191 
(1931). 
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scattered light from the arc region in the early 
experiments. 

In Table I are given values of the densities as 
measured with a’ Moll microphotometer of a 


TaBLe |. Results of microphotometric measurements of 
densities on typical spectrograms with grid floating. 


Wave- Densities at different hydrogen pressures 
length p=Omm 1.6X10-* 3.6X10-* 6.4x10-* 


2537A 0.140 0.096 0.072 0.057 0.045 
2967 071 -053 083 097 -107 
4358 071 056 071 096 lll 


few of the lines of the afterglow spectrum on 
spectrograms taken at the second window with 
grid floating for various pressures of hydrogen. 
These results are given in this form because of 
the reproduction difficulties connected with 
small differences in density. The mercury vapor 
pressure was fixed at 0.09 mm of mercury. The 
electron temperature and electron concentration 
with no hydrogen present were approximately 
the same as for exposure a above. No attempt 
was made to measure electron temperatures and 
concentrations with the gas flowing. One would 
not expect the small pressures of hydrogen to be 
effective in directly preventing the diffusion of 
positive ions and electrons down the stream. 

Exposures f, g and h show the effect of hydro- 
gen on the afterglow with a higher potential 
applied to the grid. The electron temperature and 
concentration for exposure f were of the same 
order of magnitude as in e above. The mercury 
pressure for this group of exposures was 0.08 mm 
of mercury. 

Exposures ¢ and j were taken with a slightly 
higher positive grid voltage than in f, g and h. 
Unfortunately, measurements were not made on 
the probe volt-ampere characteristics. In this 
case the visible glow around the grid becomes 
more intense and the volt-ampere characteristics 
of the accelerating grid show an abrupt increase 
in electron current. The other conditions were the 
same as in /, g and A. 

Photographs of the emitted radiation at the 
first window (upstream) indicated no perceptible 
change in intensity on the addition of hydrogen 
under the same conditions as obtained in Fig. 1 
and Table I. Apparently the small gas pressures 
used were not sufficient to alter conditions in 
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the arc region and the amount of hydrogen diffus- 
ing against the vapor stream was small. 


DiscussION 


One would expect on simple theories* of the 
recombination of positive ions and electrons that 
the intensity of the spectral lines radiated as a 
result of such processes would vary as the square 
of the electron concentration (assuming positive 
ion and electron concentrations to be equal) 
and to decrease with increase in electron tempera- 
ture. The spectra a, b, c and d are then not in- 
consistent with what one might expect on the 
assumption that we are here only concerned with 
recombination. On the other hand the marked 
enhancement of all of the arc lines in exposure ¢ 
for the highest grid voltage cannot be associated 
with recombination. In the latter case the elec- 
tron concentration is practically equal to that 
in exposure d. The small difference in electron 
concentration in exposure e¢ is offset by an in- 
crease in electron temperature. Apparently, as 
will be made more obvious later, excited atoms 
are produced in the vicinity of the grid. This 
additional excitation upon being carried down the 
stream either in the form of excited atoms or 
resonance radiation makes a contribution to the 
emitted radiation at the point of observation. 
Under such conditions the intensity of the 43025 
group which has its origin in the 7*Dyws and 
7 states is noticeably greater than 2967, 
the upper state of which is the 6*D, state. As 
will be discussed later in connection with the 
quenching by hydrogen, we are probably here 
concerned with a form of excitation in which the 
*Po.» states plays an important part. If such is 
the case, the decrease in intensity of the radia- 
tion in a, b, c and d cannot definitely be assigned 
to a decrease in recombination alone. Changes in 
electron temperature may decrease the lifetime 
of atoms in the metastable 6 *P, and in this way 
alter the importance of any process in which 
they are involved. It would probably be ex- 
tremely difficult to examine experimentally in a 
quantitative way such changes in populations of 
excited states. 

The data of Table I indicate that for small 
additions of hydrogen the line radiation of the 


*F. L. Mohler, Rev. Mod. Phys. 1, 216 (1929). 
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afterglow under these conditions decreases in 


_intensity. For the higher hydrogen pressures 


with the same grid conditions, the density of the 
lines with the exception of \2537 increases until 
at the highest pressure used (approximately 10° 
mm Hg) their densities exceed their initial 
values. The small but definite decrease in the 
intensity for small additions of hydrogen may be 
associated with the destruction of excited atoms 
by collisions of the second kind. The recovery of 
the intensity of the afterglow with grid floating 
and with higher hydrogen pressures may possibly 
be due to an increase in the probability of re- 
combination due either to a lowering of the 
electron temperature or possibly to the fact that 
such hydrogen molecules or atoms may provide 
third bodies which aid in recombination. 
Whether or not the addition of hydrogen causes 
an increase or decrease in the intensity of the 
afterglow should depend on the relative impor- 
tance of the contribution associated with excited 
atoms and the radiation due to recombination. 

As is shown in Fig. 1 f, g and h, the quenching 
effect of hydrogen is much more pronounced with 
higher grid voltages. The addition of hydrogen 
to a measured pressure as low as 4X 10~* mm of 
mercury shows a very noticeable decrease. With 
slightly less positive voltages on the grid (about 
3.5 volts), at these mercury vapor pressures, 
it is possible practically to wipe out all lines 
except 42537 which is still detectable but very 
much weaker. The decreased effectiveness of 
hydrogen in the quenching of some of the lines 
in the afterglow at the higher grid voltage (ex- 
posures i and j) may be due to the fact that the 
grid may become surrounded by a sheath and 
thus lose its ability to affect the diffusion of 
positive ions and electrons from the arc region 
to the point of observation. 

The fact that the half value pressure for the 
quenching of resonance radiation by hydrogen® 
is of the order of 0.2 mm of mercury, a value 
considerably greater than those effective here, 
would seem to indicate that the quenching of the 
metastable 6*P,) and 6*P: states may also be 
involved. In connection with the marked en- 
hancement of the \3025 group with no hydrogen 


* Ruark and Urey, Atoms, Molecules and Quanta, p. 497. 


present as in exposures e and i of Fig. 1 it is of 
interest to point out that the energies of these 
states is approximately the sum of the energies 
of the 6*P, and 6*P, states. It may therefore 
well be that these as well as other higher excited 
states are filled by collisions between two such 
excited atoms as might be suggested by the 
work of Rayleigh.'® It is also possible that these 
levels may be excited by collision between such 
excited atoms and some of the high speed elec- 
trons present in the discharge under these condi- 
tions. Under both conditions one would expect 
the emitted radiation to be quenched by hydro- 
gen. The introduction of a stream of hydrogen 
such as to build up a pressure of only 8x 10-* 
mm of mercury at the point of observation 
quenches the 43025 group as well as other lines 
completely as shown in exposure j. The line \2967 
although weaker under such conditions is defi- 
nitely present on the original spectrograms. 

As indicated above, this quenching action of 
hydrogen on the afterglow is apparent under the 
conditions of this experiment not only at high 
electron speeds but under conditions where the 
electron temperature is quite low. It is, of course, 
much more apparent under the former conditions. 

It would seem therefore that the intensity of 
the lines in the afterglow of mercury vapor, 
when excited under conditions such as in this 
experiment (arc voltage, Hg vapor pressure, 
electron temperature etc.) is in part dependent 
upon the concentration of atoms in the 6 *Po. 
states. 

The results here presented are not necessarily 
inconsistent with the assumption of Webb and 
Sinclair that they were primarily concerned in 
their work with recombination of ions and elec- 
trons, and, that excited atoms played a minor 
role in their work. Their measured arc voltage 
was about 8 volts and their arc current was ap- 
proximately 10 amperes. It is, however, the 
purpose of this note to present additional evi- 
dence that the intensity in the mercury after- 
glow of some of the lines which have their origin 
in the higher energy states may be perceptibly 
altered under certain conditions by changes in 
the concentration of atoms in the *Po1. states. 


‘Lord Rayleigh, Proc. Roy. Soc. A132, 650 (1934); 
A135, 617 (1932). 
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On the Theory of Fusion 


K. F. HerzFe_tD AND Maria Goeprert Mayer, The Johns Hopkins University 
(Received June 18, 1934) 


Although equilibrium between two phases is in general 
determined by both phases, fusion stands somewhat apart 
in that a crystal cannot be overheated. As a type of the 
simplest possible crystal, the exact equation of state of 
frozen argon has been developed here. It has been found 
that the pressure as function of the volume at a given 
temperature has a minimum, which means the breakdown 


of the crystal under these conditions. The temperature at 
which this minimum occurs at zero pressure is interpreted 
as the melting point. It has been found in accordance with 
experiment that a positive pressure is necessary for the 
existence of crystalline helium but the numerical agreement 
is bad. 


INTRODUCTION 


OME time ago, we attempted! to get a better 

understanding of the liquid state through 
investigation of a simplified one-dimensional 
model but succeeded only in recognizing the 
conditions for the existence of definite phases. 
Here the same problem is attacked in a different 
manner. The equation of state of simple crystals 
(e.g., frozen rare gases) is calculated from. the 
forces according to lattice theory. We find that 
the isotherms in which p is plotted against V 
have a minimum. Now a state in which the 
pressure increases with volume is unstable. 
Accordingly the crystal breaks down at the 
minimum pressure. Above the temperature at 
which this pressure minimum lies at zero pres- 
sure, the crystal cannot exist without external 
pressure. This temperature we interpret as the 
melting point. 

In general, equilibrium between two phases is 
not defined by any particular happening in one 
of the two but by the equality of the free 
energies of the two at given pressure and 
temperature. For example, the equilibrium 
between liquid and vapor being given by the 
equality of their free energies, one can undercool 
the vapor or overheat the liquid without any 
discontinuous change in the properties of either. 
The same is even true for equilibria between 
allotropic crystals. It is possible to overheat 
grey tin or undercool white tin without any 
discontinuity appearing anywhere. The only 
known exception is the fact that no crystal can 
be heated above the melting point at the given 


1K. F. Herzfeld and M. Goeppert Mayer, J. Chem, 


Phys. 2, 38 (1934). 


pressure. It seems therefore reasonable to assume 
that fusion is conditioned by something happen- 
ing in the crystal alone and to identify the point 
of instability of the crystal, which was described 
above, with the melting point. But then another 
question arises. At the melting point, the melt 
is-in equilibrium with the crystal, ie. the 
pressure, temperature and free energy of the 
melt must have values which are prescribed 
completely by the crystal (while in other equi- 
libria, as stated above, the pressure is selected so 
as to give the equality of free energies). The 
only variables that seem to be left to the liquid 
to adjust its free energy are the values of the 
constants that appear in the potential energy 
and frequency and these are again determined 
by the arrangement of the molecules. Accord- 
ingly, the arrangement of the molecules in the 
liquid has to adjust itself so as to give the 
correct free energy.'* What this arrangement has 
to be is not investigated here but is reserved for 
later consideration. 

If this view of fusion is correct, one would 
expect, for all normal substances, the volume 
of the liquid to be greater than that of the solid 
for the same temperature. It is then clear why, 
in the process of melting, the crystal does not 
melt in different small regions inside. A small 
molten region inside would have to contain the 
liquid in a highly compressed state, which would 
have a considerably higher free energy than the 
surrounding solid. It is, however, quite possible 
for the fusion to occur in channels communi- 


cating with the outside. 


1s However it seems that this arrangement should be 
the one to give the smallest free energy of the liquid. 
How that is possible is not yet clear. . 
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On the other hand, it seems not impossible for 
small crystals to persist in a liquid, a fact for 
which there appear to be experimental indi- 
cations.? Up to now three other theories of fusion 
exist. 

Lindemann‘ made the assumption that melting 
takes place if the amplitude of the vibrations is 
a certain fraction of the distance between 
nearest neighbors. Upon determining this frac- 
tion from the data for one crystal of each type 
and applying it to all crystals of the same type, 
one gets very good agreement. 

Braunbeck* considers the mutual vibrations 
of two parts of the lattice, taken as rigid, 
assumes a suitable form of the mutual potential 
curve and finds positions of minimum energy. 
One of these he interprets as belonging to the 
solid, the other to the liquid state. The melting 
point is then determined as the temperature at 
which the free energies corresponding to these 
two states are equal. 

Raschevsky® has come close to the ideas 
followed in this paper by pointing out that the 
heat of fusion is of the same order of magnitude 
as the elastic energy necessary to expand the 
lattice to the breaking point, i.e., the point 
where the tension has a maximum. 


THE ENERGY AND THE ELASTIC CONSTANTS FOR 
A FACE-CENTERED LATTICE 


Argon has a face centered lattice.* Call r the 
distance to the nearest neighbor, which is at the 
center of the face, so that the side of the cube is 
ry 2. The volume per atom is 


(ry 2)?/4=r°/¥ 2. (1) 


As potential energy ¢; of attraction between two 
atoms we use a van der Waals expression (A, 


a constant) 
(2) 


which, summed up over the lattice, gives’ 


*P. Othmer, Zeits. f. Allgem. und anorg. Chemie 91, 209 
(1915). A. Goetz, R. C. Hergenrother and A. B. Focke, 
Phys. Rev. 34, 546 (1929); A. Goetz, Phys. Rev. 35, 193 
(1930); A. Soroos, Phys. Rev. 41, 516 (1932). 

°F. A, Lindemann, Phys. Zeits. 11, 609 (1910). 

*W. Braunbeck, Zeits. f. Physik 38, 549 (1926). 

5 N. Raschevsky, Zeits. f. Physik 40, 214 (1927). 
om and Cl. v. Simson, Zeits. f. Physik 25, 160 

7E. Jones (Lennard-Jones) and A. E. Ingham, Proc. 
Roy. Soc. A107, 636 (1925). 


= —14.45394 (3) 
with The repulsive potential is written® 
(B, p constants) which gives for the whole 
crystal 
(3’) 


This includes the atoms at the corners of the 
cube and in the centers of all the faces on all 
four cubes surrounding the atom considered here. 

For p two values were used: p = 0.345A, which 
Born and Mayer found to fit best the behavior 
of alkali-halides, and p=0.2091A which Bleick 
found theoretically for the interaction of two 
Ne atoms. Unfortunately, it is not possible to 
determine A and B directly from the normal 
valve of r, ro>=3.83A and the normal heat of 
sublimation because the zero point energy plays 
a considerable part in the heat content and 
displaces the lattice distance by about 0.05A 
beyond the value given by the minimum of 
the potential energy alone. 

A method of trial and error leads to® (N, 
Avogadro number) 


NA =6.706 X10-" erg cm®, B=1.27010" 
for p=0.345. 


NA =4.729X10-" erg cm®, B=0.9484x 
for p=0.2091. 


These give as the heat of sublimation,’® the 
value 1655 cal. The experimental values are 
rather uncertain. 

There are three elastic constants, given by" 


5M. Born and J. E. Mayer, Zeits. f. Physik 75, 1 (1932); 
W. Bleick, Dissertation, The Johns Hopkins University 
(1933). W. E. Bleick and T. E. Mayer, J. Chem. Phys. 2, 
252 (1934). , 

®* The value of NAB=8.517X10", which we find for 
p=0.345, is to be compared with NX 10~%e~2.660/e = 13.53 
X10" which would follow from J. E. Mayer and M, 
Huggins, J. Chem. Phys. 1, 633 (1933), if one puts the 
radius of 2Av=K*+Cl-. Furthermore, the value of NA 
is to be compared with the theoretical values of 3.3 and 
3.7X10-* calculated by F. London, Zeits. f. physik. 
Chemie B11, 222 (1930). 

1A. Eucken, Verh. D. ao Ges. 18, 4 (1916); A. 
Eucken and F. Hauck, Zeits. f. physik. Chemie 134, 161 
(1928); F. Born, Ann. d. Physik 69, 473 (1922). 

1M. Born, Dynamik der Krystallgitier, p. 536, 548, 554, 
Leipzig, 1923; M. Born and M. Goeppert Mayer, Handb. 
d. Phys. 24, 2nd half, 2nd Ed., p. 630, Berlin, 1933. en 
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dV 2 rdrXrdr 
rdr\rdr 
d 
Cag =— — (4”) 


where V is the mol volume. 

At the minimum of the potential energy 
d@/dV=0 so that there are only two elastic 
constants, because ¢2=Cyu. This does not hold, 
however, for other distances, where the general 
form must be used. 

By using (2) and (2’) one finds” 


2r+wv2p 
2-Drlp 
r+p 


NA, r+rp 
i+ 
p 
3p 
+2 
r+p 


242.544 d 
=, 


2 


e 


Brt+rp 
r+p 


r® 


measures the tension in one direction due to an expansion in 
the same direction with the cross dimensions kept constant, 
C12 measures the tension in the cross direction which tries to 
contract the specimen and determines largely the Poisson 
constant, ¢4, measures the shearing stress. The summation 
is to be extended over all the lattice points with the 
coordinates x, y, 2. 

122A calculation of the numerical factors of r* was 


effected by directly summing up over the nearest 168 
atoms, i.e. up to and including a distance 3r. Then an’ 


integration was performed over a continuum reachin 
from a sphere of radius 3r to infinity and from this half 
the contribution of the atoms at distance 3r, which had 
already been taken into account in the summation, was 
subtracted. The values calculated in that way were 2.541 
and 1.135. But in a cubic crystal one must have 


(x*+-2x*y*) = /r = x27? = x2 = 


As 2.5414+2X1.135=4 of 14.433 instead of 14.454 (3), 
the constants were raised to the values above. 


THE FREQUENCY 


According to Debye the maximum frequency 
of a cubic crystal is 


v=w(3N/4rV)! (7) 


where w is the sound velocity. In dependence of 
the direction this is given by a determinant of 
the third order. The method used by Born and 
K4rman"™ which consists in a development into 
powers of 3¢2 — Cr, cannot be applied 
here because this quantity is by no means small. 

But on the other hand there is no particular 
reason why one should average just in the way 
proposed by Debye, -using in (7) the quantity 
found from 


1/w* = (1/wy*)+(1/we*) +(1/ws?), 


where w, %2, ws mean the three roots of the 
characteristic equation and the bars an averaging 
over the direction. It is much simpler to form 
one of the symmetric functions of the roots and 
average. We have used 


= D*w,*w,*w;? = 
(C11 — Cas)* — 4012" Jos? (xr) (yr) 


+cos* (yr) cos* (zr) +cos* (zr) cos* xr] 
+0 (611 — — 12 (11 — Cag) + 
Xcos* xr cos* yr cos* sr (8) 


in which D is the density. In terms of numerical 
values and the atomic weight M this becomes 
+ (C11 — Caa)* — 4012" (1/105) 
X — Caa)® — 12( — Cag) +1612" ]}. (9) 


THE EQUATION OF STATE 


The total energy at 0° abs. and any volume 
is given by the lattice energy plus the zero point 
energy } Nh» (4h» for each degree of freedom). 


N® 3 
(10) 


The limiting » calculated in (9) is for the Debye 
spectrum 4/3 of the average ¥ meant above. At 

13 Born, ref. 11, p. 647, Born and Goeppert-M 4 
p. 672; M. Born oa Th. von Karman, Ph Zeits. 14, 


15 (1913); L. Hopf and G. Lechner, Verh. d. D. phys. Ges. 
16, 643 (1914). 
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Taste I. Frequency, elastic constants and energy. r in A, energy in kilojoule, p =0.345A, © =hv/k in degree kelvin. 
t , 3.73 3.78 3.83 3.88 3.93 3.98 4.03 4.08 
ts) 94.31 81.90 72.80 60.31 $4.47 48.98 43.68 
X 107 3.901 2.742 2.015 1.573 1.207 0.9064 0.6606 0.4583 
12 X 107 2.273 1.613 1.183 0.9198 0.7028 0.5249 0.3800 0.2609 
y cu X 107 2.115 1.592 1.262 1.048 0.8657 0.7122 0.5830 0.4737 
—7.627 —7.647 —7.600 -7.501 —7.356 —7.129 —6.974 —6.756 
; 3 Nh + 0.882 0.766 0.681 0.623 0.564 0.510 0.458 0.409 
d -U 6.745 6.881 6.919 6.878 6.792 6.619 6.516 6.347 
7 
q Taste II. Pressure. (Fig. 1.) V in cc, pressure in 10° dynes /cm*, p =0.345A 
q r 3.73 3.78 3.83 3.88 3.93 3.98 4.03 4.08 
q V 22.26 23.17 24.10 25.06 26.04 27.05 28.08 29.14 
q else 0.167 0.143 0.109 0.097 0.102 0.102 0.104 0.113 
av 
| pa 15.82 2.08 —7.93 —12.82 — 16.29 —18.73 —20.30 —21.28 
15 10.9 7.42 6.04 5.74 5.18 4.78 4.62 
31 13.0 —0.51 —6.78 —10.55 —13.55 —15.52 — 16.66 
70° ps 17 15.6 12.61 11.65 12.67 13.10 13.88 15.50 
? 48 28.6 12.10 4.87 2.12 —0.45 —1.64 —1.16 
76° ps 19 i7.6 14.16 13.04 14.13 14.57 i541 17.17 
p 50 30.7 3.65 6.26 3.58 1.02 —0.10 0.51 
100° ps 29 25.8 20.46 18.67 20.07 20.55 21.56 23.84 
| ft 60 38.8 19.95 11.89 9.52 7.00 6.04 7.18 
v/dV). 


higher temperatures, the free energy pro mole 
at constant volume is 


f=U+F, (11) 
where 


rc aT 
F=f | —E. (11’) 
0 0 T 0 T? 


Here C is the specific heat and E= (\"CdT. 
The pressure is given by 


Here it has been assumed that F/T is a function 
of v/T only. F and E can be found tabulated." 

It is necessary for possible existence of a 
state that 


(ap aV) <0. (13) 
If a p, V curve for constant T shows a minimum, 


the crystal breaks down (cannot exist) beyond 
that minimum. 


p=—of/aV (12) 
| and is accordingly made up of three parts: Tue Discussion OF THE RESULTS 
P= PitPrt Ps. (12") All total pressures are underlined, the mini- 


p: is the elastic tension 


mum pressure for each temperature is doubly 
underlined. The results are shown in the Tables 


—(d/a V)(N®/2). (12’’) I-IV 
p2 is the zero point pressure If, according to this interpretation, the melting 
a3 a33 point is given by 
p=0, ap/aV=0, a%p/aV?>0, (14) 


Both these values depend only on V, not on 7. 
At higher temperature, one has to add the 
thermal pressure'* 


“F/T = — fy" (E/T?)dT = — If we 
T/v=x and assume that C is a function of x only 


t 
ke = fCdT=vfCdx. F/T=— fx*dx Cdx is a function 
of x only, say ¥(x), (dy /ax) = —(Ev/T*). Then 


(8F/8V) = /ax)(dx/aV) 
= = (E/»)(dv/AV). 


it follows that the compressibility of the crystal 
at the melting point must be infinite. The same 
should be true for the coefficient of expansion. 
In considering the tables, one sees that for 
the equilibrium distance r=3.83A ©=73 for 
p=0.345 and @=96 for p=0.2091 while the 
experimental value is 0 = 85. For p=0.2091, O is 


“ F. Simon, Handb. d. Physik 10, 368 (1926), Berlin. 
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TABLE III. Frequency, elastic constants, energy. r in A, energy in kilojoule, p =0.2091 A. 
re r 3.73 3.78 3.83 3.88 3.93 3.98 4.03 
128.5 111.2 96.46 82.81 70.56 59.06 48.44 
‘583 eu X10-" 7.165 5.169 3.660 2.523 1.673 1.0383 0.5724 
2609 cis X 107% 3.908 2.796 1.987 1.328 0.8594 0.2578 
1737 cu X 107% 3.761 2.821 2.099 1.547 1.1263 0.8087 0.5634 
xe —7.891 -7.923 —7.841 —1.666 —7.425 —7.143 —6.830 
wii ; Nh : > 1.202 1.040 0.902 0.775 0.660 0.552 0.453 
47 -U 6.689 6.883 6.939 6.891 6.765 6.591 6.377 
4 
== Taste IV. Pressure (in 10° dynes/cm*), volume in cc, p =0.2091 A. (Fig. 2.) 
08 
14 r 3.73 3.78 3.83 3.88 3.93 3.98 4.03 
113 "4 22.26 23.17 24.10 25.06 26.04 27.05 28.08 
0.168 0.156 0.156 0.160 0.169 0.185 0.198 
62 | 14.72 —2.45 —14.17 21.91 —26.69 —29.39 — 30.56 
06 pe 20 16.3 14.09 12.39 11.17 10.20 8.98 
350 0° 35 13.8 —0.08 —9.52 —15.52 —19.19 —21.58 
-16 50° ps 78 S84 10.20 11.98 1687 
42 21.6 8.76 0.68 —3.54 —4.80 —4.71 
61° ps i13 12.56 14.13 16.26 19-79 22.10 
84 ? 46 25.1 12.48 4.61 0.74 0 0.52 
65° pa 12 127 1394 15.59 17.84 20.96 24.02 
sat ? 47 26.5 13.86 6.07 2.32 1.77 2.44 
100° ps 25 7.1 26.64 28.85 32.05 36.64 41.00 
58.9 26.56 19.33 16.53 17.45 19.42 
ion 
60 
50 
13) 
im, 40 
p= 02091 A 
N =O.345 A \ \ 
™ 
ni- > \ XY 
bly _\ WANN 
v 0 76% 61C, 
ing 70C 
aid 50 
14) 4 
Z0 ie 
tal 
me 
S75 376 385 386 393 316 403 408 573 376 383 388 393 3% 403 
ad rinA rmA 
the Fic. 1. Fic, 2. 
9 is Fic. 1. pr isotherms for argon. p =0.345A, the value Born and Mayer found for the alkali-halides. 
‘ Fic. 2. pr isotherms for argon. p=0.2091A, the value found by Bleick and Mayer for He atams. 
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83 for r=3.88, which corresponds to 50°. On 
theoretical grounds, one should expect to have 
p rather alike for argon and neon and therefore 
prefer p=0.2091. 

The melting points are both too low, 76° and 
61° compared with 84° observed. It is true that 
we have, following Born, used only a first 
approximation insofar as we considered the 
vibrations as purely harmonic around the given 
lattice distance with the frequency appropriate 
for that distance. A higher approximation 
would need rather lengthy calculations, but one 
should expect that the effect would be to lower 
the melting point still further. 

The thermal pressure is determined essentially 
by @lnv/dV. A higher melting point would 
result if d In v/8V would be lower ; then a higher 
temperature would be necessary for the same 
thermal pressure. A smaller means a 
slower falling off of the repulsive forces. 

In comparing the melting points and 0's (or 
v's) for the two p-values, one might be astonished 
to find the higher 7 with the lower »(p=0.345) 
in contradiction to the Lindemann formula. 
But we are comparing here for constant heat of 
sublimation, while the variation of » in the 
Lindemann formula measures essentially the 
variation of this heat. From an inspection of the 
table, 9 In »/dV is seen to vary little with V, in 
agreement with Griineisen.'® 


HELIUM 


The properties of condensed helium have 
recently been the subject of much experimenta- 
"8 It is found that crystallized helium 
can exist only under pressure (25 atm. at 0°). 
Simon" has already interpreted this as due to a 
zero point pressure > | 

He has given as provisional data: 

Heat of sublimation: 11 cal.=46.1 joule. 


The density of solid He at 4° and under a 
pressure that provides equilibrium (~ 100 atm. ?) 
is found to be 0.23 (V =17.4 cc), the compressi- 
bility at 3.7° and 115 atm. is about 1.5x10~-° 
abs. units and the characteristic temperature 
=32.5. 


* E. Griineisen, Ann. d. Physik 39, 257 (1912). 


17 W. Keesom and co-workers, Leid. Comm. No. 184b, 
190b, 216b, 219e, 221e. Suppl. 67b, 7le, 76b. Physica 1, 
128, 161 (1934). 

18 F, Simon, Nature 133, March 24 and April 7, 1934. 


Now the potential energy of two helium 
atoms has been calculated from the model by 
Slater and Kirkwood" to 


where ao is the Bohr radius.*° With this, the 
same calculation as for argon was performed 
under the assumption that solid helium, like 
argon, has a face centered cubic lattice. The 
result is given in the Tables V and VI. 

One sees that there is a pressure minimum at 
0° at about 3.12A. The pressure there is positive, 
so that the crystalline helium cannot exist with- 
out pressure. However, the volume is much 
smaller than the experimental one (13 cc against 
17.4; at V=17.4 and r=3.45 one is already far 
beyond the maximum of ),, all the elastic 
constants are negative). The minimum pressure 
is too high, about 1250 atm. instead of 25 atm., 
finally © is too large (about 55° instead of 32°). 
On the other hand, the heat of sublimation 
seems about right. 

There are two possible explanations for these 
deviations: The repulsion is considerably larger 
than given by (15), or the assumed lattice is 
wrong. Now the only other lattice that might 
be plausible is the body centered cubic one. 
While the face-centered lattice has a mol- 
volume Nr*/ 2, the body centered has one of 
N(4/3v3)r'*, which, for the same r is larger by the 
factor 1.088. On the other hand, the summation 
of ¢ over the lattice leads to a decrease of the 
numerical factor for the repulsive energy from 
12/2 to 8/2 while for the attraction this sum- 
mation results only in a decrease of the factor 
from 4X 14.454 to? }X<12.253. That means that 
for the same ¢ the minimum of the potential 
energy and the place where — ~; has a maximum 
are shifted to smaller r. Therefore, even with 
this lattice, this maximum of — ,; would occur 
around V=13.5 cc. 

Accordingly, an attempt has been made to 
repeat the calculation with a changed potential 
energy for the face centered lattice. The repulsion 
was increased by a factor 3. This moves the 
minimum pressure out to the distance where it 


(1988) C. Slater and J. G. Kirkwood, Phys. Rev. 37, 682 
lt might be remarked that this amounts to a p=0.219, 
very close to the Mayer-Bleick value for neon. 
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Tasie V. Energy, elastic constants and frequency for He. Energy in joule, r in A. 


r 2.8 3.0 3.05 3.10 3.15 3.20 3.25 

8 169.5 80.8 68.75 57.29 46.41 34.35 

eu X10°% 164.6 26.47 16.99 9.856 4.648 0.901 — 1.666 

cu X10* 94.86 14.77 9.228 5.151 2.204 0.089 —1.344 

cou X10-* 80.36 18.29 13.35 9.577 6.715 4.540 2.940 

“ — 1586 —756.4 — 645.5 —536.1 434.3 —321. 

Nh 598.0 593.5 570.8 544.7 516.5 487.4 

-U +988 +162.9 +747 —8.6 —82.2 — 166.0 -U 
TABLE VI. Pressure. Pressure in 10° dyne/cm?, V in ce. U=Uo+(U"/2)(V— Vo)*. 

r 3.0 3.05 3.10 3.15 3.20 cannot 

4 Ford In »/8 V two alternate approximations were 

one 0.379 0.322 0.310 0.333 0677 made, either dln »/8V was assumed constant or 

—3.52 —4.12 —4.43 —4.51 —448 was put proportional to For the first 

22.5 18.4 16.9 17.2 —33 

0° p 19.0 14.3 12.5 12.7 285 case the calculations were very simple. If the 

4° Ds 0.021 0.041 0.21 


should be (actually a little further), namely to 
3.50A (V =18.4 cc). That reduces © to 50, the 
minimum pressure is reduced to about 530 atm. 
but it gives a positive energy content, that is to 
say a negative heat of sublimation. The zero 
point energy is then about 450 joules, the lattice 
energy about 250, having been cut to about 4 
of its former value by the increase of the distance 
and the repulsive term. The compressibility also 
comes out far too small. Without further calcu- 
lations it cannot be said what has to be done to 
improve the agreement. We reserve this calcula- 
tion and shall only indicate that a slower falling 
off of the repulsion might be necessary. 


SUBSTANCES WITH HIGH MELTING POoINTs 


Approximations for substances with melting 
points sufficiently high that the classical specific 
heat is reached should simplify the calculations 
considerably because in that case the zero point 
energy can be neglected and E in equation (12’’’) 
can be replaced by 3RT. However, we do not 
know the potential energy in the metals well 
enough to perform the calculation and in the 
case of ionic crystals there is still a certain 
difficulty in calculating the frequencies cor- 
rectly.” Attempts at an approximate calculation 
were made. The lattice energy U was developed 
into a power series up to the 3rd power 


*1 See, for example, Born and Goeppert-Mayer, Handb. 
d. Physik 24, 2nd half, 733, 2nd edition. 


experimentally found coefficient of expansion is 
given by 
a=ao(1+ fT). 
then the melting point should be given by 
T =1/28. 


For the second assumption the calculations are 
rather complicated but in both cases the calcu- 
lated melting points are far too low. 

The general results of this paper, however, 
should hold for all substances. The existence of 
an upper limit of temperature for the stability 
of a crystal is based on the appearance of a 
minimum in the p, V isotherm, shifting to more 
positive pressures with increased temperature. 
Now such a minimum should appear in the p, V 
curve of every crystal, becavse with increasing 
volume the elastic tension must finally go to 0. 
On the other hand, the thermal pressure should 
increase to infinity. It is namely given by (12’’’’). 
v goes to 0 at a certain volume, probably with 
the first power of the deviation from that 
volume. Accordingly @lnv/@V must go to 
positive infinity when this volume is approached 
and the combination of the negative elastic 
tension, decreasing with volume, and the positive 
thermal pressure, increasing to infinity with 
increasing volume, must give a minimum. Be- 
cause of the factor E in the thermal pressure 
this minimum will always shift to more positive 
pressures with increasing temperature. But, of 
course, for a numerical calculation the exact 
energy formula must be known. 
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The energy of interaction between free electrons in an 
electron gas is considered. The interaction energy of 
electrons with parallel spin is known to be that of the 
space charges plus the exchange integrals, and these terms 
modify the shape of the wave functions but slightly. The 
interaction of the electrons with antiparallel spin, contains, 
in addition to the interaction of uniformly distributed 
space charges, another term. This term is due to the 


fact that the electrons repell each other and try to keep 
as far apart as possible. The total energy of the system 
will be decreased through the corresponding modification 
of the wave function. In the present paper it is attempted 
to calculate this ‘correlation energy’’ by an approximation 
method which is, essentially, a development of the energy 
by means of the Rayleigh-Schrédinger perturbation theory 
in a power series of e°. 


1. 


HE attempt has been made in previous 
work! to give a more general expression 
for the wave function of free electrons in metals 
than that provided by Hartree’s method of the 
self-consistent. field?: * or Fock’s equations. The 
form of the wave function assumed in Fock's 
equations for a system of 2n electrons, occupying 
n doubly-degenerate states is 


n 


where x stands for three Cartesian coordinates 
of electrons with upward spin, and y for those 
of electrons with downward spin. The y, are 
the solutions of a Schrédinger equation in which 
the potential of the charge distribution of the 
other electrons enters as well as the potential 
arising from the ions. 

In a metal the charge distribution of all 
electrons is practically unaltered by removing 
one so that the second quantity may be replaced 
by the former and the potential for a given 


electron at the point u is given by adding to the 
Coulomb field of the ions the fields of all electrons 
with parallel and with antiparallel spin. The 
former distribution may be obtained by inserting 
u for x, in (1) and integrating over all coordinates 
except x, and u, while the latter is obtained by 
a similar operation with the exception that the 
integration should be carried out over all 
coordinates except y: and u. 

Actually, it had been shown in':‘ that the 
wave functions y, of the free electrons in a 
Na-lattice are very nearly plane waves ¢**:*/% 
where L is the cube edge of the crystal and » 
stands for a set of three integers, v-x denotes 
the scalar product of v and x. Hence the charge 
distribution of the electrons with opposite spin 
is practically uniform, that of the electrons with 
parallel spin uniform with a “‘hole’’ around u.° 

In no wave function of the type (1) is there a 
statistical correlation between the positions of 
electrons with antiparallel spin. The purpose of 
the aforementioned generalization of (1) is to 
allow for such correlations. This will lead to an 
improvement of the wave function and, therefore, 
to a lowering of the energy value. This energy 
gain will be called “correlation energy.” 


2. 
The new form of the wave function, assumed in' was 
Wily Mi) Xn) ¥i(y1) ¥i(yn) 


1 E. Wigner and F. Seitz, Phys. Rev. 46, 509 (1934). 
2D. R. Hartree, Proc. Camb. Phil. Soc. 24, 89 (1928). 


4]. C. Slater, Phys. Rev. 45, 794 (1934); A. Sommerfeld 
and H. Bethe, Geiger-Scheel’s Handbuch der Physik, 


*J. C. Slater, Phys. Rev. 35, 210, 1930; V. Fock, Vol. 24, 2nd part, 2nd edition, p. 406. 
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Zeits. f. Physik 61, 126 (1930). 


5 E. Wigner and F. Seitz, Phys. Rev. 43, 804 (1933). 
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of the symmetric group). If, however, the functions 1s in this case near to (2). 
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which contains the functions y,(y:- - -y,; x) which are different functions of x for different configura- 
tions y:, ---, Ys Of the electrons with opposite spin, instead of the y,(x). It is proposed to find 
the best wave function, i.e., that with the lowest total energy of the form® (2). The y, «--, ¥ in 
¥,(¥1"-*¥n; X) are to be viewed merely'as parameters. (Cf. Eq. (23) ref. 1.) The total energy of the 
wave function (2) was previously calculated and relative to the solution of Fock’s equations yielded 
the following energy 


n! 2m « 


where |y| denotes the second determinant of (1), the summation « runs over all 3m coordinates y 
and that over » and »’ over all occupied states, i.e., over all indices occurring in the wave functions 
in (1) or (2); ¥,(y; x) stands for ¥,(y:---y,; x) and dy for dy;- - -dy,. The quantities ¢, are the integrals 


= x)*{ (h? /2m) (Az, + Ay, + + A,,) } x, )dx, 
+ f (3a) 


where again V(y,---y,; x) is the difference between the potentials at the point x of a charge distribu- 
tion corresponding to ¥:(y), Waly), on the one hand, and point charges at y, Ya, 
on the other. 

It is necessary now to assume for y,(y:, ---, Ya; x) the form 


and that if v and »’ are both occupied states 
x1) /dy.)dx,=0 (4a) 


so that the second term in (3) vanishes. Both (4) and (4a) will turn out to be correct in the approxi- 
mation to be used. By means of (4) it is possible to transform (3a) so as to get rid of all derivatives 
with respect to the y, after which one may minimize (3) by minimizing ¢,(,, --:, ya) for every 
combination of the y. This would lead to a differential equation for the ¥,(y:---y,; x) in which the 
y would be merely parameters. (The solutions ¥,(x) of Fock’s equations which also enter into this 
equation are supposed to be known.) The result of the transformation is especially simple if one uses 
for ¥,(x) namely, 


a= x)*{ V—(h?/m)(A.—(2riv/L) grad z)} x)dx. (3b) 


In addition to the energy contribution (3) which is negative and was calculated in' there is a 
further one which is generally positive. This arises from the fact that the probabilities of the relative 
distances of electrons with upward spin are changed by the transition from the y,(x) to the 
¥.(¥1, --+, Yn; x). Since the latter will be large for x’s, which lie in regions comparatively free from 
y's, the distribution of the x's will not be uniform throughout space and they will be nearer together 
than they were under the previous assumption. 

* The form (2) of the wave function is certainly not the y,(yi---¥.; x) are not too different from the functions 
correct one. It does not belong even to one single multi- _¥,(x), they all belong to very low multiplicities. This is the 


plicity but is a linear combination of functions of different only case, anyway, in which the present approximation is 
“multiplicities” (belonging to different representations of good and it can be expected that the real wave function 
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In order to evaluate this energy change, one may first calculate the probability of two electrons 
with parallel spin being at the points x; and xz, respectively, if the complete wave function is given 
by (1). Under the approximate assumption that y,(x) =yo(x)e***"*/¥ one obtains for this in the way 
given in reference 5, 


| ¥o(x2) (5) 
j in (r/d) —(r/d) cos (r/d)\* 
4 sin cos (7 
q =1-9 5 
| is the probability of the distance |x,—2:| =r for free electrons with parallel spin, 

d= (v/3x*)'= (4/9x)'r,= 0.5211, (6) 


a is 1/(2r) times the wave-length of the fastest electron, v is the atomic volume and r, the radius of the 
a sphere with this volume. If we make the assumption y,(y, x) =y,(y, x)e?**“-”)-=/“, an expression of 

| the form (5) is valid for the wave function (2) as well as for (1) and the change of mutual potential 
energy of the electrons with upward spin arising from the transition from (1) to (2) is 


q where 

= f dx ydxo(me*/2r) {| x1) |? — | |*} g(r). (7a) 
f The } enters once again because the interaction of a pair of electrons should be counted once only. 
| In order to be able to evaluate the integral (7a), g(r) has been replaced by 

1 +-1.6r/d+1.2(r/d)?), (8) 


| which, as is shown in Fig. 1, runs rather near to g(r). 


i 3. The situation is somewhat more complicated 
q Gaction wilt We to calculate the in the calculation of the terms of O**+, O*, and 
wave functions x) which minimize O By sad. Black,’ since Hertree’s 

method has been used instead of Fock’s, and 


the sum of expressions (3b) and (7a) and to 
calculate E,+E, corresponding to these wave 

| functions. Before doing this, however, an esti- 
q mate of the order of magnitude of the effect to 
f | be expected should be given. This can be taken 
from calculations of atomic spectra by the 
method of Fock’s equation or Hartree’s field, 


also because in the latter cases more than two 
electrons play important réles. For O** the 
differences between observed values (in brackets) 


and theory are 
§P(4.050)0.074; 'D(3.868)0.090; 


1$(3.658)0.176. 


| and their comparison with experimental results. 


The best result in this connection seems to be 
that on the normal state of He, where Fock’s 
equation is identical with Hartree’s. The dis- 
crepancy here is’ 0.077Ry (Rydberg units) for 
both electrons, or 12 Cal. per electron. This 
must be the amount of correlation energy in He. 


7™Cf. D. R. Hartree and A. L. Ingman, Mem. of the 


Manchester Lit. and Phil. Soc. 77, 69, 87 (1933). 


If we denote the radial wave function for the 
electrons by P(r), these terms correspond to the 
linear combinations P(r)P(r’) multiplied by 
(x-+iy)(x’+iy’); xx’ +yy' +22’. 
In the first case the correlation energy is very 
small, since the electrons are probably far away 


*D. R. Hartree and M. M. Black, Proc. Roy. Soc. 
A139, 311 (1933). 
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from each other anyway, because their spins are 
parallel. In such cases in the theory of metals 
we do not take into account any correlation 
energy at all. The discrepancy in this case may, 
to a considerable extent, be due to the use of 
Hartree’s instead of Fock’s method, and cer- 
tainly would be further diminished by the use 
of the latter. On the other hand, in the case of 
singlet’s, Hartree’s equations correspond much 
more closely to Fock’s equation and the increased 
discrepancy of about 0-050Ry is probably due 
to the neglect of the correlation energy. It is 
not quite clear, however, why it is so much 
greater in the 'S than in the 'D term. 

A comparison of experimental and theoretical 
values in O* and O points in a similar direction, 
the correlation energy is smaller though in these 
cases by a factor of the order 2. It is evident 
that it must diminish eventually if one goes 
over to more and more loosely bound electrons, 
since because of the lower electron densities the 
total interaction energy diminishes and the 
correlation energy is only the non-appearance 
of part of this. ; 

If one goes over to a metal like Na, on first 
sight the effect could be expected (because of 
the low electron density) to be much smaller 
than in He, about as great as in O. This will 
not be quite so, however, because the fluctuations 
in the potential of the electrons with downward 
spin will be greatly increased by their great 
number. The effect to be calculated in the next 
section will be about equal, therefore, to the 
effect per electron in He. 


4. 


Although the actual wave functions y¥,(x) in 
the Na lattice are actually different from plane 
waves ¢****-*/© in the Hartree-Fock approxima- 
tion, we shall use plane waves for ¥,(x) in the 
subsequent calculation. Since only integrals over 
the unperturbed functions occur in the perturba- 
tion calculation, this will not introduce a great 
error, because the ¥,(x) are extremely near to 
plane waves in much the greatest part of the 
volume. 

The whole following calculation will be per- 
formed in the approximation which corresponds 
to the second approximation in the Rayleigh- 
Schrédinger perturbation theory. The wave 
functions will be developed into a Fourier series 


2) = + (9) 


er'-2/© being taken as unperturbed function and 
***, *) the perturbation. The a,,, 
which are functions of the y, are supposed to be 
small, so that third order terms of a and V will 
be neglected for the energy, and second order 
terms for the wave function. 

For the actual calculation it can be seen, first 
of all, that one can replace a set of ¥,(y; x) by 
any orthogonal linear combination of them with- 
out affecting the final result. The orthogonality 
condition between y,(y; x) and y,-(y; x) gives 


yy * +L a, = 0. (10) 


By Schmidt’s method one can build a set of 
orthogonal y¥,(y; x) such that a,,,.=0 for 
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not be zero for all values of the y, but they 
evidently will not change the wave functions 
and merely give a contribution to the energy. 
The average value over all configurations of this 
contribution is zero, however, since the mean 
value of the potential of different configurations 
of charges is equal to the potential of the mean 
charges. Therefore, we shall set V,,(¥1, X) 
=0. For we have 


if v’ is occupied. It then follows under omission 
of the last term in (10) that a,,-=0 if » and v’ 
are both occupied states, whence the summation 
over uw in (9) must be extended only over the 
unoccupied states. It is then seen at once that 
(4a) is satisfied in our approximation. 

In order to calculate ¢,(y;---¥,) one must first 
calculate the matrix elements of V(y:, ¥n3%)- 
The matrix elements «++, ¥n3 = Voo will 


= f V(y; x)dx (11) 


cml 
because of Poisson’s equation. 
For ¢,, (3b) yields 


€(¥1" Vn) = Vip Vig*) +L View tL (u?—p-v) | (12) 


» mL? 


In order to calculate the ¢,’(¥,---y,) by means of (7a), the charge distribution for ¥,(y; x) will first 
be found : 


Now g(r) is 1 minus the function of the ‘‘hole,”” which decreases rapidly with increasing r. The 1, 
inserted into (7a), simply gives the energy difference of a uniform charge distribution and that 


corresponding to y,(y; x), namely, 


Avy 
x R (14a) 
4n2(u—v)?/L? 


in which the terms higher than the second order in a are omitted, and R means that the real part 
of the following expression is to be taken. The second sum must be taken only over those y's for which 
2v—y is unoccupied, while u itself is an unoccupied state in a// summations. For the calculation of the 
other part of (7a), arising from the function of the hole, one sets x, +r for x2 in (7a), introduces (13), 


and performs the integration over x. The result : 
unocc. 
must be multiplied with me*/2r and the function of the hole g(r) —1, which must be taken from (8), 
and integrated over r. Setting « =2r(u—v)d/L, this yields 


1.063 3.75 
). (14b ) 


1 
2.56L' \« 1+02/2.56 (1+02/2.56)? (1+02/2.56)? 


Added to (14a) this gives with help of the relation 4r7,°/3= v9 


q 

| 
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2! 
3ixt 
1 0.415 1.465 


o?(1+02/2.56) (1-+0%/2.56)? (1-+02/2.56)* 


o*e'(c) is given graphically in Fig. 2. 

This formula for the increase of potential energy between electrons with upward spin is not 
exact and may be viewed as containing two parts: First (14a), the increase of the potential of the 
space charges, due to the less even charge distribution for ¥,(y; x) than for ¥,(x). This increase is . 
lowered by the second part (14b), caused by the greater efficiency of the Fermi hole in a non-uniform 
charge distribution. The first neglection was made in setting ¥,(y; x) =@"*%-”)*/4),(y; x) when 
calculating ¢,’. This tends to increase the ¢,’ especially for those v, for which it is large anyway, 
because it overemphasizes the unevenness in the charge distribution. The second neglection was to 
keep no terms higher than the second order terms in a. This certainly decreases ¢,’, because part 
of the uneveness in the charge density is due to the higher terms, especially to those of the fourth 
order. Finally, the normalization constant, which is smaller than 1, enters in the second power in 
(14) and only in the first in (12). Its omission again increases (14). In the whole these errors will 
about compensate. 

The final quantity to be minimized is, after omission of the higher order terms, 


where 
4x*h? 2) e s2a(u—v)d 
t= (u?—yp-v), t,, =—— (15a) 
mL? 3ixtr, L 


V,, is given in (11), e’(¢) in (14), v is an occupied state, the last term in (15) should be taken only 
if 2»—y is an unoccupied state. By setting the derivative of (15) with respect to a,,* equal to zero, 


one obtains 


if 2v—y is occupied. If it is unoccupied 
V+ (toy +t,,’) yy = 0, + = 0. (16b) 


The last equation is obtained by differentiating (15) with respect to a,,_, and considering that 
= V>,* and =t,,’. Solving (16), one finds if is occupied 


if 2v—y is unoccupied, while yu is, of course, always unoccupied. These formulas show that the 
xX) do have the form (4) with 


and 


e 
(y—x) = —— +z ). (18) 
Inserting (17) into (15) one obtains for the total energy 
e+e’ =— (19) 


occ. unoce. 
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Instead of the second term one could write half of the sum of the terms for u and 2v—4y, which makes 
it somewhat more symmetric. 
5. 
It would be rather difficult to perform the summation over yu in (19) for an arbitrary set of y:, 
+++, Ya. Fortunately only the mean value of (19) with the weight | y|?/m! is needed and this can be 
computed quite easily. Since the ¢ do not depend on the y, one finds 


4 n 


(u—v)tL? 


1 

— 
n: 

(20) 


-_—_(z i-> 
apt 


(u—v)tL?\ 


The first term comes from «=i, the second from «#\. In the second, the summation over \ can be 
carried out, and it yields 1 if »y—-y+y», is an occupied state, zero otherwise. The whole bracket is, 
therefore, equal to the number of occupied states v’, for which »+»’ —y is not occupied. Fig. 3 shows 
a cross section of the v’-space through the origin and the point »—y. The radius of the circles is 
(3n/4r)*=L/2xd. The sphere through the weak circle contains the points y—y+v’ and the hatched 
part is unoccupied. As a consequence, it is 

e'n 1 for |a| >2 


1 
n(2x(u—v)d/L), = 21 
for |o| <2 
4 1 
and hence 
1 2R )d 2R) 
n! 3x3 3x (a) 
|o+p|>1 (22) 
|o—p|>1 


Here p=2nvd/L, o=2x(u—v)d/L and the summation over uw has been replaced by an integration 
over oa. The total energy is expressed in Rydberg units. The constant c is 


c= (e®md?/h?r,) = 0.1106r,/ao (22a) 


when expressed as function of the radius of the “‘s-sphere’’ in Bohr units ao. 


Fic. 3. Cross section of »’ space. Fic. 4. Section of o space. 


| 
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e ree 

d 

a 

0.1 

Fic. 5. Note: (1—p*)* should read (1 —p*)!. Fic. 6. 


In order to perform the integration of (22), one may first introduce elliptic coordinates o; = |¢|, 
o:=|o+p|=|2rud/L| in the o-space, using as centers the origin and the point —p (cf. Fig. 4). 
The first integral is to be extended over the horizontally, the second over the obliquely hatched 
region. The first will be a sum of two integrals (a) and (b) in Fig. 5, in both of which 9(¢) =30,/4 
—o,*/16. The second integral is a sum of three integrals (c), (d) and (e), and in the first two of 
them the same expression is valid for n(¢) while n(¢) =1 in the last. For de one writes 2r0,do,09d02/ p, 
of for and p= The integration over can be carried out simply and gives the 
five integrals 


p J 2ce'(o1) +301? + 4 16 
— +e 2o0+0;7+p?—1 


where v stands for (2¢e’(o,) +¢;*)!. A calculation of this quantity for r, =4 shows that it is practically 
equal o; if o,;>2 and the last integral can be evaluated accordingly. It yields 


8 2p? 2—p 2p) 9x? 20 


The other integrals were evaluated numerically and the results are plotted against p. Fig. 6 shows 
the plots for r, = 1, 2, 4, 8. From this, the mean value of the correlation energy which is the mean value 

—4F, with the weight p*, was calculated, and plotted against r,. The } enters, because the whole 
energy correction is present only for half of the electrons, that is, those with upward spin. In Fig. 7 
the upper curve represents the values calculated in this way. The energy is given in multiplés of ¢*/r,. 


1010 E. WIGNER 


6. 

One must remember, of course, that the 
preceding calculation is only an approximate 
one. Even if one confines oneself to wave 
functions of the form (2), the upper curve of 
Fig. 7 gives the correlation energy in first 
approximation only. The neglections are due to 
three causes: first to the use of an unnormalized 
wave function, second to the neglection of the 
terms with higher than the second power of a 
in ¢,, when going over from (12) to (14a) and 
third to the non-complete orthogonality of the 
wave functions employed. Our approximation 
is good if the a’s are small. One gets an idea 
about the accuracy of the approximation by 
calculating (which is the 
reciprocal square of the normalization constant 
minus one), though an idea only, since >, | a,,/* 
should really be small for all configurations of 
the y, not only its mean value. A calculation of 
the former quantity shows® that it stays well 
below one, except for large r, and for p which 
are very near to 1. It is in these cases that our 
approximation must be expected to break down. 

The real value of the correlation energy will 
be smaller in these cases than the calculated one. 
The correction with the normalization constant 
could easily be taken into account, as has been 
shown at another place.'’ It always decreases 
the correlation energy, not very much, however, 
as the magnitude of the normalization constant 
or the formulas in'® show. The second neglection 
is probably more dangerous and also more 
laborious to correct. It has been done for one 
point (r,=4) only and for this one very roughly. 
The second neglection also increased the calcu- 
lated value of the correlation energy, since it 
amounts to taking 1—2/ instead of (1—/f)? for 
the probability of the electron being at a certain 
point, and the minima of 1—2f are much 
deeper than those of (1—2f)* (the maxima are 
lower but less important). 


® The greatest part of the numerical work has been done 
by Dr. M. Vermes of Budapest. A table of the calculated 
values is given here: 


0 0-006 0-10 0-45 
0-4 0-009 0-14 0-54 
i 0-04 0-30 0-94 


1° To appear shortly in the Bull. of the Hung. Acad. 


If the electrons had no kinetic energy, they 
would settle in configurations which correspond 
to the absolute minima of the potential energy. 
These are closed-packed lattice configurations, 
with energies very near to that of the body- 
centered lattice. Here, every electron is very 
nearly surrounded with a spherical hole of 
radius r, and the potential energy is smaller 
than in the random configuration by the amount 
0.75 =e*/r,. This would be the sum of the correla- 
tion energy and that due to the Fermi hole. Since 
the latter one is," 0.458e?/r,, the maximum 
amount of the correlation energy is 0.292e?/r,. 
This value will be attained only if the kinetic 
energy can be neglected, i.e., for r,=~, and 
represents the asymptote to the real correlation 
energy curve, which is attempted to be drawn 
into Fig. 7. It appears to run much higher than 
one would have thought without calculation. 
I believe it to be in error everywhere by less than 
20 percent. 

The dotted line at 0.142e?/r, corresponds to a 
correlation energy as great as assumed in the 
first calculation® giving 0.6e?/r, together with 
the energy of the Fermi hole. 

The calculated constants of the Na lattice with 
the correlation energy of Fig. 7 and the other 
quantities as in reference 1 are as follows: lattice 
constant 4.62A as compared with the observed 
value of 4.23A. The binding energy associated 
with this is 26.1 Calories, to be compared with 
the observed value of 26.9. The calculated value 
of the binding energy for the observed lattice 


Fic. 7. 


" F. Bloch, Zeits. f. Physik 57, 545 (1929). 
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constant is 22.3 Cal. As far as the lattice constant 
goes, one must remember, however, that both 
the correlation energy and that due to the Fermi 
hole are calculated for a flat wave function and 
the wave functions are flat only for r,=4 and 
its neighborhood. 

The magnitude of the correlation energy is 
important for questions of paramagnetism and 
ferromagnetism as well as for questions of lattice 
energy. It modifies Bloch’s original theory on 
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the ferromagnetism of free electrons" in such a 
way that it yields ferromagnetism in fewer cases 
than in its original form.” I hope to return to 
this question at another time. 

I wish to express my gratitude to Dr. F. Seitz 
for his kind help in connection with the prepara- 
tion of this manuscript. 


'? A paper of S. Schubin and S. Wonsowsky, Proc. Roy. 
Soc. A145, 159 (1934) which appeared recently, points in 
the same direction. 
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The coefficient of thermal expansion is measured at intervals of 2.5°C between 200°C and 
500°C for single and polycrystalline specimens of nickel of varying known degrees of purity. 
The data yield the value 3.24+0.1510~‘ for the ferromagnetic change in volume per unit 


volume of pure nickel. 


INTRODUCTION 


HEN the coefficient of thermal expansion 

of nickel is plotted as a function of tem- 
perature the resulting graph has a hump which 
starts at about 200°C and terminates in the 
neighborhood of the Curie point. The area under 
the graph represents a change of length (and 
hence of volume), and the area under the hump 
may properly be taken as a measure of the 
change in volume associated with those inter- 
atomic forces (or energies) in terms of which the 
ferromagnetism of the material finds its expla- 
nation. Fowler and Kapitza' were the first to 
point out that Heisenberg’s theory of ferromag- 
netism is competent to offer a quantitative 
description of this phenomenon. Their calculation 
has been extended by Powell,? who obtained an 
expression relating the change in volume per unit 
volume to the exchange energy between pairs of 
electrons belonging to neighboring atoms. The 
observations upon which Powell based his nu- 
merical estimate of the former quantity were 


1 Fowler and =o, Proc. Roy. Soc. Al24, 1 (1929). 
? Powell, Proc. P 


ys. Soc. 42, 390 (1930), 


obtained by Colby,’ who worked with poly- 
crystalline nickel of unspecified purity. It seemed 
worth while to repeat Colby’s measurements 
upon single and polycrystalline nickel of varying 
known degrees of purity, and to extend them 
over a greater range of temperature in order to 
increase the precision of the base line from which 
the hump is reckoned. The present paper is a 
report based upon these experiments. 


APPARATUS AND METHOD 


The dilatometer (Fig. 1) is constructed entirely 
of fused quartz. A knife edge, A, and a table, B, 
rest upon a flat plate, C. The top of the table is 
ground flat and polished, and carries a ground 
roller about 1 mm in diameter. The specimen, 5, 
is a cireular cylinder about 5 mm in diameter and 
6 cm to 7 cm long. A tiny lateral scratch is made 
near one end and the other end is polished. The 
scratch engages the knife edge and the pclished 
end rests on the roller. Mirrors M; and M,, of 
gold sputtered on quartz, are fused to the ends of 
the roller. These mirrors are rotated slightly with 


* Colby, Phys. Rev, 30, 506 (1910), 


A B 


Fic. 1. The quartz dilatometer. 


Fic. 2. The dilatometer and furnace assembly. 


respect to each other about the roller axis to 
increase the range of the instrument. A mirror, 
M,, is fused to the knife edge, and another, M;, 
to the table. 

Fig. 2 is a diagram of the dilatometer and 

furnace assembly. The dilatometer rests in the 
center of the furnace on a quartz plate fused to 
one end of a quartz tube 2.5 inches in diameter 
and 2 feet long. The other end of this tube is held 
in a mounting insulated against building vibra- 
tion. The heating element of the furnace is a 
bifilar winding of No. 18 Constantan wire on a 
quartz tube 3 inches in diameter and 26 inches 
long, together with an auxiliary winding over a 
length of 6 inches at each end. The open end of 
the furnace is closed with a microscope cover 
glass. . 
The junction of a platinum vs. platinum—10 
percent rhodium thermocouple nearly touches 
the specimen at its center. The thermal e.m.f. is 
measured with a Leeds and Northrup Type K 
potentiometer. 

The mirrors are so aligned that they reflect 
into a theodolite light from four vertical scales, 
one for each mirror. The scales are four meters 
from the dilatometer, and can be independently 
illuminated. Readings on the scales associated 
with mirrors M, and M; reveal any rotation of 
the dilatometer as a whole. 

Prior to each set of observations the specimen 
is demagnetized in an alternating magnetic field, 
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mounted in the dilatometer, and carried through 
several cycles of heating and cooling over the 
temperature range 300°C to 400°C. Thermal 
equilibrium is then established and observations 
are made at 5° intervals over a complete cycle of 
heating and cooling. The initial temperature is 
then altered 2.5 degrees and the procedure re- 
peated. The data thus yield mean values of the 
coefficient of thermal expansion over 5° intervals 
corresponding to mean temperatures which 
differ by 2.5°. 


RESULTS 


Observations are here reported upon five dif- 
ferent specimens of nickel, namely: 


I. A single crystal of Mond nickel grown in an atmos- 
phere of hydrogen at a pressure of about 1 mm;* 

II. A single crystal grown of nickel purified in the 
manner described by Fink and Rohrman,* but contami- 
nated with carbon, and with molybdenum from the winding 
of the furnace in which the crystal is grown; 

III. A polycrystalline specimen of the same material as 
II, but containing more molybdenum; 

IV. A specimen of hard drawn commercial nickel; 

V. Specimen IV, annealed 1 hr. in hydrogen at 1100°C 
and cooled slowly. 


Table I shows the purity of the several speci- 
mens. The analysis of specimen IV was supplied 
by the International Nickel Co. The carbon 
content of specimen III was determined chemi- 
cally by Dr. L. A. Wooten, of the Bell Telephone 
Laboratories. The remainder of the data repre- 
sent quantitative spectroscopic analyses made by 
Lucius Pitkin, Inc. of New York. 

Laue photographs taken through thin sections 
sawn from the ends of specimens I and II 
indicate that the orientations of the respective 


TABLE I. Percent impurily in specimens. 


I II Ill IV 


Mo 0.075 0.5-0.6 0.75-0.9 0.11 
0.02-0.06 0.06-0.1 0.05-0.09 0.3-0.4 
Cu 0.01-0.02 0.005-0.01 0.00 1-0.005 0.41 
Fe 0.001-0.005 0.005-0.01 0.01-0.02 
Mg trace trace trace 
Al 
Si trace trace trace 0.09 
Sn s trace trace 

* Not found. 


‘ Dingwall, Zacharias and Siegel, Trans. Electrochem. 


Soc. 63, 395 (1933). 
5 Fink and Rohrman, Trans. Electrochem. Soc. 57, 


325 (1930); 58, 403 (1930) ; 59, 359 (1931). 
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Fic. 3. Variation of coefficient of thermal fre gy with temperature for various specimens 
nickel. 


cylinder axes relative to the crystal axes are 
quite dissimilar. 

The data obtained with the several specimens 
are shown graphically in Fig. 3. The crosses on 
curve IV represent observations made on speci- 
men V. Each point in the range 320°-380° 
represents the mean of four to six observations. 
All the observations are corrected for the ex- 
pansion of the quartz dilatometer. 

It is to be noted that there is no evidence of 
hysteresis in this phenomenon; the heating and 
cooling curves are indistinguishable. Again, the 
annealing process to which specimen IV was 
subjected leaves the phenomenon unaltered. 
Rough measurements of the magnetization as a 
function of temperature in specimen I indicate 
that the Curie point of this material is in the 
neighborhood of 370°C. 


PRECISION 


The sources of systematic error in these meas- 
urements are (1), errors in measurement of the 
specimen length and of the quantities which 
determine the dilatometer magnification; (2), a 


longitudinal temperature gradient in the furnace ; 
(3), the thermocouple calibration; and (4), a 
difference in temperature between the thermo- 
couple junction and the specimen. 

The distance between the roller and knife edge 
is measured with a travelling microscope to 0.1 
percent, and the dilatometer magnification to 0.2 
percent. The maximum allowed temperature 
gradient near the center of the furnace,measured 
with the specimen removed, is 0.13°C per cm. 
Accordingly, the actual temperature of the 
specimen must be uniform to better than 0.5°C. 
The thermocouple calibration used is that given 
by the Bureau of Standards.* This is checked and 
the remaining sources of error estimated by 
replacing the specimen with a glass encased 
Leeds and Northrup platinum resistance ther- 
mometer calibrated by the Bureau of Standards. 
Under these circumstances the maximum dis- 
crepancy between the two thermometers over a 
period of eighteen months was 0.8°C in the range 
300°-400°. Furthermore, the 5° intervals as 


* Bur. Standards J. Research 10, 280 (1933). 


| 


1014 
i 
g 320° 340° 360° 380° 
“Temperature 


Fic. 4. Illustrating the behavior of the dilatometer. The 
substance is fine silver. 


measured by the two thermometers agree to 
better than 0.09°C. 

The magnitude of the accidental errors inher- 
ent in the method are exhibited in Fig. 4, which 
depicts the result of a series of measurements on 
a bar of fine silver. The lengths of the lines in- 
dicate twice the mean deviation of four sets of 
observations. The crosses show the mean values 
of two sets of observations taken after the 
dilatometer had been completely dismantled and 
reassembled. 


THE FERROMAGNETIC CHANGE OF 
VOLUME 


An evaluation of the ferromagnetic change in 
length requires an authentic base line from which 
to calculate the area under the hump. This is 
secured by combining the data obtained in these 
experiments with those reported by Hidnert.’ 
Hidnert worked with a specimen of electrolytic 
nickel which contained 0.03 percent iron and 0.03 
percent of other metallic impurity, and reported 
mean values of the thermal expansion coefficient 
over temperature intervals of 100°C between 
100°C and 900°C. 

Fig. 5 is a graph of the composite data. The 
crosses indicate Hidnert’s values and the circles 
in the temperature range 400°-500° indicate 
values obtained in these experiments. The 
humps are reproductions, on the altered scale, 
of curves I and II of Fig. 3. With the exception 
of the single cbservation of Hidnert over the 
range 400°-500°, the data are in excellent accord, 
and yield a fairly satisfactory base line shown by 
the dotted portion of the curve. It appears that 
0.6 percent or less of impurity leaves the base line 


* Hidnert, Bur. Standards J. Research 5, 1305 (1930), 
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Fic. 5. The base line for measuring the ferromagnetic 


on in length. The humps are curves I and II of 
ig. 3. 
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Fic. 6. Variation of the conse pee change in length 
with total impurity for nickel. 


and ascending part of the hump unaltered, and 
that the value of the change in length per unit 
length diminishes with increasing amount of 
impurity. Accordingly it is suggested that the 
value of the ferromagnetic change in length for 
pure nickel be obtained by extrapolation. 

Fig. 6 is a curve for which the ordinates of the 
plotted points are the areas under the four curves 
of Fig. 3 with respect to the base line shown in 
Fig. 5, while the abscissae are the percentages of 
total impurity as given in Table I. The extra- 
polated value of the ferromagnetic change in 
length per unit length for pure nickel is 


1.08+0.05 10~. 


The precision measure is obtained from the un- 
certainty in the total impurity as given in Table I. 

In conclusion the author desires to express his 
thanks to the Physics Department of Columbia 
University for the facilities placed at his disposal 
to Dr. L. A. Wooten for the determination of the 
carbon content of specimen III, and to Dr. S. L. 
Quimby, who suggested the problem and fol- 
lowed the progress of its solution with helpful 
counsel and encouragement. 
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A Survey of Time Lag of Sparkover in a Uniform Field 


Ase TILLEs, Department of Electrical Engineering, University of California 
(Received July 14, 1934) 


An original instrument for measuring time lag of spark- 
over has been developed. This is capable of measuring lags 
over a range from one one-hundred-thousandth of a second 
up to any longer values. This entire time range is available 
for any individual measurement without readjustments 
and independently of the method or wave shape of voltage 
application. Spark lags are found to fall into three distinct 
domains. First,.long initiatory lags (10** to 10~* second or 
less) of random distribution spent in awaiting a favorable 
fortuitous space and time arrangement of the somewhat 
insufficient supply of initiatory electrons. Second, inter- 


mediate formative lags (of the order 10~* second) dis- 
tributed about a definite peak value of the order of time 
required for a positive ion to cross the gap and assist in 
building up a space charge gradient. Third, short formative 
lags (of the order 10-7 second) distributed about a definite 
peak value of the order of time required for an electron 
to cross the gap. A picture of the mechanism of sparkover 
in an initially uniform field is presented accounting for 
the existence of these three domains. The available experi- 
mental reports are correlated and shown to be consistent 
with each other and with the proposed explanation. 


INTRODUCTION 


S is well known the theory of sparkover has 
long been based on the original formu- 
lations of Townsend. Increasing knowledge 
brought to light serious difficulties in this formu- 
lation. First, at the field strengths envisaged, the 
positive ion could not possibly attain an energy 
requisite to produce ionization. Second, because 
of the difference in electron and positive ion 
mobilities, a space charge is produced which 
must render the field at the time of sparkover 
decidedly non-uniform. From theoretical con- 
sideration of these factors, Rogowski' and Loeb? 
independently concluded that the Townsend 
mechanism did occur but at field strengths 
greatly multiplied locally by distortion caused 
by the space charge. They pointed out that this 
modified Townsend mechanism demanded a 
definite formative time for the sparkover of the 
order of time required for a positive ion to cross 
the gap, e.g., 10~* second, and that a measure- 
ment of the time lag of sparkover under con- 
ditions which eliminated the random variations 
found by Zuber*® would therefore provide a 
crucial verification of the conclusions reached. 
Shortly thereafter much shorter time lags, e.g., 
10-* second were observed for impulse sparkover 
at substantial overvoltages by Beams,‘ Tamm,’ 


et al. This apparent direct contradiction left the 


subject in confusion. 


* Rogowski, Archiv f. Elektrotechnik 16, 761 (1926). 
? Loeb, J. Frank. Inst. 205, 305 (1928). 

5 Zuber, Ann. d. Physik 76, 231 (1925). 

* Beams, J. Frank. Inst. 206, 809 (1928). 

* Tamm, Archiv f. Elektrotechnik 19, 235 (1928). 


Because of the possible variety of phenomena 
occurring under different conditions this diver- 
gence did not in fact invalidate the above men- 
tioned theory. This was specifically pointed out 
by Loeb in an ensuing paper* and in a seminar 
which it was the author's good fortune to attend 
in 1931. Loeb pointed out that it still remained 
necessary to check the theory by means of a 
time lag measurement under those conditions of 
stati¢ sparkover for which the theory had been 
formulated, namely, small overvoltage and initial 
ionization adequate to eliminate random vari- 
ations. 

It was, accordingly, the object of this present 
work first to ascertain under which, conditions 
the spark lag distribution was random and ander 
which the distribution curve was peaked about 
some definite value. Secondly, it was desired to 
measure the time lag for the important case, 
hitherto unexplored, of adequate initial ioniza- 
tion and small overvoltage. Third, by means of 
these results it was desired, if possible, to corre- 
late the available data on spark lag and on the 
nature of its dependence on initial ionization and 
on overvoltage and thereby, perhaps, to clarify 
the mechanism of sparkover in a uniform field. 


APPARATUS: AND PROCEDURE 


The timing instrument evolved used common 
three electrode vacuum tubes in a resistance- 
coupled amplifier so biased that a small in- 
crement of gap voltage at minimum sparking 


* Loéb, Phys. Rev. 38, 1891 (1931). . 
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Fic. 1. An instrument to measure time lag of sparkover. 
K, kenetron rectifier tube; S, safety gap and discharge 
switch; V, électrostatic voltmeter; C, storage condenser 
7.2X 107" f; G, test gap; Q, quartz mercury arc; A, bias 
adjustment; 7, synchronous timer; B, ballistic galvanom- 
eters. 


Plate Plate 
battery circuit 
Tube voltage resistor 
10M 5 156 1.92 MQ 
2CX 340 180 0S MQ 
3 CX 345 180 0.1 MQ 
4CX 340 90 0.25 MQ 
5 CX 345 180 500 


voltage was steeply amplified and applied to a 
tube whose output current was thereby swept 
from zero to a constant maximum value. This 
constant current flowed until the sparkover 
occurred. The current was passed through a 
ballistic galvanometer whose throw was directly 
calibrated in seconds, thus measuring the spark 
lag. 

The wiring of the instrument is fully shown in 
Fig. 1. Tube 1 is directly across the spark gap, 
acting as an inverted vacuum tube voltmeter 
and reducing the voltage to a conveniently small 
value. This reduced voltage in series with an 
adjustable bias is applied to the input of tube 2. 
Tubes 2 and 3 act as a two-stage, resistance- 
coupled, linear amplifier. The output of tube 3 
applied to the grid of tube 4 swings this grid 
from a substantial negative voltage to a sub- 
stantial positive voltage causing the output 
current of tube 4 to sweep from zero to the 
maximum obtainable (i.e., plate supply voltage 
divided by external plate resistance) for a small 
rise in voltage at the gap. The adjustable bias 
determines at which gap voltage increment this 
sharp response takes place and is set so that the 
minimum sparking voltage of the gap gives half 
maximum output current. It is evident that for 
the approximately trapezoidal wave shape of 


current the total charge measured by the galva- 
nometer divided by the known maximum value 
of current gives the time that the current has 
remained above half maximum value, hence the 
bias setting mentioned. The galvanometer throw 
is calibrated to give directly in seconds the time 
that the voltage at the gap has remained higher 
than the value for which the bias was set, i.e., 
higher than the minimum sparking voltage of the 
gap. Thus the galvanometer reads directly the 
time lag of sparkover of the gap. The output 
current is passed through two ballistic galva- 
nometers of different sensitivities in series. It is 
feasible manually to short circuit one or both 
of these after the beginning of its swing when 
necessary to prevent its being damaged by an 
excessive swing. Tube 3 also supplies the input 
of tube 5 which operates a synchronous timer 
from which any time intervals greater than 0.05 
second are easily read. Theoretically computed 
transients in simple series circuits were applied 
to the instrument and these gave a calibration 
check to better than 4 percent. The instrument 
reads times from 10~* second up. 

The timing circuit itself serves also for meas- 
urement of voltage. It is to be noted that the 
curves of output current of tubes 2, 3 and 4 may 
be calibrated and used to read the relative 
voltage at the gap, i.e., the difference between 
the voltage at the gap and the minimum sparking 
voltage for which the bias has been adjusted. 
Tube 2 allows readings over a wider range before 
saturation is reached but with a lesser precision 
than tube 3; similarly in comparing tube 3 to 
tube 4. In most of the tests the plate current of 
tube 2 was used to read directly with a precision 
of about 1/4 percent overvoltages less than 7 
percent of the sparking potential (which was 
near 4000 volts). 

When the resistances of the instrument were 
adjusted for maximum amplification a shift of 
six volts at the gap, i.e., 0.15 percent of 4000 volts, 
produced a change in the output current from 
zero to the maximum. Thus even if the voltage 
were raised slowly, e.g., at the rate of 100 
percent of sparking potential per 0.01 second, 
the time required for the output current to 
rise from zero to maximum would still be only 
1.5X10-* second. Thus the instrument can be 
readily used to measure lags of this or any larger 
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Fic. 2. Voltage wave shapes. a, “static” sparkover; 
b, “impulse” sparkover. (The time for the voltage changes 
drawn vertically is about 1/1000 of the smallest spark 
lags measured.) 


order even for such slow rates of voltage applica- 
tion and, of course, for any more rapid rates. 
(In the present test series the tube 4 plate re- 
sistance was readjusted to give a larger current 
than at the above mentioned setting, necessi- 
tating the use of larger voltage increments at the 
gap than mentioned above. There was no 
difficulty on this score, however, as in this test 
series increments of more than 100 volts were 
used throughout.) 

In order that the effect on the spark lag of the 
magnitude of the overvoltage might be clearly 
evident, it was decided to use a rectangular 
voltage wave shape. That is, the circuit was 
arranged so that the time duration of the voltage 
rise from minimum sparking potential to the con- 
stant maximum voltage applied to the gap 
throughout the time lag interval, and also the 
time duration of the drop of voltage from its 
full value to the minimum sparking potential, 
were each very much less than the time lag 
being measured. The voltage was swept past the 
sparking voltage very rapidly by applying all 
or an increment of the voltage suddenly, i.e., 
within some 10~* second. In most cases this was 
accomplished by applying the voltage to the gap 
by closing a switch having higher insulation 
resistance than the gap. An auxiliary spark was 
thus obtained at the closing switch which served, 
therefore, to apply the voltage to the gap very 
quickly. Two wave shapes were utilized (Fig. 2). 
For the first, which conveniently approximated 


“static” sparkover, the main d.c. voltage on the 
gap was set at 96 percent of sparking potential 
and a bank of “B” batteries was then thrown 
across a resistancé in the circuit to suddenly 
raise the gap voltage to sdme percent above 
sparking potential. The second wave shape, 
constant d.c. voltage suddenly applied, approxi- 
mated “‘impulse’’ sparkover but at small over- 
voltages. The d.c. voltage was obtained from a 
synchronously driven self-excited motor-gener- 
ator set, transformer, kenetron tube and con- 
denser. Time lags of sparkover were obtained for 
both these voltage shapes, for various voltages, 
mostly in the neighborhood of 105 percent of 
static sparking potential, and for ultraviolet 
illumination of measured intensity varied over a 
very large range. 

Ultraviolet light from a quartz mercury arc 
illuminated the test gap and also a plate of clean 
nickel in vacuum—both of these being illumi- 
nated through identical quartz windows. The 
intensity of illumination was measured by meas- 
uring the saturation current liberated from the 
nickel photo-cell. This was allowed to flow 
through a specially constructed resistor con- 
sisting of a micarta strip buried in paraffin; 
the voltage drop across the resistor was measured 
by a Dolazelek electrometer, the resistance being 
measured after each voltage measurement by 
means of the voltage time decay curve when the 
illumination was intercepted. The spark gap 
used consisted of copper spheres at small spacing 
in dehydrated air at atmospheric pressure. 

The coherence of the distribution curves in- 
creases decidedly as the number of sparkovers 
per run is increased, but so too, does the work in- 
volved. Accordingly, as a compromise, a run of 
50 points for the random distribution and of 25 
for the peaked was selected as a standard pro- 
cedure. It may be noted that the adherence of 
the points to the curve is always best at the top 
where the points are most numerous and poorest 
at the bottom where the scale is stretched and a 
small number of points is used to locate a 
statistical distribution. For this reason the curve 
is placed so that the number of points above it 
equals the number below, no weight being 
allowed to the added ‘‘moment”’ of the straggling 
points at the lower end of the curve. . 
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DATA AND RESULTS 


Fig. 3a shows three curves of spark lag dis- 
tribution at “static” sparkover at three different 
voltages and at one intensity of initial ionization. 
As labelled this intensity is measured by the 
current J, in micro-micro-amperes per square 
centimeter, ejected from a plate of clean nickel 
in vacuum by the ultraviolet illumination util- 
ized. The lower left hand corners of the light 
curves give the actual data points which are 
plotted by listing all the spark lags of the set 
in order of magnitude and listing also the percent 
order of magnitude. Thus the fourth largest of a 
set of 50 is 8 percent etc. This percentage is 
plotted against the time duration of the par- 
ticular lag. The equation of the resultant curve is 
n,=100e~?' wherein n, is the percentage of the 
total number of lags which are greater than the 
time ¢ and ) is the probability of a lag terminating 
within unit time. It is evident since the curves 
are linear on semi-logarithmic paper that / is 
constant and the distribution truly random. As 
is well known,’ the average time lag 7 =1/p and 
can be taken from the curve at ,= 36.8. 

The average time of these distributions is 
plotted against percent overvoltage in Fig. 30. 
From this curve the average time at 5 percent 
overvoltage is readily interpolated, thus giving a 
point on the master curve 1 of Fig. 5. Each of 
the curves similar to 3d, i.e., average time vs. 
overvoltage, was roughly exponential for over- 
voltages from 1.to 6 percent and the several 
exponents were approximately the same. 

As the illumination was further increased, 
however, a different transitional distribution 
appeared wherein some 70 percent of the lags 
formed a true random distribution and one, 
moreover, fitting the same master curve 1 of 
Fig. 5, as shown by its dotted portion. The other 
30 percent of the lags were grouped into an 
entirely different type of distribution. As the 
illumination was further increased larger per- 
centages of the total were in the second group, 
and for somewhat greater illuminations yet all 
the lags fell on the second type of distribution 
curve. This was a highly peaked distribution 
containing a few large and many small deviations 
from a definite peak value. In fact, some 40 


von Laue, Ann. d. Physik 76, 261 (1925). 
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Fic. 3. Domain one spark lags. a, spark lag distributions, 

rcent of lags greater than given by abscissa; b, average 
ag versus percent overvoltage. Approach voltage=96 
percent; overvoltage, curve 1=0.7 percent; overvoltage 
curve 2=3.5 percent; overvoltage, curve 3=4.9 percent ; 
J=0,00183. 
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Fic. 4. Spark lag distributions for different illuminations; 
a, random at J=0,00183; b, random at J=0.0258; c, 
transitional at J=0.365; d, peaked at J=5.75. 


percent of the total lags fell within a single time 
increment only 10 percent as large as the 
maximum time intervals represented in the dis- 
tribution. For the peaked distributions no appre- 
ciable shift in the time lag was produced by 
slight changes of overvoltage nor by extensive 
further increase in illumination. This also is 
shown on the curves of Fig. 5. Fig. 4 illustrates 
the two distributions and the transitional dis- 
tribution by curves at different illuminations and 
near 5 percent overvoltage. 

A series of tests similar to the above was re- 
peated for the “impulse” wave shape with very 
similar results which are also plotted on Fig. 5 as 
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Fic. 5. Master curve of average or peak values. Time 
lag of sparkover versus illumination 


Approach 
Curve voltage Voltage 
1 96% 105% 
2 96% 103°) 
3 0% 105% 


curve 3. Throughout the random distributions 
the lags at impulse sparkover are longer by a 
factor of two or more than those at static spark- 
over. Apparently the only difference is that the 
96 percent approach voltage, due to cumulative 
ionization, causes a greater concentration of 
initial ions to exist in the gap than exists with 
the 0 approach voltage. At the smallest ioniza- 
tions studied, i.e., at the upper end of the impulse 
sparkover curve 3 of Fig. 5, the curve takes on a 
greater slope. This is perhaps to be expected as 
at the lower ionization density one would expect 
a unit slope. That is, a tenfold increase in ioniza- 
tion should effect a tenfold shortening of the 
time lag. For higher initial ion densities the lesser 
slope observed indicates that the increasing 
ultraviolet illumination becomes in some way 
progressively less efficient in precipitating spark- 


over. 

It is apparent that if this 96 percent approach 
voltage is lowered, the random spark Jags will be 
lengthened and that at some low approach 
voltage, for which the cleansing action entirely 
outweighs the ionization caused, the random 


lags must reach a maximum duration greater 
than that at zero appréach voltage. This was 
checked by a brief test series at 50 percent 
approach voltage which gave nearly three times 
as great an average lag as that for zero approach 
voltage. A few brief tests were taken, also, with 
overvoltages as high as 50 percent. As was 
expected increasing the voltage, thus, sub- 
stantially decreased the observed spark lag. 
At 50 percent overvoltage, again, a peaked dis- 
tribution was obtained at a high illumination 
I=9, but a somewhat lower illumination, J 
=1.75, which had been adequate to give a 
peaked distribution at 5 percent overvoltage was 
inadequate to give a peaked distribution at the 
higher overvoltage. That is for J=1.75, at 5 
percent overvoltage the initiatory lags were 
negligible compared to the formative lags, but 
at 50 percent overvoltage the greatly shortened 
formative lags were negligible compared to the 
initiatory lags. 

The results of these tests may be summarized 
as follows : Spark lags are found to fall into three 
distinct domains. First, long imitiatory lags 
(10*+* to 10-* second or less) of random dis- 
tribution spent in awaiting a favorable fortuitous 
space and time arrangement of the somewhat 
insufficient supply of initiatory electrons. Second, 
intermediate formative lags (of the order 10~* 
second) distributed about a definite peak value 
of the order of time required for a positive ion 
to cross the gap and assist in building up a 
space charge gradient. Third, short formative 
lags (of the order 10-7 second) distributed about 
a definite peak value of the order of time required 
for an electron to cross the gap. 

In the first domain with small overvoltages for 
the gap studied (copper spheres of 0.952 cm 
radius at 0.0683 cm separation, and 3820 volts 
static sparking potential) the average time lag 
is given by: 

t=f(v)I-*-*, (1) 


where t= average time lag in seconds of a group of spark- 
overs taken under identical conditions. 
I=intensity of photoelecric illumination as meas- 
ured by the current, in micro-micro-amperes per 
square centimeter, exepelled from a clean nickel 
plate in vacuum. 
f(v) =a function of the applied voltage wave shape and 
magnitude. 
Although f(e) varies somewhat with the illumination, it 
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may be inclusively stated for the ‘‘static’’ sparkover (i.e., 
voltage suddenly shifted from 96 percent of static sparking 
potential to a constant value more than 100 percent) as: 
f(v) =0.0037¢-0 0.8)» (2) 
and for the “impulse” sparkover (i.e., a constant maximum 
voltage suddenly applied) as: 
=0.0078e— © #0 43)" (3) 
where v=the percent overvoltage 
max. applied voltage —static sparking voltage 
static sparking voltage 


100. 


Eqs. (2) and (3) hold for from 1 to 6 percent 
overvoltage ; outside of this range the function is 
not exponential. The central value cf exponent 
indicated is correct at 5 percent overvoltage in 
every case. For static sparkover Eq. (1) holds 
for illumination from J=0.002, the lowest 
studied, to / =0.5 where transition to the second 
domain occurs. For impulse sparkover the curve 
deviates slightly, the exponent becoming greater 
as I decreases; Eq. (1) is adequate, however, 
from J=0.02 to J=0.5. 

At the transition to domain two, ie., at 
I=0.5, 70 percent of the sparkovers occur in a 
distribution whose average exactly fits equation 
(1) of domain one and 30 percent fit into the 
domain two distribution. For successively larger 
illuminations a successively sma!ler percentage 
of sparkovers come into domain one and a 
larger percentage into domain two. 

The existence of the second domain, which had 
been anticipated by Rogowski,' and specifically 
predicted by Loeb,?:* is now experimentally 
verified for the first time. From J=0.5 to J=10, 
the highest studied, for this gap and either of 
the voltage wave shapes mentioned and for over- 
voltages of from 2 to 6 percent, the lags have a 
highly peaked distribution, some 40 percent of 
them falling within one 10 percent increment 
of the maximum time intervals encountered. 
Throughout this region the formative time lag, 
as given by the peak of the distribution curve, 
remains at (0.70+0.15) x 10~ second. The above 
mentioned prediction*® applied to the present gap 
length calls for a formative time lag of 0.5 x 10~ 
second, an excellent agreement. 

When the voltage is raised to 150 percent of 
static sparking potential at J=9, the lags be- 
come less than 1.5X10-* second, this being 
about the smallest interval that could be meas- 


ured. At 150 pércent voltage and J=1.75 a 
fraction of the Jags are larger, falling in domain 
one and illustrating roughly the transition from 
domain one directly to domain three at higher 
voltages. This transition, as well as the existence 
of the domain one distribution for quite brief 
spark lags, is more comprehensively illustrated 
by the excellent experimental work of Strigel,*: * 
although his interpretation of domain three is 
confused. Domain three at higher overvoltages 
has received considerable investigation by means 
of the cathode ray oscillograph from Rogowski's 
students’: and from various engineers” 
studying impulse sparkover at high overvoltages. 
It has also been studied optically.‘ 

The time of propagation of the breakdown 
path and the completion of sparkover once it 
has been adequately produced involves still 
shorter time intervals."* These may add appre- 
ciably to time lag only in the case of long gaps 
(i.e., initially non-uniform field) and do not 
concern us here. 


SPARKOVER MECHANISM 


The mechanism of sparkover in an initially 
uniform field is envisaged as follows. 


Domain one 

With somewhat inadequate initial ionization 
the imitiatory lag in awaiting the fortuitous 
occurrence of electrons so located as to be 
capable of commencing sparkover is very much 
greater than the formative lag. In general, since 
this fortuitous occurrence is as likely within one 
time interval as within any other, distribution of 
these initiatory lags is truly random. The order 
of magnitude of the initiatory lag depends on 
the external ionizing forces and may be extremely 
large if these are small. In general, with a 
probable exception at voltages barely large 
enough to produce breakdown upon the appear- 
ance of a single electron at the cathode, it is the 
density of ionization in the gap that controls the 
initiatory lag rather than the rate of production 
of mew electrons. As the initial ion density is 


*Strigel, Wissenschaftliche Veroffentlichungen a.d. 
Siemens-Konzern 11, 52 (1932). 
* Strigel, Archiv f. Elektrotechnik 27, 137 (1933). 
(193 — and Masch, Archiv f. Elektrotechnik 25, 787 
1). 
4 Viehmann, Archiv f. Elektrotechnik 25, 253 (1931). 
1 Torok, A. I. E. E. Trans. 47, 349 (1928). 
144 Dunnington, Phys. Rev. 38, 1535 (1931). 
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sufficiently increased the entire initiatory lag 
becomes negligibly small and the formative lag 
of either domain two or three as the case may be, 
prevails. 

At one level of ionization density we may 
visualize the variation of initiatory lag with 
voltage variation as follows. The initiatory time 
lag is the random time spent in awaiting the 
drrival of the suitable initiatory electrons within 


a specific effective vclume in the gap. The extent 


of this volume effective for the initiation of 
adequate electron avalanches obviously increases 
with increased voltage. Thus for sparkover at 
the very lowest overvoltages, e.g., 0.1 percent, 
it is certainly necessary for the initiatory elec- 
trons to appear at the very surface of the cathode 
and it is also necessary—in all probability—for 
several to appear almost simultaneously, each 
one taking advantage of the gradient produced 
by the space charge of its predecessor. For some- 
what higher voltages it is certainly sufficient for 
sparking that the initiatory electrons appear 
anywhere within an “effective” shell of appropri- 
ate depth extending from the cathode; also, a 
smaller number of electrons probably suffices. 
As the voltage increases the thickness of this 
effective shell increases. For very high voltages 
certainly it is sufficient that a single initiatory 
electron exist anywhere within an effective vol- 
ume including all the gap except a thin shell at 
the anode. Further voltage increases will then 
have little effect in shortening this initiatory lag. 


Domain two 

For a sufficient initial ionization the formative 
lag prevails. Each initial electron avalanche 
leaves behind it an increasing phalanx of positive 
ions pointed toward the cathode and hence a 
space potential which increases about expo- 
nentially as one leaves the cathode and starts 
towards the anode. The slope of this curve gives 
the field strength at each point. Using low over- 
voltages, e.g., 3 percent and a small gap length, 
this field strength is at no point great enough at 
the first instance to cause positive ion ionization. 
As the whole pialanx of positive ions, however, 
proceeds relatively slowly toward and into the 
cathode, the resultant field: between the cathode 
and the next row of positive ions becomes 
successively steeper. As the last row of positive 
ions approaches the cathode, the field there 


becomes great enough to cause positive ion 
ionization giving a tremendous self-accelerating 
increase in the number of new electfons formed. 
As soon as this burst of electrons reaches the 
anode, the voltage begins its complete drop and 
sparkover is very rapidly completed. The total 


elapsed time is about that required for a positive 


ion to cross the gap. 

With slightly higher overvoltages (e.g., 15 
percent) since, as Sanders“ has shown, the 
Townsend coefficient for the initial electron 
ionization is greatly increased as the voltage in- 
creases, the exponential potential curve men- 
tioned above is much steeper. Hence a sufficient 
gradient for positiye ion ionization is created as 
soon as a portion of the positive ion phalanx is 
swept into the cathode. Thus sparkover ensues 
after a lag corresponding to time for a positive 
ion to cross a fraction of the gap length. As the 
overvoltage is increased this fraction becomes 
smaller. 


Domain three 

For an adequate initial ionization and a higher 
overvoltage on the same gap (or the same percent 
overvoltage in a longer gap”) another effect 
intervenes. Owing to the intense electron ioniza- 
tion at increased voltage the exponential curve 
mentioned is itself steep enough near the anode 
so that positive ion ionization commences there 
at once. (If this secondary ionization is extremely 
close to the anode it might not be sufficiently 
cumulative until this active rear portion of the 
positive ion phalanx has moved an adequate 
distance away from the anode.) For still higher 
oveivoltages (e.g., 60 percent), however, the po- 
tential curve caused by the passage of the 
initial electron avalanche is steep enough to 
give positive ion ionization at some point in the 
middle of the gap at once before the positive 
ions have moved any appreciable distance. The 
mammoth flood of electrons released by this 
positive ion ionization starts for the anode before 
the original avalanche has reached there and 
arrives, starting the voltage drop, one step 
behind the initial avalanche, as it were. Thus the 
time lag is of the order of time for an electron to 
cross the gap. The above explanations for domains 
two and three are worded, for simplicity, as if a 
single electron avalanche sufficed for spark- 


4 Sanders, Phys. Rev. 41, 667 (1932). 


1022 ABE TILLES 


over. It may well be that the space and time co- 
incidence of several is essential under most con- 
ditions; this does not, however, alter the general 
argument. (Because of such coincidence a marked 
increase in initial ionization may tend slightly 
to shift the sparkover from the second domain 
to the third but not to the same degree as an 
increase in voltage.) 
CORRELATION OF THE LITERATURE 

It is particularly gratifying that in the light 
of the above explanation the various available 
experimental reports are seen to be consistent 
with each other and with the spark mechanism 
envisaged, whereas hitherto they had appeared 
isolated and even contradictory. The more im- 
portant of these are briefly presented. 

Zuber* and von Laue,’ working at low initial 
ionization showed that the spark lags were very 
long (running into seconds) and of random dis- 
tribution and that they decreased as the ionizing 
illumination was increased. They established the 
nature of, but not the determining conditions 
for, the first domain. Strigel,®: ® using fairly high 
ionization and fairly high overvoltages, also 
found a random distribution of spark lags at 
much smaller times, illustrating that for test 
conditions giving a short formative time domain 
one extended into quite short intervals. Strigel 
established that the magnitude of the random 
time lags varied greatly with the composition 
and condition of the illuminated cathode. Strigel 
also observed, but was at a loss to interpret, 
the transition from the random distribution to a 
specific, formative time lag. 

Because no instrument existed combining the 
necessary swiftness with sufficient precision as to 
voltage to allow measurements of the formative 
lags of sparkover at low overvoltages, adequate 
measurements in the second domain have not 
hitherto been obtained. Perhaps the only pre- 
vious experimental work bearing on this matter 
is Reukema’s" study of the shift of a.c. sparking 
voltage at high frequency, which gives strong 
indirect verification of the necessity of positive 
ion motion for completion of sparkover at small 
overvoltages. 

Among others, Buss and Masch'® using the 
cathode-ray oscillograph and Beams‘ using the 
Kerr cell, working at high overvoltages, found 


18 Reukema, A. I. E. E. Trans. 47, 38 (1927). 


that the spark lag was very brief, e.g., 3 10-° 
second, and had the definite value of the third 
domain. Viehmann" used the oscillograph in 
conjunction with a saw tooth wave shape of 
voltage whose slope was changed from test to 
test. This, together with the small number of 
sparkovers per series, entirely obscured the differ- 
ence in kind between the various lags he ob- 
served. He did, however, observe lags in all 
three domains and demonstrated roughly the 
tremendous range of spark lags existing and 
their very great shortening as the ionizing illumi- 
nation or overvoltage is increased. 

Observations of the initiation of the visible 
spark directly establish that spark initiation 
may occur by one of two distinct mechanisms. 
This was first observed, in a qualitative way by 
Torok.” Dunnington’s" precise Kerr cell obser- 
vations at various pressures and gap lengths and 
at small overvoltages, supplemented by the 
similar work of von Hamos" at high over- 
voltages, establish two types of sparkover. First, 
one in which the visible spark initiates exclusively 
at the cathode. Second, one in which the spark 
commences somewhere in mid gap. It is further 
established that this second sparkover type is 
caused to replace the first, either by increasing 
the gas pressure or the gap length while using 
the same percent overvoltage, or by increas- 
ing the percent overvoltage. It is to be noted 
that any one of these three increases corresponds 
merely to an increased magnitude of electron 
ionization. Thus for the test conditions of domain 
two, whose time lag indicates a travel of positive 
ions to the cathode, the spark is seen, literally, 
to initiate at the cathode. For the conditions of 
domain three at which more intensive electron 
ionization must occur, the time lag does not 
allow for positive ion migration but calls for 
positive ion ionization at once and the spark is 
seen to commence at once in mid-gap. 

In closing, the author wishes to express his 
indebtedness to Professors L. F. Fuller, L. E. 
Reukema and L. B. Loeb for their kind and un- 
failing guidance throughout the course of this 
investigation; to Professor D. D. Davis and Mr. 
L. E. Evans for their invaluable cooperation in 
furnishing equipment; and to his brother, Mr. 
I. H. Tilles, for assistance in the laboratory. 


* yon Hamos, Ann. d. Physik 7, 857 (1930). 
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-8 The Arc Spectrum of Europium 


d HENRY Norris RUSSELL AND ARTHUR S. KING, Carnegie Institution of Washington, Mount Wilson Observatory 
: (Received November 1, 1934) 
o 7 A N extension of the temperature classifica- TABLE II. Lines in Arc Spectrum of Eu. 
of tion' for longer wave-lengths (now in prog- Siti 
i> ress) confirms the structure of the spectrum J ALA. Int. Cl. a { 
a predicted by Albertson.? This resembles that of s 
Mn I except that the multiplicities are greater sp ans 
by two. The lowest level is 4/76s*(8S°). Combina- 
i i wp s 8 
tions with 4f'6s6p("P *P and P) give the great ssp 
blue triplet and strong furnace lines in the red. =f 
Higher and D terms form series going to d 
e of Eu II as limit. | =i} 7040.20 + 
The terms so far identified are given in Table _ ace 
3 4)-4 8226.81 3 IV + 01 
y I. Terms in Eu I. 3 6470.75 + 08 
6685.27 5 il + 01 
iP 24-3 7336.28 2 Il .00 
d Cus 7369.69 2 il — 04 
uration Designation Level uration Designation Level al —3 7746 10 H Il 20 
“Pq $4505.26 ap — ase 24-3 6522.75 + 13 
pod 6549.14 3 itl + 
5271.93 8 Ill 00 
d The two *P and terms of the same configura- 6822.65 + 
Ss tion belong to the *S and 7S limits of Eu II. 3} -3 6898.27 IV 00 
° 35-4 6841.05 4 IV 
n Their mutual repulsion separates them by more 4-4 6914.83 3 1v + 01 
n than the difference of the limits (1669.3). The = aca 
e terms of multiplicity 10 are the first to be de- J ALA. in. ch a» 
f The classified lines are listed in Table II. 24-3 4947.39 3 IV + 01 
34-2 4968.73 3 IV — 10s 
n To save space they are given by multiplets. 3-3 4962.55 5 it + 2 
t The wave-lengths and intensities are from Eder‘ 4 =4 $098.73 : WV + At 
r supplemented to the red of 47400, and in a few 
— 34-2) 4543.17 1 Vv + 16 
other cases, by measures on Mount Wilson 34-3 4538.06 1 Vv — 0 
4513.20 1 Vv + 
plates—the temperature class from determina- 44-4 4616.52 2 Vv o4 
5 tions by King from these plates. The column 
headed Av gives the residual (O—C) of the ob- 4322.56 2 + 
served wave number from the values resulting =i} 3 ill ‘Os 
: from Table I. The average value of Av, regardless pgp? 3-2 3964.48 1 IV 0.00 
34-3 3963.63 1 IV 00 
34-4 3961.15 1 Iv 
1A. S. King, Astrophys. J. 72, 221 (1930). Seas ; 
. 2 W. Albertson, Phys. Rev. 45, 499 (1934). 54-5 4209.09 2 Ill + .10 
| *Compare Miss Moore's Multiplet Table (Princeton 54-6 4202.65 1 ul 0 
1933). 
*Eder, Wien. Ber. 126, [la, 473 (1917). * Close double. A 
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of sign, is + 0.06 corresponding to about +0.02A 
——a testimony to Eder’s measures. 

Rydberg formulae applied to the S° and D° 
terms give 46,058 and 45,803 for the limit *S° 
of Eu IT. As the denominator m* usually increases 
by slightly more than unity in similar series, the 
value 45,750 has been adopted. This gives for 
the components of highest / the values of n* 


eS 2.542 1.549 f¥D 3.129 2°P 1.908 
3.582 2.600 g®D 4.147 
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which run nearly parallel with the corresponding 
values for Mn I. The terms f*S and f*D appear 
to belong to the 7S limit. The principal ionization 
potential is 5.64 volts. The observed "D terms 
come from 6d and 7d electrons. The 5d term 
should be at about 20,000 (referred to a®S°) and 
its combinations with the P terms are far in the 
infrared. Many strong lines, including some of 
Class I, are still unclassified, and may arise from 
other électron configurations. 
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Proceedings of the Metropolitan Section 


of the 
American Physical Society 


MEETING OF OCTOBER 26, 1934 


HE first meeting of the Metropolitan Section of the American Physical 
Society for the season 1934-1935 was held on Friday, October 26, 1934, 

in Room 301 of the Physical Laboratory of Columbia University in New York 
City. An afternoon and an evening session were held. The presiding officer was 


Professor H. C. Urey. 


The afternoon session was devoted for the most part to a symposium on the 
physics of the solid state; the papers were as follows: 


Free Electrons as the Source of Metallic Binding Energy, F. Seitz AND 
E. WIGNER, Princeton University. 


The Internal Dissipation of Solids, R. L. WEGEL, Bell Telephone Labora- 


tortes. 


Measurement of the Elastic Constants of Crystals at Low Temperatures, 
S. L. Qutmsy, Columbia University. 


Some New Magnetic Properties of Alloys, R. M. Bozortu, Bell Telephone 


Laboratories. 


In the latter part of this session Professor Ladenburg of Princeton University 
gave a brief account of the sessions devoted to nuclear physics at the conference 
of the International Union of Pure and Applied Physics recently held at 


London. 


The evening session was devoted to a lecture by Dr. J. A. Becker of the Bell 
Telephone Laboratories. His subject was “The Life History of Adsorbed 


Atoms and Ions.” 


W. S. Gorton 
Secretary-Treasurer 
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LETTERS TO THE EDITOR 


Prompt publication of brief reports of important 
discoveries in physics may be secured by addressing 
them to this department. Closing dates for this 
department are, for the first issue of the month, the 


On the Velocity of Light 


The determination of the velocity of light with the aid 
of mirrors furnishes some average between the velocities to 
and from the most remote mirror. This average may, of 
course, be a complicated function depending on various 
factors, although it seems universally assumed that the 
velocity either to or from the remote mirror is exactly equal 
to the average velocity. The recent anomalies observed by 
Pease and Pearson' in their experiments in California 
emphasize the need of a closer examination into these 
factors, among which is the structure of the light émployed 
for such measurements. 

All previous experiments have made use of light in 
which the incident ray and the reflected ray were essen- 
tially of the same kind, the difference consisting in the 
change in phase introduced on reflection. Because of this 
close relationship between the rays, I shall call them 
“related.” Here, I wish to urge an experiment in which 
the rays are as unrelated in nature as possible and thus it 
is probable that the velocity obtained can be set equal to 
that of an unidirectional beam. The aim is to compare in 
a single measurement the average velocity furnished by 
such “unrelated” rays with the velocity of “related” rays. 
Should these velocities be found unequal, a greater dis- 
crimination in the definition of the velocity of light will be 
necessary. The possible repercussion on the theory of 
relativity which identifies the velocity of light with the 
fundamental? velocity of the Lorentz transformations 
would make the following experiment especially worth 
while. 

A beam of light in the shape of a slit will pass through 
the apparatus (such as that of Foucault-Michelson) in the 
usual way but instead of being reflected at the distant 
mirror throughout its entire length, the image of the slit is 
to be reflected in one-half .of its length at the mirror and 
the other half is to excite fluorescence of another color in 
an adjoining vessel. The comparison of the velocities 
consists in observing in what measure the fluorescent color 
is shifted with respect to the incident color in the return 
image. 

Comments on the experiment: The time-lag in the 
fluorescence of numerous substances is so short that the 
image radiated back by the fluorescing medium will be 
practically as sharp as the one reflected by the mirror. 
The single bi-convex lens employed by Michelson in his 
early work may, for example, be replaced by two plane- 
convex lens so that one plane-convex lens can be situated 
close to the fluorescing material. A resonance lamp with 
its high efficiency may serve as the vessel adjoining the 
remote mirror. However, the incident and reflected rays 
under these conditions may not be as completely “‘inde- 


twentieth of the preceding month; for the second 
issue, the fifth of the month. The Board of Editors 
does not hold itself responsible for the opinions ex- 
pressed by the correspondents. 


pendent” of each other as when the rays are of different 
wave-lengths. The time-lag of the fluorescence and the 
depth of penetration into the vessel, if the fluorescing 
medium be a transparent substance rather than a fluoresc- 
ing surface, must of course be taken into account. Only 
such fluorescence need be used whose time-lag has been 
found (with the aid of the Abraham-Lemoine shutter for 
example) to be comparable with the time required by the 
mirror to rotate through the smallest detectable angular 
displacement and also it is to be very short as compared 
with the duration of the journey of the light. However, 
the latter errors may be minimized by carrying out the 
experiment with the fluorescing maferial at various dis- 
tances from the rotating mirror. 

The methods now available for constructing efficient 
fluorescing lamps would seem to make this experiment 
feasible without much further technical improvement. On 
the other hand, should the experiment confirm the genera! 
belief in the equality of the velocities, an experiment of 
this kind could determine the time-lag of various radiation 
processes. 

FREED 

Department of Chemistry, 

University of Chicago, 
October 19, 1934. 


' To be publi ned bp the ournal. 
The Mathematic 19, 1923. 


Remarks on the Band Spectrum of Sulfur and the 
Statistics of the Sulfur Nucleus 

A few yearssago Naudé and Christy' examined the fine 
structure of five of the bands of diatomic sulfur in the near 
ultraviolet and made an analysis of them from which they 
calculated that in the normal state of the sulfur molecule 
the internuclear distance is 1.603A. If this result were 
correct it would appear that the sulfur molecule were very 
abnormal! since a distance somewhat greater than 1.8A is 
predicted by several rules, including the covalent bond 
radii and an empirical rule of the author* by which inter- 
nuclear distances can be calculated from vibrational 
frequencies. 

Since no other deviation from these rules of this mag- 
nitude is known, it occurred to the writer to examine the 
data to determine whether the analysis of Naudé and 
Christy were unique. Hé soon came to the conclusion that 
it is not and, for several reasons, is probably incorrect. 
He noticed, for example, that their 4;F’ values show 
certain peculiarities. They appear to be identical within 
experimental! error for the bands 8-1 and 9-1, and also 
for 10-1* up to the point where this band runs into the 
head of 8-0. Now wheh the Naudé-Christy 4:7" values 
for 10-1 and 10-2 are plotted against K they fall on straight 
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lines up to the point where these bands run into the heads 
of 8-0 and 8-1, respectively. The succeeding points deviate 
from these lines and follow lines with a considerably 
greater slope. This might be the result of a perturbation as 
these authors believed, but since it might also indicate that 
they failed to follow the branches of 10-1 and 10-2 through 
the regions where there is a great overlapping of lines, the 
data on these regions were carefully examined. It was 
found that the measured lines can be reassigned to the 
various branches in a very plausible way so tliat the devi- 
ation just mentioned vanishes and we obtain a set of com- 
binations which agree for the bands 8-1, 9-1 and 10-1, 
and are slightly different for 8-0 and 10-2. This, together 
with the fact that the rotational constants obtained from 
these combinations are what we should expect for the 
ground level, is rather conclusive evidence that these com- 
binations are A,F” rather than A;F’ as Naudé and Christy 
believed. 

A new set of combinations was then made to find the 
true 4;F’ values. From these one finds that the rotational 
constants vary in the vibrational states of the upper elec- 
tronic level im such a fashion that it is difficult to make an 
extrapolation to the vibrationless state which has great 
significance. Furthermore some perturbations are found in 
the level v’=8. These duplicate in the bands 8-0 and 8-1 
and are additional confirmation of the correctness of the 
new combinations, Large perturbations are to be expected 
in the upper state since it has been observed‘ that many of 
the band heads are considerably displaced from their 
expected positions. 

The results of the new analysis are given in Table I. An 
important feature is that it shows that the even rotational 


TABLE I. Constants of the sulfur molecule. 


B(em™) D,{cale), ae, 
v’=0 0.309 Df! = -2.2X107 re’ =1.840A 
v’=1 0.307 =0.0016 
v’=2 0.306 Be =0.30%. 
v= 8 0.230 Df = -2.4X107 

(0.211 for K >21) 
v= 9 0.229 
v=10 0.239 


levels are present in the ground state of the sulfur molecule 
instead of the odd levels as found in oxygen. Since the 
electronic states of these two molecules are presumably of 
similar character, this shows that if the oxygen molecule is 
symmetric in the nuclei, as has been generally accepted, 
the sulfur molecule is antisymmetric. However, in case the 
character of the electronic states has been incorrectly 
assigned, it is possible that the symmetry character should 
be reversed in both cases. 
RicHARD M. BapDGER 
Gates Chemical Laboratory, 
California Institute of Technology, 
Pasadena, California, 
October 29, 1934. 


1S. M. Naudé and A. Christy, Phys. Rev. 37, 490 (1931). 

2R. M. Badger, J. Chem. Phys. 2, 128 (1934). 

3 The v’ numbering of Naudé and Christy should be increased by one 
oye to the results of Fowler and Vaidya, Proc. Roy. Soc. 4132, 310 

*A. Christy and S. M. Naudé, Phys. Rev. 37, 903 (1931). 
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Concerning the Recovery Time of Geiger-Mueller Counters 


We have equipped a Braun tube with a timing “sweep” 
which is initiated by the phenomenon to be studied. Thus, 
pulses recurring at random will give retracing curves on the 
fluorescent screen and so may be photographed. Figs. 1-3 
show three oscillograms with a schematic diagram at the 
right of each. 

The ordinates of the curves represent the negative 
potential v of the wire of a Geiger-Mueller counter con- 
nected in the usual way, namely, with the outer cylinder 
at constant (negative) potential V, and the wire grounded 
through a high resistance R. Figs. 1 and 2 show that when 
a discharge occurs, v rises suddenly and then falls following 
an exponential curve. As one would expect its time constant 
is equal to RC where C is the total capacity between the 
wire and ground. Curves 1 and 2 in Fig. 1 show discharges 
for two values V; and V; of the applied voltage. The 
horizontal lines V;— Vo and V2— Vo show the excess of V 
over the starting potential Vo. In other words V—V»o 
represents that value of v for which the voltage across the 
counter equals the starting voltage. 

By suddenly applying x-radiation we have seen that a 
counter will discharge after v has fallen below V—Vo 
but not before. Therefore, the recovery times for the re- 
spective systems are defined by the associated time inter- 
vals denoted by the r’s in the oscillograms. 

In Fig. 1 V;is greater than Vj, so one sees that the excess 
of the initial value of yover V— Vo increases with increasing 
V. Therefore the time +r is larger for higher operating 
voltages. 

The two curves in Fig. 2 were taken with the same 
applied voltage V. In curve 6 the value of C was about 
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twice that in curve a. This increase of C has increased the 
time constant of the curve but decreased the initial ordi- 
nate. The former change tends to increase r while the latter 
tends to diminish it. In this particular case, doubling C 
has increased r by about 27 percent. It seems quite possible 
that in some cases increasing C may actually decrease r. 

If a pulse occurs before the v of a preceding one has fallen 
to zero, it will be smaller than the preceding pulse and 
have a smaller r. So for very rapid counting rates one has 
a statistical distribution of r’s. 

For higher values of V or C, the discharges may take the 
form shown in Fig. 3. This is a photograph of a single 
sweep during which two discharges occurred. (Some 60 
cycle induction is present.) Here rises to V— Vo» and 
remains there for a time T before decaying exponentially. 
The time + is of course given by T and may vary widely 
im successive pulses. The average T is an increasing function 
of C and decreases as R increases. 

The observations described here are quite uniform 
throughout a set of counters made in the same “batch.” 
The wire in all counters was 3 mil tungsten and the cyl- 
inders were oxidized copper. The pressure in all cases was 
7 or 8 cm. The qualitative phenomena are independent of 
counter size, of whether the gas is air or argon and oxygen, 
and of change in R between 5 X 10° and 5 X 10° ohms. 

We hope to submit a report of quantitative nature in 
the near future. 

It is a pleasure to acknowledge the helpful support of 
Dr. W. F. G. Swann. 

W. E. Danrorta 

The Bartol Research Foundation 

of the Franklin Institute, 
Swarthmore, Pennsylvania, 
November 6, 1934. 


On the Statistical Theory of Errors 


Professor R. A. Fisher has most kindly responded to our 
request for criticisms of the article that appeared under the 
above title.! The material in his letter is much too valuable 
to be filed away, so with his consent we here present the 
substance of his comments, together with some additions 
here and there of our own. 

It is doubtful if on page 135 it was made sufficiently clear 
that in the absence of a reliable estimate of o, the u test 
cannot be used, and that the z test (which is equivalent to 
Fisher and Student's ¢ test) is the only recourse. (By a 
reliable estimate of « we mean an estimate that is consider- 
ably more reliable than can be obtained from the single 
sample under test.) The z test is not inherently misleading; 
it tests objectively a proposed value of z, and for this 
purpose it is of course perfectly valid (as we say). Like any 
statistical test, the s test lays down and accepts a perfectly 
definite hazard. Misinterpretations of the z test may be 
common, but the blame should be placed, not on the test 
itself, but on misunderstandings of the nature that we 
point out on page 135. What is more to be feared than 
over-confidence in the z test is the use of the normal prob- 
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ability integral (the u test) with an estimate of « based 
on the single sample under test. 

The separation of the parameters of the parent popula- 
tion from estimates of these parameters has been a gradual 
process. Many writers have been extremely careless in 
confusing that which is estimated with an estimate of it. 
Thinking to avoid any such looseness, we systematically 
used Greek and Latin letters to distinguish the two classes 
of quantity. It is perhaps well to go even further and use 
distinguishing names for the two classes. For this purpose 
there are in use today the terms “parameter” of the parent 
population and “statistic” of the sample, the work 
“statistic” having been introduced by Fisher (footnote 4 
of our article) in 1921 to fill the need of a term antithetical 
to “‘parameter."’ A parent population is completely specified 
by its one or two or more parameters, but a sample of n 
would require n different statistics for its specification. To 
each of these statistics there corresponds a particular 
parameter or parametric function toward which the value 
of the statistic tends as the sample is indefinitely increased; 
but to each parameter there “corresponds,” in this sense, 
as many statistics as there can be of samples from a given 
parent population, to which number there is no limit. 
For these reasons it would doubtless have been better to 
have written “corresponding statistic of a sample” on page 
142, 7 lines below section (3e), to avoid giving the im- 
pression that there is a one to one correspondence between 
the two quantities s and e. 

In further connection with fiducial probability it should 
be mentioned that fiducial values can be taken only from 
distributions of statistics that contain the whole of the 
information that can be ebtained from the sample. The 
distribution of s fulfills this requirement, and our discussion 
of fiducially related values of « and s is therefore valid, but 
it is worth while to note that the distribution of, for 
example, the arithmetic mean deviation, from which 
Peters’ formula (see any text on least squares) is derived, 
could not be so used. There is not room here, and neither 
was there in the original article, to discuss the criteria of 
“efficiency” and “sufficiency,” but they might at least be 
mentioned with a reference. The reader will find them 
discussed in the papers cited in footnotes 4 and 31. 

W. Epwarps DEMING, 
Bureau of Chemistry and Soils, 
Washington, D. C. 
RayMonp T. BirGe, 
University of California 
Berkeley. 

November 9, 1934. 

1 Deming and Birge, Rev. Mod. Phys. 6, 119-161 (1934). 


A Low-Power Positive-Ion Source of High Intensity 


Numerous tests have been made in this laboratory of 
various low voltage and high voltage ion-sources, all of 
which presented serious limitations in use. Dr. F. L. 
Mohler of the National Bureau of Standards kindly de- 
scribed to us last summer some experiments having a 
similar objective, in which with low power he had obtained 
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very large current-densities to a negative probe inserted 
from one side into a hot cathode, low voltage arc confined 
in a quarts capillary. We duplicated this result, but rapid 
deterioration of the quartz and the difficulty of its re- 
placement with the necessary probe-canal accurately in 
alignment led us to consider other materials for the arc- 
capillary. Since every positive ion diffusing to the quartz 
wall is neutralized on contact it appeared that a metal 
capillary, allowed to “‘float’’ in potential, should serve 
equally well. This was promptly confirmed, and a series of 
subsequent experiments, concerned chiefly with the em- 
pirical space-charge problem of focussing such high ion- 
densities into a beam, led to the ion-source shown in Fig. 1. 
Difficulty has been encountered in brief efforts to focus 
these intense ion-beams on the target of the high voltage 
tube, but description of the ion-source may be of immediate 
value to other workers. 

The important features of the ion-source are shown to 
scale in Fig. 1. The steel arc-body A is bored out, forming 
the arc-capillary B joining the arc-spaces which contain 
the oxide-coated filament C and anode D (110 volts with 
resistance). A side volume E for ions diffusing to the probe 
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Fic. 1, Ion-source detail. (Department of 


Magnetism, Car- 


Terrestrial 
negie Institution of Washington, No. 3274.) 
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F is provided with the important ‘‘focussing-diaphragm"’ 
G which gives very high current-density to the small canal 
Hi in the center of the probe. Conical insert J provides flat 
gas-seal faces. The probe is supported on the Lavite in- 
sulator K, and is held at 500 to 7000 volts filtered direct 
current by a Variac-controlled '66 rectifier-set. Focussing 
electrode L is held at 2 to 15 kilovolts direct current by a 
ZP-85 kenotron set with condensers. Protective resistors 
of high value (50,000) are essential to prevent probe- and 
focus-voltages from “blowing out"’ the arc (possibly by 
oscillations). A self-striking “auxiliary arc” (0.1 ampere) 
to the steel arc-body is usually used to strike the main arc 
through the capillary (0.2 to 2 amperes); otherwise the 
application of a small spark-coil to the (floating) arc-body 
starts the capillary-arc. Pressures are roughly 2 in 
the arc and 2 10~* in the tube; speed of Apiezon pumping 
system is about 30 liters/sec. (air); hydrogen flow is 
perhaps 10 cc per hour at N P T, allowing the use of 
deuterium gas without recovery at nominal cost. The ratio 
of atomic to molecular ions is high and is unimportant 
since magnetic analysis is used at the target. 

With a total power-consumption at maximum under 
250 watts, and this gas flow which can be handled by 
modest pumping arrangements, true ion-currents (no 
secondary electrons) up to 1.5 milliamperes are readily 
obtained through the probe-canal 1 mm in diameter and 4 
mm long shown in Fig. 1 (probe-voltage 7 kilovolts; probe- 
current 3 milliamperes; arc 2 amperes). Higher ion-currents 
undoubtedly can be obtained by using a larger probe-canal 
and higher current-density in the capillary. Application of 
6 kilovolts to the focussing electrodes shown (we have 
made tests on various electrode-arrangements) will com- 
pletely focus the probe-canal output-current of 25 micro- 
amperes through a 5-mm diameter hole at a distance of 
11 cm (position of next electrode-gap) when the probe is 
held at 700 volts; 9 kilovolts are required with the probe 
at 1000 volts (80 microamperes), and 15 kilovolts are 
required to focus the probe-output of 250 microamperes at 
1500 volts. The simple “‘lens’’ FL does not bring the whole 
ion-beam into good focus at large distances without proper 
intermediate focussing-gaps. 

M. A. TvuveE 
O. 
C. M. Van Arta! 
Department of Terrestrial Magnetism, 
Carnegié Institution of Washington, 
Washington, D. C., 
November 14, 1934. 
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